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detected  power  at  double-  and  sum-frequencies  is  found  to  be  multiplied  by 
a  coefficient,  which  is  less  than  or  equal  to  unity,  and  which  depends  on  the 
incident  photon  energy  and  on  the  effective  temperature  of  the  system. 

The  Neyman  Type  A  and  Thomas  counting  distributions  provide  a  good 
description  for  the  counting  of  photons  generated  by  multiplied  Poisson 
processes,  if  the  time  course  of  the  multiplication  is  short  in  comparison 
with  the  counting  time.  We  have  obtained  analytic  expressions  for  the 
probability  distributions,  moment  generating  functions,  moments  and  variance- 
to-mean  ratios.  Sums  of  Neyman  Type-A  and  Thomas  random  variables  have  been 
shown  to  retain  their  form  under  the  constraint  of  constant  multiplication 
parameter.  The  doubly  stochastic  Neyman  Type-A,  Thomas  and  fixed  multiplica¬ 
tive  Poisson  distributions  have  also  been  considered.  Applications  include 
the  photon-counting  scintillation  detection  of  nuclear  particles  when  the 
particle  flux  is  low,  the  photon-counting  detection  of  weak  optical  signals 
in  ionizing  radiation,  and  the  design  of  star-scanner  spacecraft  guidance 
systems.  Conditions  under  which  the  Neyman  Type-A  and  Thomas  converge  in 
distribution  to  the  fixed  multiplied  Poisson  and  to  the  Gaussian  have  been 
determined.  This  latter  result  allows  approximate  solution  of  the  related 
detection  and  estimation  problems.  A  more  exact  characterization  of  these 
systems  is  provided  by  a  new  generalized  performance  parameter  for  systems 
with  arbitrary  counting  statistics.  It  was  shown  explicitly  that  this 
parameter  varies  monotonically  with  system  performance  whereas  the  signal-to- 
noise  ratio  does  not.  The  analysis,  which  utilizes  normalizing  transforms 
is  also  useful  in  the  direct  calculation  of  the  Bayes  risk,  probability  of 
error,  and  system  sensitivity.  A  technique  for  evaluating  system  performance 
directly  from  photocounting  cumulants  was  also  derived.  It  is  most  useful 
in  situations  where  exact,  closed  form,  explicit  solutions  exist  for  the 
cumulants,  but  not  for  the  distribution,  such  as  in  the  general  case  of 
communication  through  atmospheric  turbulence. 

General  expressions  are  obtained  for  the  mean  and  variance  of  the 
number  of  events  in  a  fixed  but  arbitrary  sampling  time  for  a  nonparalyzable 
dead-time  counter.  The  input  is  assumed  to  be  a  Poisson  point  process  whose 
rate  is  a  stochastic  process,  and  the  dead  time  is  assumed  to  be  small  in 
comparison  with  the  fluctuation  time  of  the  driving  process.  The  mean  is 
shown  to  depend  only  on  the  first-order  statistics  of  the  rate,  whereas  the 
variance  is  formally  shown  to  depend  on  both  the  first-  and  the  second-order 
statistics  of  the  rate.  For  the  particular  process  arising  in  the  detection 
of  chaotic  light,  an  explicit  expression  is  obtained  for  the  dependence  of 
the  dead-time-modified  variance  on  the  power  spectrum  of  the  radiation.  It 
is  demonstrated  that  the  variance  takes  on  a  particularly  simple  form  for 
chaotic  light  with  Lorentzian  and  Gaussian  spectra.  In  the  region  in  which 
our  study  is  valid,  it  turns  out  that  the  dead-time  dependence  of  the  varianc 
is  contained  in  a  multiplicative  function  that  is  essentially  independent  of 
the  spectral  properties  of  the  light. 

Current  suppression  in  AIq. 35Gao. 65As-GaAs  N-p  heterojunctions  is 
observed  experimentally.  The  data  are  explained  by  examining  the  carrier 
transport  modified  by  electron  accumulation  in  the  interface  potential  notch 
resulting  from  the  conduction-band  discontinuity.  In  addition,  it  is  found 
that  the  temperature  dependent  I-V  measurements  can  be  used  to  deduce  the 
magnitude  of  AEC  at  the  junction  interface. 
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A  new  method  has  been  developed  to  form  a  p  -n  polycrystalline- 
crystalline  junction  at  a  low  temperature.  The  technique  involves  the  use 
of  an  E-beam  evaporated  polysilicon-aluminum-polysilicon  multi-layer 
structure  to  obtain  an  Al-rich  polycrystalline  film  by  heat  treatment  at 
60noc.  Experimental  evidence  indicates  that  the  current  conduction  is 
dominated  by  electron  recombination  at  the  interface  between  the  polvsilicon 
and  the  substrate.  The  best  fit  between  a  simple  model  and  experimental 
data  yields  an  effective  recombination  velocity  of  5  x  10^  cm/sec,  a  value 
consistent  with  carrier  recombination  at  the  grain  boundaries  of  poiv- 
crvstalline  silicon. 

A  Schottky  barrier  is  fabricated  with  a  silicon-aluminum-silicon 
structure  by  the  irradiation  of  either  a  pulsed  or  cw  laser.  The  completed 
diodes  have  a  large  barrier  height  (1  eV)  and  good  rectifying  properties. 

Thermoluminescent  (the  glow  curve)  and  phosphorescent  decay  experiments 
have  been  performed  for  a  range  of  manganese  activated  zinc  silicate 
phosphors.  It  is  observed  that  the  phosphorescent  decay  time  is  a  non¬ 
monotonic  function  of  temperature  when  the  light  output  is  contributed 
simultaneously  by  two  energetically  adjacent  traps.  By  representing  the 
decay  curve  with  exponential  terms,  the  calculated  trap  depths  can  be 
correlated  with  the  data  obtained  from  the  glow-curve  experiment.  Thus, 
the  apparent  contradiction  between  the  relative  insensitivity  of  the  decay 
time  with  respect  to  temperature  and  the  thermal  activation  of  traps  can 
be  reconciled.  In  isolated  craps,  the  decay  time  is  always  a  monotonic 
function  of  temperature  although  the  decay  curve  may  not  be  exponential. 

The  dependence  on  the  energy  and  electric  field  of  Si-Si02  interface 
state  capture  cross  sections  has  been  measured  for  Al-Si02-nSi  structures 
with  an  oxide  thickness  -  800  &.  The  capture  cross  section  for  holes. 

Op,  has  been  observed  to  be  several  orders  of  magnitude  larger  than  that 
for  electron,  an,  and  further  it  has  been  observed  to  show  a  strong 
dependence  on  energy  in  the  upper  half  of  the  silicon  energy  gap.  This 
apparent  depenaence  of  a p  on  energy  has  been  shown  to  arise  primarily  from 
a  dependence  on  the  electric  field  at  the  semiconductor  surface.  For  the 
same  electric  field,  interface  states  at  different  energies  in  the  gap  are 
observed  to  have  very  nearly  the  same  magnitude  of  capture  cross  section. 

The  minimum  dimensions  of  bipolar  transistors  with  ion  implanted 
impurity  profiles  for  predictable  device  and  system  behavior  is  calculated 
on  the  basis  of  limitations  arising  from  the  random  positions  of  the 
implanted  impurities.  Doubly-stochastic  effects  arising  from  uncertainty 
in  exact  implantation  parameters  are  considered.  For  a  chip  yield  of  99%, 
the  minimum  emitter  size  is  of  the  order  of  1  £  L  £  5  pm  when  a  straight¬ 
forward  scaling  theory  is  applied  to  dimensions  and  impurity  concentrations 
from  present  device  designs  and  a  parity  check  on  a  bit  word  is  applied. 

We  have  determined  the  mechanism  and  rate  for  energy  transfer  between 

the  v„  and  v,  modes  of  C0F„.  A  novel  relaxation  behavior  has  been  observed 
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for  this  process.  Temperature  dependent  relaxation  rates  have  been 
determined  for  vibrational  energy  transfer  in  SO2.  The  overall  translational- 
vibrational  energy  transfer  process  in  oxalyl  fluoride  has  been  observed. 

This  molecule,  which  has  a  pseudo-quasi-continuum  of  levels  due  to  internal 
rotation,  is  a  model  for  larger  molecules  which  undergo  multiquantum 
infrared  absorption. 

Collision  assisted  laser  pumping  of  SFg  has  produced  molecular  beams 
with  internal  energies  exceeding  3500  cm~l/molecule .  Such  beams  are 
particularly  useful  in  studies  of  surface  chemical  reactions  and  crossed 
molecular  beam  reactions. 

An  excimer  laser  has  been  used  to  produce  high  velocity  Cl  atoms  via 
photodissociation.  Upon  collision  with  virtually  any  small  molecule,  these 
fast  atom  fragments  cause  vibrational  excitation  of  the  molecule. 

We  have  discovered  several  novel  features  of  the  dynamic  behavior  of 
Br2  molecules  trapped  in  argon,  krypton,  and  xenon  matrices  at  10-30°K 
when  irradiated  by  a  pulsed  dye  laser.  Infrared  emission  from  several 
electronic  states  is  observed  and  evidence  for  Br2~Xe  molecular  exciplex 
species  has  been  found. 

129 

We  have  successfully  polarized  Xe  by  spin  exchange  with  optically 
pumped  rubidium.  The  polarization  was  detected  by  adiabatic  reversal  of 
the  xenon  spins,  as  described  in  the  preceding  progress  report. 

We  have  discovered  that  the  cesium  atoms  within  our  glass  sample  cells 
lase  continuously  on  many  transitions  in  the  infrared  spectral  region, 
e.g.  7P  to  6S  at  3.095p  and  2.93y;  7S  to  6P  at  1.469u  and  1.359u;  7P  to 
5D  at  1.376u  and  1.360u.  The  atoms  are  excited  with  cw  laser  radiation 
at  4593X  and  4555$  and  the  lasing  occurs  without  any  feedback  mirrors.  We 
believe  this  is  the  first  reported  example  of  mirror less  cw  lasing. 

We  have  completed  our  studies  of  the  newly  discovered  infrared  absorption 
bands  of  dense  alkali  vapors  of  Cs,  Rb,  K  and  Na.  All  elements  show  very 
similar  structure  between  lu  and  2p,  structure  which  may  be  due  to  absorp¬ 
tion  from  triplet  molecules.  Preliminary  attempts  to  observe  these  bands 
in  fluorescence  failed  to  show  any  signal. 

We  have  detected  strong  two  photon  excitation  of  rubidium  according 
to  the  scheme  discussed  in  the  previous  progress  report.  Substantial 
(i.e.  10%  to  20%)  spin  polarization  is  also  detected.  However,  our 
analysis  of  the  data  is  complicated  by  the  occurrence  of  strong  excimer 
absorption  which  also  leads  to  spin  polarization  of  a  comparable  order  of 
magnitude. 

Extensive  measurements  of  the  magnetic  circular  dichroism.of  the 
CsAr  excimer  bands  and  CS2  dimer  bands  have  been  made  during  the  past 
year.  We  have  been  able  to  gain  for  the  first  time  unambiguous  information 
about  the  angular  momentum  and  angular  momentum  coupling  schemes  of  the 
electronic  states  of  excimer  molecules. 
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We  have  used  spin  exchange  magnetic  resonance  spectroscopy  to  unambiguously 
identify  free  hydrogen  atoms,  free  deuterium  atoms  and  free  electrons  in 
cells  containing  mixtures  of  hydrogen  or  deuterium  and  cesium  vapor  which 
is  optically  pumped  by  the  4593$  second  resonance  line  of  the  cesium  atom. 

The  spin  destruction  rates  of  sodium  in  xenon  have  been  measured 
and  found  to  be  even  greater  than  the  rates  of  rubidium  in  xenon. 

Van  der  Waals  molecules,  NaXe,  are  also  found  to  play  an  important  role 
in  the  relaxation. 

The  dynamics  of  photodissociation  of  diphenyl-diazomethane  and  energy 
relaxation  in  the  diphenylcarbene  fragment  are  investigated  using  picosecond 
laser  methods.  The  diphenylcarbene  singlet-ground  triplet  equilibrium  con¬ 
stant,  free  energy  difference,  energy  separation,  and  a  maximum  value  for 
the  carbon-nitrogen  bond  energy  in  diphenyldiazomethane  are  obtained. 

A  method  for  determining  the  dynamics  of  intramolecular  chain  motion 
has  been  developed  which  employs  the  rate  of  exciplex  formation.  The 
kinetics  of  the  end  to  end  motions  of  the  terminal  donor  and  acceptor  mole¬ 
cules  determines  the  exciplex  formation  rate.  The  decay  of  the  excited 
acceptor  and  the  rise  time  of  the  exciplex  emission,  following  picosecond 
laser  excitation  of  the  acceptor,  is  used  to  determine  the  chain  dynamics. 

The  motion  is  found  to  be  exponential  in  nonpolar  solvents  for  the  mole¬ 
cule  A-(CH2)3~DMA  and  is  the  first  experimental  observation  of  the  end  to 
end  dynamics  of  short  chain  molecular  motion. 

Photoexcitation  of  extended  conformers  A-(CH2)3~D  in  the  highly  polar 
medium  acetonitrile  is  followed  by  very  rapid  electron  transfer  yielding 
the  solvated  ion  pair  (7  ps).  The  rapid  rate  of  ion  pair  formation  indi¬ 
cates  that  there  are  no  strong  geometric  requirements  for  this  process 
and  furthermore  that  this  electron  transfer  process  occurs  directly  with¬ 
out  passing  through  an  excited  CT  complex  state.  This  is  unlike  the 
strong  geometric  requirements  for  excited  CT  complex  formation  both  in 
polar  and  nonpolar  media.  These  studies  were  carried  out  using  a  Nd : 
phosphate  glass  picosecond  laser  for  excitation  and  a  streak  camera- 
optical  multichannel  analyzer  for  fluorescence  detection. 

Two-pulse  photon  echo  experiments  have  been  performed  in  Na  and  Li  vapor 
which  probe  the  sub-doppler  regime  of  the  resonance  line.  In  Li  we  have 
made  the  first  observation  of  the  effects  of  velocity  changing  collisions 
in  the  asymptotic  regime  where  the  total  scattering  cross  section  obtains. 

A  method  for  observing  echo  signals  has  been  developed  which  exploits  the 
dependence  of  the  echo  polarization  on  the  magnetic  substates  which  are 
excited.  A  diagrammatic  method  was  developed  to  provide  a  simple  unified 
procedure  for  analyzing  coherent  reordering  effects  in  gases. 

Photon  Echo  Nuclear  Double  Resonance  (PENDOR)  experiments  on  the 
3h4  transition  of  PR^+  in  LaF^  have  unambiguously  determined  the  exci¬ 
ted  state  splitting  parameters.  Improvements  in  experimental  techniques 
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have  allowed  us  to  obtain  echo  modulation  data  for  the  transition 

which  when  analyzed  will  provide  a  rather  complete  description  of  its 
spectroscopic  and  relaxation  character.  Initial  steps  are  being  taken  to 
perform  excited  state  and  tri-level  photon  echo  experiments  in  Pr3+:LaF3. 
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I.  QUANTUM  DETECTION  AND  SENSING  OF  RADIATION 


A.  COHERENT  DETECTION  AND  SENSING  IN  THE  INFRARED 
(M.  C.  Teich,  L.  Edelshteyn) 

(Principal  Investigator:  M.  C.  Teich  (212)  280-3117) 

Heterodyne  detection  has  a  long  and  august  history  in  the  annals  of 
electrical  engineering,  reaching  back  to  the  earliest  years  of  the  century. 

The  term  has  its  roots  in  the  Greek  words  "heteros"  (other)  and  "dynamics" 
(force) . 

In  1902,  Reginald  Fessenden  was  awarded  a  patent ^  "relating  to  cer¬ 
tain  improvements  ...  in  systems  where  the  signal  is  transmitted  by  [radio] 
waves  differing  in  period,  and  to  the  generation  of  beats  by  the  waves  and 
the  employment  of  suitable  receiving  apparatus  responsive  only  to  the  combined 
action  of  waves  corresponding  in  period  to  those  generated  ..."  The  sub¬ 
sequent  realization  that  one  of  these  waves  could  be  locally  generated  (the 
development  of  the  local  oscillator)  provided  a  substantial  improvement  in 
system  performance.  Practical  demonstrations  of  the  usefulness  of  the  tech¬ 
nique  were  carried  out  between  the  "Fessenden  stations"  of  the  U.  S.  Navy 
at  Arlington  (Virginia)  and  the  Scout  Cruiser  Salem,  between  the  Salem  and 
the  Birmingham  (1910),  and  at  the  National  Electric  Signaling  Company.  In 
1913,  John  Hogan  provided  a  thoroughly  enjoyable  accout  of  the  development, 

use,  and  performance  of  the  Fessenden  heterodyne  signaling  system  in  the 

(2) 

first  volume  of  the  Proceedings  of  the  IRE. 

It  was  not  long  thereafter,  in  1917,  that  Edwin  H.  Armstrong  of  the 
Department  of  Electrical  Engineering  at  Columbia  University  carried  out  a 
thorough  investigation  of  the  heterodyne  phenomenon  occurring  in  the  oscillation 

(3) 

state  of  the  "regenerative  electron  relay".  A  major  breakthrough  in  the 

field,  the  development  of  the  superheterodyne  receiver,  was  achieved  by 
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(4) 

Armstrong  in  1921,  and  this  famous  invention  is  now  used  in  systems  as 
diverse  as  household  AM  radio  receivers  and  microwave  Doppler  radars.  The 
prefix  "super"  refers  to  "the  super-audible  frequency  that  could  be  readily 
amplified  .  " 

In  the  succeeding  years,  the  application  of  heterodyne  and  superheterodyne 
principles  followed  the  incessant  march  toward  higher  frequencies  that  cul¬ 
minated  in  the  remarkable  developments  in  microwave  electronics  about  the 
time  of  World  War  II. 

The  marriage  of  heterodyning  and  the  optical  region  took  place  in  1955. 

In  a  now  classic  experiment,  Forrester,  Gudmundsen,  and  Johnson^^  observed 

the  mixing  of  two  Zeeman  components  of  a  visible  (incoherent)  spectral  line 

in  a  specially  constructed  photomultiplier  tube.  With  the  development  of  the 

laser,  optical  heterodyning  became  considerably  easier  to  observe  and  was 

studied  in  1962  by  Javan,  Ballik,  and  Bond^  at  1.15  pm  using  a  He-Ne  laser, 

and  by  McMurtry  and  Siegman^^  at  6943  £  using  a  ruby  laser. 

The  development  of  new  transmitting  and  receiving  components  in  the 

middle  infrared  region  of  the  electromagnetic  spectrum  led  Teich,  Keyes,  and 

(8) 

Kingston^  in  1966  to  perform  a  heterodyne  experiment  using  a  C0?  laser  at 

10.6  pm  in  conjunction  with  a  copper-doped  germanium  photoconductive  detector 

operated  at  4°K.  Subsequent  experiments  with  lead-tin  selenide  photovoltaic 

detectors  confirmed  the  optimal  nature  of  the  detection  process.  ^  (H)  The 

state  of  the  art  in  applying  these  results  to  coherent  infrared  radar  was 

reviewed  by  Kingston  in  1977.^^  Elbaum  and  Teich^^  have  recently  shown 

the  importance  of  properly  choosing  the  performance  measure  for  a  heterodyne 

system.  This  can  sometimes  be  critical  to  avoid  faulty  estimates  of  signal 

and  noise  levels  leading  to  inconsistent  or  incorrect  results,  as  pointed 

(8), (14) 


out  by  a  number  of  authors. 


The  foregoing  introduction  has  dealt 


with  more-or-less  idealized  systems;  it  must  be  kept  in  mind  that  there  are 
a  variety  of  effects  that  can  alter  heterodvne  performance  in  important  ways. 

We  point  out  also  that  a  number  of  system  configurations  employing  dif¬ 
ferent  forms  of  nonlinear  heterodyne  detection  have  been  proposed  for  various 
applications .  ’  (1^)  •yhggg  maw.e  use  of  multiple  frequencies,  nonlinear 

detectors,  and/or  correlation  schemes.  All  are  aimed  at  increasing  the  signal- 
to-noise  ratio  in  situations  for  which  usual  operating  conditions  are  relaxed 
in  particular  ways.  One  of  these  schemes,  using  a  two-frequency 

transmitter,  is  currently  under  consideration  for  possible  application  in 
over-water  pilot  rescue  missions  where  the  use  of  a  light-chopping  mechanism 
is  interdicted  because  of  its  interruptive  nature  (private  communication  with 
W.  Tanaka,  Naval  Weapons  Center). 

/  1  Q  \ 

We  have  carried  out  a  study'  in  which  we  relaxed  the  assumption  that 

heterodyne  detection  cakes  place  by  means  of  photon  absorption,  as  it  does  in 

(19) 

the  middle  infrared  and  optical  regions,  or  by  means  of  electric-field 

detection,  as  it  does  in  the  radiowave  and  microwave  regions.  This  has  en¬ 
abled  us  to  examine  the  heterodyne  detection  process  in  the  submillimeter/ 
far-infrared  region  of  the  spectrum  in  a  simple  way.  Caution  must  be  exercised, 
however,  in  that  the  treatment  we  provide  is  intended  to  provide  only  a  quali¬ 
tative  description  of  the  underlying  detection  process. 

The  increasing  importance  of  heterodyning  and  electric-field  detection  in 
the  submillimeter,  infrared,  and  optical  has  encouraged  us  to  examine 

the  operation  of  such  systems  in  relation  to  the  more  conventional  photon- 

/  o  n  \  /  /  \ 

absorption  detector. 

For  heterodyne  mixing  with  coherent  signals  in  a  two-level  detector,  a 

simple  argument  indicates  that  double-  and  sum-freauencv  terms  ir.  the  detected 

power  are  multiplied  bv  the  factor  sech(hv/2kT  ),  which  varies  smoothlv  from 

e 
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unity  in  the  electric-field  detection  regime  to  zero  in  the  photon-absorption 

detection  regime.  This  factor  depends  on  both  the  incident  photon  energy 

hv  and  on  the  effective  excitation  energy  of  the  (two-level)  detector,  kT^. 

It  is  observed  that  these  terms  only  appear  when  it  cannot  be  determined 

whether  photon  absorption  or  photon  emission  has  taken  place.  Difference- 

frequency  signals,  on  the  other  hand,  arise  from  our  inability  to  determine 

(19) 

from  which  beam  a  photon  is  absorbed  in  a  heterodyne  experiment. 

We  consider  a  simple  hypothetical  two-level  system  with  an  effective 
temperature  T  .  It  is  assumed  that  the  system  responds  linearly  to  the  inci¬ 
dent  radiation  intensity,  and  that  its  interaction  with  the  field  is  suf¬ 
ficiently  weak  such  that  the  state  of  the  field  is  not  perturbed  by  the  pre¬ 
sence  of  the  detector. 

We  label  the  initial  and  final  states  of  the  system  as  jot)  and  and 
of  the  radiation  field  as  |i)  and  |f),  respectively .  For  an  electric-dipole 
transition,  the  transition  probability  Wfi  (which  is  related  to  the  detected 
power)  is  given  approximately  by 

Wf .  =  |<f  n|eq(E+  +  E')|ai>|2,  (1) 

where  e  is  the  electronic  charge,  q  is  the  detector  coordinate,  and  E+  and 
E  are  the  positive-  and  negative-portions  of  the  electric-field  operator, 
respectively.  Since  E+  corresponds  to  photon  absorption  or  annihilation,  and 
E  corresponds  to  photon  emission  or  creation,  the  transition  probability  may 
be  written  as 

Wj^  =  |<fu|eqa  E+|(i^  +  (fileqa^E  jui)]2,  (2) 

where  jf.^  and  |  u^>  represent  the  lower  and  upper  states  of  the  two-level  system, 
respectively.  This  equation  assumes  that  the  atomic  system  is,  in  general,  in 


\ 


a  superposition  state.  Using  microscopic  reversibility,  the  quantity 
2 

|<^u|eqj£)|  ,  which  represents  the  quantum  efficiency  k  ,  may  be  factored  out 

(23) 

of  Eq.  (2).  We  then  sum  over  the  final  states  of  the  field,  which  are  not 
observed,  to  obtain 


W  .  <i|»„**„EVn> 

+ /i|a„*a  EE  +  a„a  *E+E+|i\  (3) 

'  1  £  u  £  u  1  / 


The  normally  ordered  first  term  corresponds  to  stimulated  absorption,  the 

(25) 

antinormally  ordered  second  term  corresponds  to  photon  emission,  and  the 
third  is  an  interference  term. 

We  now  assume  that  before  the  interaction,  the  probability  amplitudes  of 
the  two  possible  states,  a.  and  a  ,  were  related  by  the  Boltzmann  factor,  with 

U 

lower  and  upper  level  energies  represented  by  E.  and  E  ,  respectively,  and 

X  u 

with  excitation  energy  kTg  defining  the  effective  temperature  of  the  detector. 
Thus , 

|au|2/(a£|2  =  exp[-(Eu-Ep/kTe],  (A) 


yielding 


=  (1  +  e-V1/2 


i<t> 


and 


(5) 


a 

u 


«'x/2d 


♦  ,-v1/2 


ie 

e  , 


(6) 


with  x  =  hv/kTg  and  e*®,  e^  representing  phase  factors.  Then,  generalizing 
to  an  arbitrary  radiation  field  represented  by  the  density  operator  p,  we 
obtain  the  expression 
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W  “  (^[eX,/‘"  sech(x/2) ] tr{pE  E+}  +  [e  X^~  sech(x/2)  ]  tr(oE+E  } 

+  [sech(x/2) ] tr{p(e  iYE  E  +  e^E+E+)}j  ,  (7) 

with  Y  =  <t>  -  9 . 

Considering  ideal  heterodyne  detection,  i.e.,  two  parallel,  monochromatic, 

and  coherent  waves  of  frequencies  and  v,  impinging  normally  on  the  detector, 

and  neglecting  spontaneous  emission  so  that  the  antinormally  ordered  and  the 

(25) 

normally  ordered  terms  are  equal  in  magnitude,  the  first  two  terms  above 

generate  dc  and  difference-frequency  signals,  while  the  third  term  contributes 
double-  and  sum-frequency  signals.  This  may  be  clearly  seen  by  explicitly 
rewriting  Eq.  (7)  as 

w  -  |  e°!  2  +  \t°2\2  +  2|e°|  |  I  cos  [2r(v1-v2)  t  +  (6-a)] 

+  [sech(hv/2kTe)  ]  {  |  e°  ]  2  cos(Am;1t-2ci-Y)  +  |  e°  |  2cos  (47rv2t-2S-Y) 

+  2|e°|  |e°|cos[2Tr(v1+v2)t  -  (a+6)-y] }  (8) 

where  =  |e°|e*a  represents  the  complex  electric-field  amplitude  of  the 

constituent  field  with  frequency  and  phase  a.  Double-  and  sum-frequency 

terms  in  the  heterodyne  signal  are  therefore  multiplied  by  the  factor 

sech(hv/ 2kT  ) . 

e 

For  hv/kT^  -*•  0,  this  factor  approaches  1  and  the  classical  electric- 
field  heterodyne  signal  obtains 

Welec  “  n  !  c°s  ( 2ttv;l  t-a  -  +  j  e°  |  cos(2rv2t-B-  ]  2 .  (9) 

For  hv/kT^  -*  °°,  sech(hv/2kTe)  -*■  0  and  the  photon-absorption  (optical) 
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SECH  (hv/2kT) 


heterodyne  signal  obtains, 

Wabs  “  ^  I  I  2  +  \€°2\2  +  2le°!  |e°  |c°s[2ir(v1-v2)t  +  (B-a)]}.  (10) 

A  graphical  presentation  of  the  function  sech(hv/2kT  )vs.  (hv/kT  )  is  provided 

e  e 

in  Fig.  1. 


hv/kT 


Figure  1.  Factor  sech(hv/2kT)  vs.  hv/kT.  This  quantity  appears  as 

a  coefficient  in  the  absorption/emission  interference  term 


The  result  for  the  heterodyne  case  is  easily  reduced  to  the  direct 


detection  (video)  case  for  a  coherent  signal  by  setting  | £ ^ I  =  0,  which  yields 

W  «  |e°|2  {1  +  [sech(hv/2kT  ) ]  [cos (4ttv.  t  -  2a  -  y)]}.  (11) 

dir  1  e  1 

Thus,  double-frequency  intensity  fluctuations  are  discerned  for  electric- 
field  direct  detectors,  while  they  are  suppressed  for  photon-absorption  direct 
detectors  which  respond  simply  as  |e°|2. 

The  photon-absorption  detector  is,  by  definition,  initially  in  its  ground 
state  and  functions  by  the  annihilation  of  a  single  (in  general  norunono- 
chromatic)  photon.  The  presence  of  the  difference-frequency  signal  is  under¬ 
stood  to  arise  from  our  inability  to  determine  from  which  of  the  two  constit- 

(19) 

uent  beams  the  single  photon  is  absorbed.  The  two-level  electric-field 

detector,  on  the  other  hand,  has  equal  probability  of  being  in  the  lower  and 
in  the  upper  state,  so  that  the  pure  processes  of  photon  annihilation  and 
photon  absorption  occur  with  equal  likelihood,  and  we  must  add  the  effects 
of  both.  When  we  are  unable  to  determine  which  of  these  processes  is  occur¬ 
ring,  we  must  add  amplitudes  rather  than  squares  of  amplitudes,  thereby  allow¬ 
ing  interference  to  occur.  Any  attempt,  in  this  case,  to  determine  whether 
photon  emission  or  photon  absorption  takes  place  would  randomize  the  phase 
Y,  and  thereby  wash  out  the  sum-  and  double-frequency  components.  In  general, 
then,  a  photon  incident  on  a  video  detector  induces  upward  and  downward  transi¬ 
tions  with  different  probabilities.  This,  in  turn,  creates  a  quantum-mechanical 
probability  density  (and  charge  distribution)  that  varies  in  time,  producing 
a  current  at  the  incident  radiation  frequency.  The  power  absorbed  then  con¬ 
tains  a  double-frequency  signal.  For  heterodyne  detection,  in  the  general 
case,  sum-frequency  signals  are  observed  as  well. 

The  foregoing  heuristic  model  yields  a  simple  result  for  an  idealized 
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two-level  system.  Replacing  the  operator  eq  by  the  non-relativistic  Hamiltonian 

„  2  /v  ^ 

(p  -  eA)  /2m,  where  p  and  A  represent  the  momentum  and  vector-potential  opera¬ 
tors,  respectively,  would  allow  transitions  more  general  than  electric-dipole , 
and  absorptions  of  more  than  one  quantum  to  occur.  A  rigorous  treatment  should 
consider  a  collection  of  such  systems  (as  a  model  for  a  bulk  photodetector 
or  metal  antenna),  in  the  presence  of  a  surrounding  reservoir,  and  should  be 
carried  out  using  the  density  matrix  formalism.  It  is  expected  that  the  ef¬ 
fects  described  here  are  important  in  a  broad  range  of  systems,  including 

( 26^  ( 27^ 

the  Josephson  detector.  ''  ' 

(  28  1 

Tucker'  '  has  recently  carried  out  a  rigorous  analysis  of  quantum-limited 
detection  in  tunnel  junction  mixers.  He  demonstrates  that  nonlinear  tunneling 
devices  are  predicted  to  undergo  a  transition  from  energy  detectors  to  photon 
counters  at  frequencies  where  the  photon  energy  becomes  comparable  to  the 
voltage  scale  of  the  <lc  nonlinearity.  It  is  apparent  that  the  character  of 
this  result  is  closely  related  to  the  discussion  presented  here. 
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B.  THE  EFFECTS  OF  WAVE  NOISE  ON  THE  ABSORPTION  AND  DETECTION  OF 
INFRARED  AND  OPTICAL  RADIATION* 

(M.  C.  Teich,  P.  R.  Prucnal,  B.  E.  A.  Saleh) 

(Principal  Investigator:  M.  C.  Teich  (212)  280-3117) 

It  is  just  over  AO  years  since  Neyman  introduced  a  family  of  straight¬ 
forward  but  intriguing  generalizations  of  the  binomial  and  Poisson  distribu¬ 
tions.^  He  called  the  simplest  and  perhaps  most  useful  of  these  "the 
contagious  distribution  of  Type-A  with  two  parameters"  (see  Fig.  1).  Neyman's 
work  was  motivated  by  a  variety  of  experimental  observations  in  entomology 
and  bacteriology  with  which  calculations  based  on  the  ordinary  Poisson  failed 
to  agree.  He  reasoned  that  the  effects  of  contagion  were  important  in  these 
studies,  and  he  succeeded  in  introducing  this  property  in  a  remarkably  ele¬ 
gant  yet  simple  way.  The  description  "contagious"  implies  that  each  favor¬ 
able  event  increases  (or  decreases)  the  probability  of  succeeding  favorable 
(2) 

events.  Feller  and  others  have  argued  that  there  are  essentially  two  kinds 

of  contagion:  "true  contagion"  as  described  above,  and  "apparent  contagion", 

where  there  is  an  inhomogeneity  of  the  population.  It  has  long  been  known  that 

certain  distributions,  such  as  the  negative  binomial  and  the  Neyman  Type-A, 

(2) 

could  be  derived  in  terms  of  both  types  of  contagion.  Indeed,  the  negative 

(3) 

binomial  distribution  was  obtained  by  Greenwood  and  Yule  in  terms  of  ap¬ 
parent  contagion  and  subsequently  by  Eggenberger  and  Polya  in  an  independent 
study  in  terms  of  true  contagion. 

A  closely  related  counting  distribution  has  also  proved  useful  in  our 
study.  In  1949,  Marjorie  Thomas^  introduced  a  two-parameter  counting  dis¬ 
tribution  distinct  from  the  Neyman  Type-A  only  in  that  primary  pulses  ap¬ 
peared  in  the  final  process  along  with  secondary  pulses.  Though  she  originally 
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PROBABILITY  DENSITY 


NUMBER  OF  COUNTS 


Figure  1:  Semilogarithmic  plot  of  the  Neyman  Type-A  counting  distribution 

p(njw)  vs  the  count  number  n.  The  overall  mean  (n)  =  10.  Curves 
are  labeled  according  to  the  value  of  the  parameter  a.  In  the 
limit  as  a  4  0,  the  Neyman  Type-A  approaches  the  Poisson  (dashed 
curve).  The  multimodal  character  of  the  distribution  is  evident. 


referred  to  this  as  the  double-Poisson  distribution,  it  has  since  come  to  be 


called  the  Thomas  distribution.  Like  the  Neyman  Type-A,  it  is  obtainable  from 
a  Neyman-Scott  cluster  process. 

We  have  developed  various  properties  and  limits  of  the  simple  and  doubly 

stochastic  versions  of  the  two-parameter  Neyman  Type-A,  Thomas,  and  fixed 

(8) 

multiplicative  Poisson  distribution.  We  have  discovered  many  applications 

to  photon,  particle,  and  pulse  counting  in  optics  and  have  begun  to  investigate 

the  performance  of  systems  containing  such  counting  detectors.  Finally,  we 

have  looked  into  a  number  of  related  contagious  distributions  pertinent  to 

photon  counting  and  optical  communications. 

The  Neyman  Type-A  and  Thomas  counting  distributions  turn  out  to  provide 

a  good  description  for  the  counting  of  photons  generated  by  multiplied  Poisson 

processes,  as  long  as  the  time  course  of  the  multiplication  is  short  in  com- 

(9)  (9) 

parison  with  the  counting  time.  We  have  obtained  analytic  expressions 
for  the  probability  distributions,  moment  generating  functions,  moments,  and 
variance- to-mean  ratios.  Sums  of  Neyman  Type-A  and  Thomas  random  variables 
have  been  shown  to  retain  their  form  under  the  constraint  of  constant  multi¬ 
plication  parameter.  Conditions  under  which  the  Neyman  Type-A  and  Thomas 
converge  in  distribution  to  the  fixed  multiplicative  Poisson  and  to  the  Gaussian 
have  been  determined.  This  latter  result  is  most  important  for  it  provides 
a  ready  solution  to  likelihood-ratio  detection,  estimation,  and  discrimina¬ 
tion  problems  in  the  presence  of  many  kinds  of  signal  and  noise.  The  doubly 
stochastic  Neyman  Type-A,  Thomas,  and  fixed  multiplicative  Poisson  distribu¬ 
tions  have  also  been  considered. 

(9) 

A  number  of  explicit  applications  are  important.  These  include 
(1)  the  photon-counting  scintillation  detection  of  nuclear  particles,  when 
the  particle  flux  is  low,  (2)  the  photon-counting  detection  of  weak  optical 
signals  in  the  presence  of  ionizing  radiation,  and  (3)  the  design  of  a  star- 
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scanner  spacecraft  guidance  system  for  the  hostile  environment  of  space. 

A  number  of  more  complex  "contagious"  distributions  arising  from  multiplica¬ 
tive  processes  have  also  been  investigated  with  particular  emphasis  on  photon 
counting  and  direct-detection  optical  communications. 

In  the  nuclear  counting  domain  we  often  detect  ionizing  radiation  opti¬ 
cally  through  a  radiation-matter  interaction  in  which  a  single  high-energy 
particle  produces  a  shower  of  lower  energy  particles.  A  case  in  point  is  the 
scintillation  detector  which  is  a  combination  of  a  scintillation  crystal 
(e.g.,  NaI:Tl,  plastic)  with  a  photomultiplier  tube.^^  When  the  incident 
high-energy  particles  (e.g.,  8,  p,  y)  are  Poisson  distributed,  and  when  each 
particle  produces  a  Poisson  distribution  of  luminescence  photons  with  a 
specified  efficiency,  such  as  occurs  in  fluorescence  and  phosphorescence, 
the  resulting  photon-counting  distribution  representing  the  signal  is  the 
Neyman  Type-A.  If  Cerenkov  radiation  is  also  present,  primary  pulses  will  be 
registered  and  the  outcome  is  the  Thomas,  assuming  again  that  we  are  perform¬ 
ing  photon  counting.  It  is  usually  assumed  in  the  literature,  generally  tac¬ 
itly,  that  the  fixed  multiplicative  Poisson  distribution  describes  the  lumines¬ 
cence  photon  statistics  and  indeed,  when  the  counting  time  is  long  and  a  is 
large,  this  is  a  good  approximation.  But  if  photon  counting  is  used,  and  the 
luminescence  is  weak,  it  will  be  necessary  to  use  the  more  accurate  forms 
described  above.  Similar  statistics  will  apply  to  the  detection  of  cathodo- 
luminescence  and  photoluminescence. 

In  certain  applications  where  we  wish  to  count  the  photons  in  an  optical 
signal  (e.g.,  in  space),  it  sometimes  happens  that  the  distributions  discussed 
in  the  previous  paragraph  are  characteristics  of  the  noise  rather  than  of  the 
signal.  Viehmann  and  Eubanks  have  discussed  sources  of  noise  in 

photomultiplier  tubes  in  the  radiation  environment  of  space.  Such  noise  may 
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arise  from  several  mechanisms  such  as  luminescence  and  Cerenkov  radiation  in 
the  photomultiplier  window;  secondary  electron  emission  from  the  window,  photo¬ 
cathode,  and  dynodes;  Bremss trahlung  in  turn  causing  such  secondary  electron 
emission;  cosmic-ray  bursts;  and,  of  course,  thermionic  emission  dark  current. 
These  effects  clearly  degrade  both  the  dynamic  range  and  the  photometric  ac¬ 
curacy  of  low-light-level  measurements,  and  therefore  must  be  clearly  under¬ 
stood.  Photon  counting  can  be  particularly  advantageous  in  such  situations: 

<✓ 

Even  if  a  large  number  of  photoelectrons  are  produced  by  the  Cerenkov  radia- 
tion  arising  from  a  single  charged  particle  (which  constitutes  noise  in  this 
case),  they  will  be  counted  as  only  a  single  pulse  since  the  Cerenkov  radia¬ 
tion  emission  time  is  much  shorter  than  the  transit  time  in  the  photomultiplier. 

The  mathematical  properties  of  the  Neyman  Type-A,  Thomas,  and  fixed  multi¬ 
plicative  Poisson  distributions  are  very  useful  for  the  study  of  the  performance 
of  systems.  The  permanence  under  convolution  of  distributions  with  identical 

parameters  a,  a,  or  A  facilitates  their  use  in  statistical  detection  and  es- 

(13) 

timation  problems  for  signal  and  noise  distributions  of  arbitrary  means. 

But  it  is  the  convergence  in  distribution  of  the  Neyman  Type-A  and  Thomas  to 

the  Gaussian  that  can  simplify  calculations  quite  a  bit  in  many  neural-, 

nuclear-,  and  photon-counting  applications.  In  that  limit  (a,  a  finite, 

W  -*•  oo)  ,  the  Neyman  Type-A  (Thomas)  random  variable  with  arbitrary  a(a)  may 

be  added  together  with  any  number  of  arbitrary  random  variables  that  converge 

in  distribution  to  the  Gaussian  (e.g.,  Poisson,  binomial,  Neyman  Type-A, 

Thomas)  to  yield  an  overall  Gaussian  random  variable  with  mean  (n)  *  I^(n^) 

2  2 

and  variance  ((An)  )  =  I^((AnJ  ). 

This  is  an  important  result  because  it  means  that  the  statistical  deci¬ 
sion  theory  and  estimation  theory  literature  relating  to  Gaussian  random  vari- 
(14) 

ables  can  be  brought  to  bear  on  the  problem.  Some  results  that  follow 
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are  single- threshold  detection,  well-known  solutions  for  estimation  prob- 

(13) 

lems,  and  a  simple  form  for  the  just-noticeable  difference  or  detection 

law.(16) 

The  conveyance  in  distribution  to  the  Gaussian  is  a  very  helpful  tool, 
but  a  far  more  potent  characterization  is  needed  to  properly  deal  with 
performance  specifications  (e.g.,  error  rate)  for  a  system  incorporating  a 
signal  or  noise  specified  by  the  Neyman  Type-A  or  Thomas  distribution. 

We  have  recently  developed  a  generalized  performance  parameter  for 
single -threshold  detection  systems^^  and  this  should  prove  very  useful  in 
just  such  a  specification. 

The  performance  of  single-threshold  communication  systems  can  be  charac¬ 
terized  in  a  number  of  ways,  including  the  probabilities  of  error,  detection, 
and  false  alarm,  the  sensitivity  (detectability),  and  the  traditional  signal- 

to-noise  ratio  (SNR).  The  SNR,  which  is  a  ratio  of  signal  power  to  noise 

.  .  ,  .........  ( 18) -( 24) 

power,  has  a  variety  of  specific  definitions,  depending  upon  the 

particular  application,  is  easy  to  measure  experimentally  in  terms  of  peak 

or  root  mean-square  values,  and  provides  an  accurate  measure  of  performance 

for  constant  signals  in  Gaussian  noise.  In  the  Gaussian  case,  the  SNR  varies 

monotonically  with  the  probability  of  error  t,  and  is  therefore  sufficient 

to  describe  system  performance.  In  general,  however,  the  SNR  might  not  vary 

monotonically  with  e,  and  is  not  necessarily  sufficient  as  the  sole  measure 

(22) (25) (26) 

of  system  performance.  The  inadequancy  of  the  SNR  is  evident 

where  the  noise  is  non-Gaussian  or  signal  fluctuations  are  important,  such  as 

wave  noise  (photon  bunching)  in  photocounting  optical  communication  systems, 

(22)  (251 

and  has  been  observed  in  both  optical  homodvne  and  heterodyne  systems. 

From  the  preceding  discussion  it  is  evident  that  it  would  be  useful  to 
obtain  a  generalized  parameter  which  varies  monotonically  with  system  per- 
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formance.  We  have  developed^^  such  a  new  generalized  performance  parameter 
that  provides  a  more  consistant  measure  of  system  performance  than  the  conven¬ 
tional  SNR.  The  analysis  relies  on  a  transformation  of  random  variables  which 
converts  the  Bayes  risk  for  arbitrary  statistics  into  an  expression  involving 
error  functions  only.  A  class  of  parameters  p  which  vary  monotonically  with 
the  Bayes  risk  was  derived.  It  was  shown  explicitly  that  this  parameter  decreases 
monotonically  as  the  Bayes  risk  and  probability  of  error  increase,  and  that  it 
reduces  to  the  SNR  under  appropriate  conditions. 

For  example,  for  Neyman  Type-A  noise  of  mean  <n>  and  Neyman  Type-A  signal- 
plus-noise  of  mean  <s>  +  <n>  ,  the  traditional  SNR  has  the  form 

SNR  =  - — - r-rw  (1) 

[ (a+l)<n>]  L  1 

where  a  is  assumed  to  be  the  multiplication  parameter  of  the  signal  and  noise 
and  the  denominator  corresponds  to  the  noise  energy.  The  generalized  perfor¬ 
mance  parameter  for  this  case  is 


P 


_ 2_  /t  -  <(s+n)^2> 

0+1  l  1/2 


/t  -  </n> 

1/2 


<n> 


(2) 


which  is  substantially  different  from  the  conventional  SNR.  In  particular, 
a  situation  was  investigated  in  which  p  was  monotonic  with  system  performance, 
but  the  SNR  was  not.  The  result  is  shown  in  Fig.  2,  in  which  p  indeed  varies 
monotonically  with  system  performance,  whereas  the  SNR  does  not. 

The  above  analysis,  which  utilizes  normalizing  transformations,  is  also 

useful  in  the  direct  calculation  of  the  Bayes  risk,  probability  of  error,  and 

■  •  •  (27) 

sensitivity.  It  has  led  to  a  new  technique  “  for  evaluating  optical  communi¬ 
cation  system  performance.  This  technique  is  most  useful  in  situations  where 
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Linear  plot  of  probability  of  error  e  vs  the 
SNR  (dashed  curve)  and  the  generalized  performance 
parameter  (solid  curve)  for  the  case  of  a  single 
tion  of  a  Gaussian  signal  in  additive,  independent 
Gaussian  noise. 
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exact,  explicit,  closed  form  solutions  exist  for  the  photocounting  cumulants, 
but  do  not  exist  for  the  photocounting  distribution,  such  as  in  the  general¬ 
ized  case  of  communication  through  lognormal  atmospheric  turbulence.  Using 
this  technique,  theoretical  probability  of  error  curves  were  generated  for  com¬ 
munication  through  lognormal  atmospheric  turbulence,  for  a  superposed  coherent- 
in-chaotic  signal,  embedded  in  additive,  independent  Poisson  noise,  with  ar¬ 
bitrary  ratio  of  sampling  time  to  source  coherence  time,  arbitrary  ratio  of 
coherent  to  chaotic  component,  arbitrary  mean  frequencies  of  the  coherent  and 
chaotic  components,  and  where  the  chaotic  component  need  not  be  stationary 
and  may  have  arbitrary  spectral  distribution.  Since  no  solution  exists  for  the 
photocounting  distribution  itself  in  this  generalized  case,  the  corresponding 
performance  calculation  has  not  previously  been  possible.  As  an  example. 

Fig.  3  illustrates  the  probability  of  error  versus  SNR  for  a  lognormally  mod¬ 
ulated  superposed  coherent-in-chaotic  signal  embedded  in  additive  independent 
Poisson  noise.  Here  S  is  the  ratio  of  the  detector  counting  interval  to  twice 
the  coherence  time  of  the  source,  and  v  is  the  ratio  of  coherent  to  chaotic 
component.  The  case  described  applies  to  the  detection  of  radiation  originat¬ 
ing  from  a  multimode  laser  or  scattered  from  a  rough  target,  and  passing  through 
atmospheric  turbulence. 

*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  NSF-ENG78-26498  and  by  the  Jet  Propulsion  Laboratory,  California 
Institute  of  Technology,  under  Contract  Nos.  955335  and  955336. 


(1) 

J.  Neyman,  Ann.  Math. 

Stat.  10,  35  (1939)  . 

(2) 

W.  Feller,  Ann.  Math. 

Stat.  14,  389  (1943). 

(3) 

M.  Greenwood  and  G.  U. 

,  Yule,  J.  Royal  Stat.  Soc .  Ser. 

A. 

83, 

255  (1920). 

(4) 

F.  Eggenberger,  Mitteilungen  der  Vereinigung  Schweizerischer 
Mathematiker ,  31  (1924) . 

Versicherungs- 

(5) 

F.  Eggenberger  and  G. 

Polya,  Zeitschrift  fllr  Angewand 

te 

Mathematik  und 

Mechanik,  1,  279  (1923) . 


19 


PROBABILITY  OF  ERROR 


Figure  3:  Probability  of  error  Pr(e)  vs  SNK  y  for  single-detector 
Neyman-Pearson  photocounting  receiver.  A  lognormal ly 
modulated  superposed  coherent-in- chaotic  signal  (of 
Lorentzian  spectrum  with  w  =  is  embedded  in  additive 

independent  noninterf eringCPoisson  noise.  Curves  labelled 
“  correspond  to  6»1,  curves  labelled  0.4  corresponded  to 
6  <<  1  and  y  =  4,  and  curves  labelled  0.0  correspond  to 
8  <<  1  and  y  =  0.  Here  <n)  =  1,  t  *  6,  and  Pp  '  5.94  x 
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C.  PHOTON  COUNT  ANTIBUNCHING  FOR  OPTICAL  COMMUNICATIONS,  RADAR, 

IMAGING,  AND  SPECTROSCOPIC  APPLICATIONS  * 

(M.  C.  Teich,  G.  Vannucci,  N.  Pierce,  R.  O'Connor) 

(Principal  Investigator:  M.  C.  Teich  (212)  280-3117) 

The  probability  distribution  for  a  dead-time-modif  ied  pulse  counter^ 
has  been  studied  by  a  number  of  researchers  in  a  broad  variety  of  disciplines, 
such  as  photon  counting, nuclear  counting, (12)  an(j  neural  count¬ 
ing.  Many  cases  have  been  studied  in  detail,  including  par- 

alyzable  and  nonparalyzable  counting  under  blocked,  unblocked,  and  equilibrium 
conditions.  Attention  has  also  been  given  to  the  gradual  recovery  (sick-time) 
system  and  to  the  variable  dead-time  case/^  ’  Muller  has  summarized  the 

results  of  a  number  of  authors  an(j  has  compiled  a  comprehensive  bib¬ 

liography  on  dead-time  effects. 

Although  most  of  the  work  cited  above  is  applicable  only  when  the  input 
to  the  counter  is  a  Poisson  point  process  with  constant  rate,  a  few  results 
are  also  available  for  the  case  in  which  the  rate  is  not  constant.  Cantor 
and  Teich, Teich  and  McGill,  ^  and  Bedard^^  present  expressions  for  the 
photocounting  distribution  when  the  intensity  cf  the  light  is  a  random  pro¬ 
cess,  with  the  sampling  time  much  smaller  than  the  coherence  time  of  the  light. 

In  a  prior  work,^^  we  obtained  expressions  for  the  mean  and  the  vari¬ 
ance  of  the  number  of  events  in  a  fixed  sampling  time  for  a  nonparalyzable 
dead-time  counter  when  the  input  process  is  Poisson  with  a  rate  that  is  a 
known  function  of  time.  We  have  now  extended  those  results  to  the  case  in 
which  the  rate  of  the  input  process  is  an  arbitrary  stochastic  process,  under 
the  constraint  that  it  vary  slowly  with  respect  to  the  duration  of  the  dead 
time.  No  constraints  on  the  sampling  time  are  imposed.  We  have  shown 
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how  Che  dead-time-modif ied  mean  and  variance  depend  on  the  first-  and  second- 
order  statistics  of  the  rate.  In  the  context  of  the  detection  of  light,  we 
have  obpained  explicit  expressions  for  the  photocount  mean  and  variance  for 
chaotic  radiation  of  arbitrary  spectral  properties.  Although  the  magnitude 
of  the  variance  depends  explicitly  on  the  power  spectrum  of  the  radiation  in 
the  absence  of  dead  time,  it  has  been  shown  that  the  decrease  in  variance  with 
increasing  dead  time  is  basically  spectrum  independent.  Analytical  and  gra¬ 
phical  results  for  the  dead-time-modified  mean  and  variance  are  therefore  now 
available  for  chaotic  light  of  Lorentzian  and  Gaussian  spectra. 

The  principal  results  of  this  effect  comprise  analytic  expressions,  com¬ 
plemented  by  computer  results,  for  the  mean  and  variance  of  a  dead-time- 
modified  Poisson  counting  process  driven  by  a  continuous  rate  process  repre¬ 
sentative  of  chaotic  radiation.  (It  has  been  assumed  that  the  nonparalyzable 
dead  time  is  small  in  comparison  with  the  coherence  time  (inverse  bandwidth) 
of  the  light.)  In  this  limit,  as  well  as  in  the  absence  of  dead  time,  the 
photocount  mean  is  explicitly  shown  to  depend  only  on  the  first-order  statis¬ 
tics  of  the  radiation  (see  Fig.  1).  As  was  expected,  the  counting  efficiency 
lies  below  that  for  a  constant  intensity  source  for  all  values  of  the  dead¬ 


time  parameter  At.,  where  A  is  the  mean  dn  g  rate  and  T,  is  the  dead  time, 
d  d 

The  photocount  variance,  on  the  other  hand,  exhibits  a  strong  dependence 
on  the  spectral  properties  of  the  radiation  even  in  the  absence  of  dead  time 
(see  Fig.  2).  The  forms  of  the  dead-time-modif ied  analytic  expressions  for 
the  cases  of  Lorentzian  and  Gaussian  spectra  are  therefore  substantially  dif¬ 
ferent.  Nevertheless,  numerical  computer  evaluation  of  these  expressions 

(variance  as  a  function  of  At.)  demonstrates  that  the  effect  of  dead  time  on 

d 

the  variance  is  essentially  independent  of  the  spectral  shape  of  the  light. 


(see  Fig.  3)  although  a  dependence  can  be  discerned  to  develop  at  higher 
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COUNTING  EFFICIENCY 


Figure 


1:  Counting  efficiency  (E[n]/AT)  versus  Aid,  where  A  is  the 

driving  rate  and  is  the  dead  time.  Curves  are  for  a  Poisson 
process  in  which  the  rate  is  constant  (solid  curve)  and  for 
chaotic  radiation  of  arbitrary  spectrum  (dashed  curve).  It  is 
clear  that  the  efficiency  is  significantly  reduced  (up  to  24%) 
when  the  rate  is  not  constant.  The  results  are  analogous  to 
those  presented  in  Fig.  1  of  Ref.  16  for  a  rate  that  is  a 
known  function  of  time. 


vsr  fn)  /XT 


100 


CONSTANT  INTENSITY 


LORENTZIAN  SPECTRUM 


GAUSSIAN  SPECTRUM 


Figure  2:  Ratio  of_dead-time-modif ied  variance  to  unmodified  mean 

(var[n]/XT)  versus  Xt^  for  a  Poisson  process  with  constant  rate 
(solid  curve),  for  chaotic  radiation  with  Gaussian  spectrum 
(dashed  curves),  and  for  chaotic  radiation  with  Lorentzian 
spectrum  (dotted  curves).  The  values  of  the  parameters  for 
the  chaotic  radiation  curves  are  FT  *=  100,  XT  =  10,000 
(upper  curves)  and  TT  *  2,  XT  «  40  (lower  curves).  The 
vertical  separation  between  any  of  the  curves  for  chaotic 
radiation  is  essentially  independent  of  Xtj,  an  indication 
that  the  different  curves  are  scaled  versions  of  one  another. 
Thus  the  normalized  variance  will  be  essentially  the  same  for 
all  four  curves,  as  is  demonstrated  in  Fig.  3. 
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NORMALIZED  VARIANCE 


CONSTANT  INTENSITY 
CHAOTIC  RADIATION 


Figure  3:  Normalized  variance  (ratio  o_f  dead-time-modified  variance  to 
unmodified  variance)  versus  Curves  are  for  a  Poisson 

process  with  constant  rate  (solid  curve)  and  for  chaotic 
radiation  (dashed  curve).  All  results  were  obtained  under  the 
assumption  that  the  dead  time  is  much  shorter  than  the  coherence 
time  (tjj  <<  1/D,  in  which  case  the  effect  of  dead  time  turns 
out  to  be  virtually  independent  of  spectral  shape  (see  texO . 

The  normalized  variance  is  therefore  the  same  function  of 
for  all  types  of  chaotic  radiation.  The  results  are  similar 
in  character  to  those  presented  in  Fig.  2  of  Ref.  16  for  a 
rate  that  is  a  known  function  of  time. 


values  of  at^.  This  can  be  understood  by  noting  that  our  results  are 

applicable  only  when  the  dead  time  is  much  shorter  than  the  coherence  time 

(t  <<  1/r) .  For  sufficiently  small  t,,  the  rate  of  the  Poisson  process  can 
d  d 

be  locally  approximated  as  being  constant  so  that  the  dead-time-modified 
point  process  has  a  rate  that  is  locally  a  function  of  the  unmodified  rate 
only  and  as  such  is  not  influenced  by  the  second-order  statistical  properties 
of  the  radiation. 

In  light  of  our  results,  we  may  expect  that  the  dead-time-modif ied 

count  variance  will  exhibit  an  increasing  dependence  on  the  spectrum  as  t. 

d 

becomes  larger;  indeed,  this  dependence  should  become  substantial  as 

approaches  1/T.  An  analysis  in  that  regime,  although  difficult  to  carry  out, 

might  therefore  lead  to  the  intriguing  possibility  of  using  a  dead-time 

counter  to  extract  spectroscopic  information  about  a  radiation  source.  The 

notion  of  dead-time  spectroscopy  is  not  unlike  photocounting  spectroscopy,^^ 

but  the  upper  limit  of  bandwidths  that  could  be  measured  by  the  technique 

would  be  'v  10  GHz.  This  value  is  higher  than  the  upper  limit  for  photocounting 

spectroscopy,  n.  100  MHz,  because  the  dead  time  can  in  principle  be  made  about 

a  factor  of  100  smaller  than  the  10-nsec  minimum  counting  time  in  which 

variations  can  be  detected  with  commercially  available  fast  electronic 

circuits.  The  region  of  usefulness  of  dead-time  spectroscopy  would  probably 

be  similar  to  that  of  the  Fabry-Perot  etalon. 

Finally,  we  point  out  that  the  dead-time-modif ied  count  mean  may  also 

depend  on  the  spectral  properties  of  the  light  for  t  n,  1/r,  in  which  case 

d 

(18) (19) 

dead-time  spectroscopy  could  be  implemented  quite  easily.  Kikkawa  et  al. 
have  calculated  the  efficiency  (normalized  mean)  of  a  dead-time-modif ied 
photon  counter  for  Lorentzian  light  and  have  discussed  the  dependence  of 
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this  quantity  on  spectral  parameters.  Their  results  are  valid  only  in  the 


limit  where  the  coherence  time  is  much  smaller  than  the  mean  time  interval 

between  pulses,  however,  and  this  is  not  a  useful  region  [see  remarks 

following  Eq.  (10)  in  Ref.  18]. 

*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  NSF-ENC78-26498 . 
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II.  PHYSICAL  PROPERTIES  AND  EFFECTS  IN  ELECTRONIC  MATERIALS 

A.  CARRIER  TRANSPORT  ACROSS  HETEROJUNCTION  INTERFACE  * 

(C.  M.  Wu,  E.  S.  Yang) 

(Principal  Investigator:  E.  S.  Yang  (212)  280-3120) 

In  our  previous  study,  we  have  calculated  the  quantum  transmission  co- 

eff icient^^  by  solving  Schrodinger 's  equation  across  the  heterojunction 

(2) 

potential  barrier  in  which  the  perpendicular  momentum  is  conserved.  Ap¬ 
plication  of  our  model  to  an  N-p  hetero jenction  shows  that  its  forward  cur¬ 
rent  may  be  dominated  by  either  thermionic  emission  or  tunneling  current 
if  the  effective  masses  of  both  sides  of  the  junction  are  much  smaller  than 

(3) 

the  free  electron  mass.  In  most  cases,  however,  the  thermionic  emission 

current  is  large  compared  with  the  diffusion  current  so  that  the  forward 
current  is  controlled  by  minority  carrier  diffusion.  As  for  a  p-P  hetero¬ 
junction,  rectifying  characteristics  are  obtained.  However,  the  reverse 
current  appears  to  be  non-saturated ,  but  the  forward  current  remains  con¬ 
stant  at  large  forward  voltage.  These  results  are  in  agreement  with  our 

(4) 

experiment  data  and  have  been  confirmed  by  others. 

An  important  consideration  rtat  included  in  the  previous  calculation  is 
the  charge  in  the  heterojunction  interface  Kroemer  et  showed  that  the 

charge  of  the  interface  state  should  be  taken  into  account  for  the  band 
diagram  calculation.  Since  the  potential  notch 'formed  at  the  narrow  band- 
gap  side  can  trap  injected  carriers  under  forward  bias,  we  find  it  necessary 
to  include  the  electron  accumulation  there.  In  our  work,  an  A1  -,,-Ga  -As- 
GaAs  N-p  hetero junction  is  fabricated  by  the  LPE  technique.  The  N-type  epitaxial 
layer  is  grown  on  the  p-tvpe  substrate  at  789°Cwith  a  cooling  rate  of  0.2°C/min. 
Multilayer  ohmic  contacts,  Sn-Au-Ni  and  7n-Au-Ni,  are  deposited  on  the  n-sideand  p-side 
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( ev ) 


respectively.  The  experimental  devices 's  current-voltage  characteristics 
deviate  significantly  from  the  exponential  form,  particularly  a  high  bias 
voltage  as  if  a  large  series  resistance  were  in  the  structure.  The  series 
resistance  thus  obtained  is,  however,  too  large  for  the  given  bulk  doping  ana 
does  not  agree  with  that  of  the  control  sample.  By  taking  into  account  the 
carrier  accumulation  in  the  potential  notch,  we  are  able  to  explain  the  na¬ 
ture  of  current  suppression  through  the  modulation  of  the  effective  barrier 
height.  The  result  explains  the  discrepancy  between  the  theoretical  predic¬ 
tion  and  experimental  data  reported  by  Kroemer  et  al . ^  In  addition,  the 
activation  energy  from  the  temperature  data  as  a  function  of  bias  is  found 
to  be  a  convenient  way  of  measuring  the  conduction  band  discontinuity  at  the 
heterojunction  interface  (Fig.  1). 
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B.  ELECTRONIC  PROPERTIES  OF  POLY  Si-Si  DIODES* 

(C.  M.  Wu,  E.  S.  Yang) 

(Principal  Investigator:  E.  S.  Yang  (212)  280-3120) 


Low  temperature  processes  for  forming  a  rectifying  junction  are  attrac¬ 
tive  in  the  production  of  low-cost  solar  cells^^  ^  and  in  possible  applica¬ 
tions  in  low-performance  active  structures  .  In  this  letter,  we  report 

a  new  low- temperature  method  of  preparing  an  aluminum-diffused  p+-n  junction 
by  making  use  of  the  high  diffusivity  of  A1  in  poly-silicon .  Experimental 

I-V  and  C-V  curves  are  shown  along  with  a  simple  model  to  explain  the  data. 

The  fabrication  procedure  of  the  aluminum-diffused  polysilicon  p+-n 
junction  starts  with  E-beam  evaporated  multi-layers  (100  X  Si-2000  A1-1000X  Si) 
on  an  n-type  (111)  silicon  substrate  (1^2  cm)  as  shown  in  the  insert  of 
Fig.  1,  It  is  followed  by  heat-treatment  (HT)  at  600°C  in  flowing  hydrogen 
for  one  hour.  Thus,  an  Al-riched  p-tvpe  polycrystalline  layer  is  formed, 
which  may  be  used  as  an  Al-diffusion  source  if  a  subsequent  drive-in  process 
is  employed. ^  It  is  important  to  note  that  the  cap  silicon  layer  and  the 
hydrogen  gas  are  employed  to  prevent  the  aluminum  film  from  reacting  with 

( g ) 

the  environment  at  the  HT.  The  HT  temperature  should  not  be  higher  than 

the  A1  melting  point  (660°C)  .  Here,  the  p+-n  junction  under  investigation 
is  a  device  without  the  drive-in  step  and  the  p+-layer  is  a  polycrystalline 
silicon,  layer.  After  making  ohmic  contact  on  the  back  side,  we  measure  the 
I-V  and  C-V  characteristics  at  different  temperatures.  We  also  find  that 
a  diode  fabricated  with  a  thin  A1  film  (100  X)  gives  similar  electric  pro¬ 


perties  with  good  transparency,  but  with  a  large  series  resistance. 

The  room  temperature  I-V  curves  are  shown  in  Fig.  1.  The  forward  current 

magnitude  is  much  higher  than  that  obtained  in  the  single  crystal  p+-n  junc- 
-10  2 

tion  (J^  v  10  amp/cm  ),  A  plausible  interpretation  is  that  the  p-tvpe 


polycrystalline  layer  is  full  of  localized  states  caused  by  grain  boundaries 
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Figure  1:  The  I-V  characteristics  of  the  p  -n  junction  at  room 

temperature.  The  insert  is  a  structure  of  the  diode  before 
the  heat-treatment. 
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and  dislocation  networks  where  electrons  recombine  after  injection  from  the 


single-crystal  n-tvpe  region. 

Since  the  silicon  film  is  evaporated  without  substrate  heating,  the  grain 
size  is  very  small,  and  the  mobility  is  low.  As  a  result,  the  diffusion  length 
of  electrons  in  Al-doped  polvsilicon  layer  is  short,  and  carrier  recombina¬ 
tion  mostly  occurs  at  the  interface  between  the  polvsilicon  and  the  single 
crystal  substrate.  If  we  assume  a  constant  recombination  velocity  S,  the  total 
forward  current  is  equal  to  the  electron  current,  which  is  given  by 

q4> 

J  =  qSNd  exp(-  -j^)  (1) 


where  N,  is  the  donor  density,  and  <J>_  is  the  built-in  potential  only  in  the 
d  D 

n-side.  However,  is  a  function  of  the  bias  voltage  and  is  given  by 


*  l+N./N 
d  a 

where  V  is  the  zero-bias  built-in  potential  across  the  whole  junction, 

D  1 

is  the  bias  voltage,  and  N  is  the  doping  in  the  polysilicon  layer. 

a 

Therefore,  the  forward  current  can  be  rewritten  as 


(2) 

V 


J  =  Jo  eXpd}  (3) 

and 

Jo  =  qSNd  exp  [-  <4) 

where  n  is  equal  to  1  +  N./N  .  The  C-V  and  I-V  characteristics  are  measured 

d  a 

at  different  temperatures  from  -60°C  to  60°C.  The  slope  of  the  1/C3  curve 

versus  the  bias  voltage  does  not  vary  with  temperature  and  gives  the  impurity 

density  of  the  lightly-doping  side,^  which  is  6.7  x  1033  cm  3  and  is  con- 

2 

sistent  with  the  impurity  of  the  n-side.  The  extrapolation  of  the  1/C  curve 
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intersects  the  voltage  axis  at  <J>.  and  yields  the  built-in  potential  V,  .  = 

mt  b  1 

(9) 

+2  kT/q  at  zero  bias.  The  zero-bias  current  J  can  be  obtained  by 
int  o 

the  extrapolation  of  the  I-V  curve  to  the  current  axis.  Table  I  lists  the 

experimental  values  of  J  and  n  from  the  I-V  data  and  <f> .  from  the  C-V  data. 

o  mt 

It  is  noted  that  the  ideality  factor  (n-value)  decreases  as  the  temperature 

increases.  According  to  the  relation  of  the  doping  density  to  the  n-value, 

the  impurity  density  of  the  p-tvpe  region  can  be  deduced  and  is  listed  in 

Table  I.  It  is  observed  that  the  active  impurity  density  in  the  polysilicon 

increases  with  the  temperature.  This  is  probably  due  to  the  high  ionization 

energy  of  aluminum  in  silicon  (0.057  eV) ,  As  a  result,  some  of  the  aluminum 

atoms  are  not  electrically  active  in  the  polysilicon. 

By  considering  the  variation  of  the  Fermi  potentials  with  the  temperature 

on  both  sides,  we  obtain  the  sum  of  the  built-in  potential  v^  and  the  Fermi 

potentials  which  is  equal  to  the  energy  gap  with  an  error  of  2%.  Therefore, 

the  zero-bias  current  density  J  has  a  temperature-dependent  activation  energy 

o 

4 

which  is  predicted  by  Eq.  (4).  Using  a  recombination  velocity  of  5  x  10  cm/sec 

at  0°C  and  assuming  S  to  be  proportional  to  the  thermal  velocity,  we  obtain 

the  theoretical  values  of  J  by  Eq.  (4),  as  shown  by  the  solid  line  in  Fig.  2. 

o 

We  find  that  our  simple  model  is  in  agreement  with  the  experimental  data  ob- 
tained  from  I-V  data  when  the  temperature  is  higher  than  -40°C.  In  Eq .  (1), 
we  consider  the  electron  current  as  the  dominant  component  and  this  might  not 
be  true  below  -40°C .  For  example,  the  recombination  in  the  depletion  region 
may  be  higher  than  the  electron  current.  However,  the  conduction  mechanism 
does  not  conform  to  the  simple  model  at  temperatures  below  -40°C. 

According  to  the  above  discussion,  the  injected  electrons  from  the  n-side 
are  recombined  at  the  interface  between  the  p+-polvsil icon  layer  and  the  crys¬ 
talline  substrate  and  the  recombination  velocity  S  here  should  be  in  the  order 
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Table  I.  Experimental  data  from  I-V  and  C-V  curves 


and  the  Fermi-potential  at  the  n-side. 
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of  that  found  at  the  grain  boundaries  of  polysilicon.  Our  empirical  value  of 
S  is  consistent  with  Ghosh*s  result  (3  x  lO^4  cm/sec) 

In  summary,  we  have  developed  a  method  using  aluminum  as  the  p-type 
impurity  to  form  a  polysilicon-silicon  p+-n  junction  at  low  temperatures. 
Because  of  the  simplicity  in  the  fabrication  process,  this  rectifying  junc¬ 
tion  may  provide  a  low-cost  method  of  producing  solar  cells.  However,  thinner 
aluminum  (  <  50  X)  and  polysilicon  ( <  100  X)  layers  will  be  necessary  to  im¬ 
prove  the  transparency  and  to  lower  the  series  resistance.  Further  studies 
are  required  to  obtain  the  optimum  structure, 

*This  research  was  also  supported  by  the  Department  of  Energy  under  Contract 
DOE/SERI  xz-4>- 9226. 
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c.  LASER  PROCESSED  SCHOTTKY  BARRIERS* 

(E.  S.  Yang,  C.  M.  Wu,  H.  J.  Vollmer,7  T.  0.  Sedgwick,  +  R.  T.  Hodgson4”) 
(Principal  Investigator:  E.  S.  Yang  (212)  280-3120) 

We  report  a  novel  technique  for  fabricating  Schottky  barrier  junctions 
by  laser-irradiation  process  (LIP) .  Both  pulsed  and  cw  irradiation  have  been 
employed.  The  electrical  characteristics  of  the  LIP  diodes  are  presented. 

The  devices  w>re  prepared  by  depositing  2000-X-thick  aluminum  dots  on 

n/n+  (100)  silicon  wafers.  The  resistivity  of  the  n-epilayer  was  measured 

to  be  0.4  ftcm.  Half  of  the  A1  dots  were  covered  by  a  100-X-thick  e-beam- 

evaporated  silicon  film.  Gold  was  evaporated  on  the  back  side  of  the  wafer 

to  provide  a  good  ohmic  contact  to  the  n+  substrate.  The  current-voltage 

characteristics  of  these  devices  were  tested  before  LIP  and  were  found  to  be 

ohmic  with  a  low  resistance.  The  devices  were  then  subjected  to  either  the 

pulsed  or  cw  laser  irradiation.  In  the  pulsed  LIP,  we  used  a  Q-switched 

Nd:YAG  laser  (X  =  1.06  urn)  with  a  pulse  width  of  60  ns  and  energy  density  of 
2 

approximately  1  J/cm  .  The  substrate  is  not  heated.  The  aluminum  dots  without 
the  silicon  top  layer  did  not  show  any  change  in  their  electrical  character¬ 
istics  although  damage  in  aluminum  was  seen.  However,  dots  with  a  silicon 
cover  layer  showed  striking  transformation.  In  the  cw  LIP,  we  used  an  argon 
laser  (4880  and  5145  X)  operated  between  10  and  13  W  and  focused  to  a  spot 
size  of  50  pm.  The  spacing  between  adjacent  scanning  lines  is  20  pm  so  that 
complete  coverage  of  the  surface  by  the  laser  beam  is  ensured.  The  substrate 
temperature  was  set  at  350°C  (the  minimum  temperature  in  the  experimental 
setup)  and  the  stepping  motor  movement  (with  a  speed  of  2  cm/sec)  gives  rise 
to  overlapping  spots  of  laser  irradiation.  Again,  the  dots  without  the  sil¬ 
icon  cover  layer  appeared  to  be  unaffected  bv  the  laser  beam.  In  contrast, 
rectifying  characteristics  are  plotted  in  Fig.  1  for  diodes  obtained  by  both 
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Figure  1:  Current-voltage  characteristics  of  Schottkv  diodes:  (a)  before 
laser  irradiation,  (b)  after  pulse  Nd:YAG  laser  LIP,  and  (c) 
after  cw  argon  laser  LIP. 

cw  and  pulsed  LIP,  The  calculated  barrier  height  form  capacitance-voltage 
data  is  1  eV,  which  is  very  close  to  the  band-gap  energy.  The  cw  LIP  pro¬ 
duces  diodes  with  a  lower  current  at  the  same  bias  voltage.  The  higher  cur¬ 
rent  in  the  pulsed  LIP  diodes  is  caused  by  incomplete  treatment  of  the  dots 
since  manual  scanning  was  employed.  The  aluminum  layer  appears  to  have 
melted,  and  the  surface  damage  is  severe  and  not  uniform. 

For  aluminum  dots  without  a  silicon  cover,  it  appears  that  the  laser 
light  is  not  absorbed  because  of  the  reflection  of  the  aluminum  surface. 

The  silicon  cover  facilitates  strong  absorption  so  that  heating  of  the  Al-Si 
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interface  is  effected.  It  is  reasonable  to  assume  that  the  surface  tempera¬ 
ture  is  higher  than  the  Al-Si  eutectic  so  that  the  aluminum-silicon  interface 
melts  and  recrvstallizes.  Furthermore,  part  of  the  aluminum  may  have  diffused 
into  silicon  to  form  a  p+  region.  The  resulting  structure  is  an  Al-P+  -n  Si 
Schottky  barrier  with  a  thin  p+  layer  located  within  the  junction  space-charge 
region  underneath  the  surface.  The  barrier  height  of  the  untreated  Al-n  Si 
contact  is  very  low  because  of  the  low  work  function  of  aluminum;  therefore 
the  contact  is  ohmic,  as  shown  in  curve  (a)  of  Fig.  1.  The  introduction  of 
aluminum  forms  a  thin  p+  layer  and  gives  rise  to  p-n  junction-like  charac¬ 
teristics  Figs.  1(b)  and  1(c).  However,  the  new  rectifier  differs  from  the 
p-n  junction  in  that  the  p+  region  is  completely  depleted.  This  leads  to 
the  features  of  high  barrier  height  and  an  ideality  factor  substantially 
larger  than  unity.  These  results  are  similar  to  those  for  thermally  annealed 

Schottky  diodes ^  and  are  in  agreement  with  the  simple  model  proposed  re- 

(2) 

cently  by  Van  der  Ziel.  The  large  current  in  the  low-bias  regime  with 
n*3  is  likely  to  be  the  result  of  surface  recombination. 

In  conclusion,  we  have  found  that  a  rectifying  junction  can  be  formed  by 
laser  irradiation.  The  ohmic  contact  of  A1  on  n-Si  has  been  converted  to  a 
Schottky  barrier  with  a  large  barrier  height.  The  characteristics  of  these 
diodes  are  comparable  to  those  of  thermally  annealed  Al-n-Si  Schottky  bar¬ 
riers.  Because  of  its  inherent  advantage  of  low  temperature  (except  in  the 

(3  j 

junction  region),  LIP  technology  may  be  important  in  device  applications. 

*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  NSF-ENG-78-19426 . 
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D.  ELECTRON'  TRAPS  IN  ZINC  SILICATE  PHOSPHORS 

(E.  S.  Yang,  J.  M.  Brownlos  ) 

(Principal  Investigator:  E.  S.  Yang  (212)  280-3120) 

In  the  study  of  Zn^SiO^:  Mn  (willemite)  Phosphors,  it  has  been  observed 
that  the  phosphorescent  decay  can  be  described  by  three  or  more  time  constants, 
the  shortest  of  which  is  typically  between  10  and  20  ms.  This  short  time 

|  j 

constant  is  believed  to  be  limited  by  the  intrinsic  Mn  radiative  transi¬ 
tion  rate.  The  long  decay  time  constants  are,  however,  related  to  the  release 
of  electrons  from  deep  traps.  Although  deep  traps  have  been  determined  by 
thermo luminescence  experiments,  there  has  been  little  published  work  attempt¬ 
ing  to  correlate  these  trap  distributions  to  the  phosphorescent  decay  data. 

The  difficulty  is  that  the  phosphorescent  decay  time  in  willemite  crystals 
is  found  to  be  weakly  dependent  on  temperature.^  This  is  contrary  to  the 
strong  temperature  dependence  of  decav  time  expected  from  thermally  activated 
traps.  Here,  we  present  measurements  on  the  thermal  dependence  of  decay 
curves  of  willemite  phosphors  for  correlation  with  trap  distribution  mea¬ 
surements  made  on  the  same  samples. 

In  our  experiments,  we  have  used  different  types  of  willemite  samples 
including  commercial  PI  and  P39  as  well  as  specially  prepared  phosphores  pre¬ 
pared.  The  powder  phosphor  is  placed  in  a  sample  holder  within  an  evacu¬ 
ated  chamber  and  its  temperature  can  be  varied  between  100  K  and  670  K.  The 
excitation  source  is  a  dc  mercury  lamp  whose  light  passes  through  an  input 

monochromater  set  at  2536  X,  Details  of  the  experimental  setup  have  been 

(2) 

given  in  a  previous  paper.  With  the  same  experimental  arrangement  we  ob¬ 

tain  the  thermoluminescent  (glow  curve)  and  phosphorescent  decay  data.  A 
typical  glow  curve  of  a  zinc  silicate  phosphor  has  four  glow  peaks,  and  the 
corresponding  trap  depths  are  calculated  to  be  0.13,  0.24,  0.53,  and  0.87  eV . 

We  then  measure  the  phosphorescent  decay  at  the  temperatures  corresponding 
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to  Che  minima  in  the  glow  curve  so  as  to  determine  the  decay  time  constant  for 
the  trap  next  to  the  minimum  on  the  high  temperature  side.  The  results  are 
shown  in  Fig.  1.  It  is  noted  that  the  decay  time  cannot  be  described  by  a 
single  exponential  term  and  the  initial  decay  rate  does  not  change  by  more 
than  a  factor  of  2  for  a  wide  range  of  temperatures.  In  addition,  the  decay 
time  increases  between  -170°C  and  -48°C,  then  drops  at  0°C  and  eventually 
goes  back  up  at  46°C.  This  fluctuating  behavior  is  observed  in  all  samples 
showing  multiple  peaks  in  the  thermoluminescent  experiment. 

The  B725  sample  has  very  different  characteristics.  This  phosphor  is 
prepared  by  using  a  mixture  of  ZnO,  SiO^,  and  MnO  fired  at  1150°C  for  15 
minutes.  It  is  subsequently  ground  to  powder  and  refired.  The  glow  curve 
of  B725  has  a  single  peak  which  is  unique  among  willemite  crystals,  and  it 
has  an  activation  energy  of  0.433  eV ,  The  decay  curves  depicted  in  Fig.  2 
show  a  monotonic  decrease  of  the  first  time  constant  with  increasing  tempera¬ 
ture.  In  other  words,  when  there  is  only  one  trap,  the  decay  time  constant 
does  not  increase  at  any  temperature  in  contrast  with  the  multiple-trap  be¬ 
havior  of  PI  shown  in  Fig,  1. 

When  arsenic  is  introduced  in  small  dosages  ('v  0.1%),  the  willemite  is 
known  as  P39  and  the  characteristics  are  substantially  modified.  The  glow 
curve  of  P39  is  found  to  have  one  main  peak  near  0°C  and  two  other  widely 
separated  peaks.  The  trap  energies  are  calculated  to  be  0.334,  0.63,  and 
1.78  eV,  and  the  temperature  dependence  of  the  decay  curves  is  similar  to 
that  shown  in  Fig.  1, 

By  examining  the  change  of  intensity  and  of  half-width  of  the  glow 
curves,  we  find  that  the  data  of  PI  follow  closely  the  pattern  described  by 
Cher  and  Kivits.^^  If  the  thermal  release  of  electrons  from  the  trap 
is  raonomolecular,  then  the  decay  should  have  a  single  time  constant. 
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Figure  2:  Phosphorescent  decay  curves  of  B625  showing  monoto 


However,  Che  data  shown  in  Fig.  1  clearly  indicate  non-exponential  decay. 

Let  us  consider  a  model  where  the  thermally  released  electrons  are  contri¬ 
buted  by  more  than  one  trap  at  any  time  such  that  the  decay  curve  is  the  re¬ 
sult  of  composite  luminescence.  Then  we  can  argue  that  the  fluctuating  na¬ 
ture  of  the  decay  time  constant  has  the  appearance  of  temperature  indepen¬ 
dent  characteristic  and  is  caused  by  the  overlapping  of  thermal  release  of 
electrons  from  energetically  adjacent  traps.  As  the  temperature  is  increased, 
the  most  shallow  trap  is  the  first  to  give  up  its  electrons,  and  its  decay 
time  decreases  with  increasing  temperature.  However,  if  the  next  deeper 
trap  starts  to  contribute  before  the  first  trap  is  completely  empty,  the 
overall  decay  time  will  increase  rather  than  decrease  since  the  deeper  trap 
has  a  longer  decay  time  constant.  If  the  traps  are  energetically  close, 
the  decay  time  is  related  to  the  integrated  luminescence  from  the  traps 
involved  such  that  the  contribution  of  each  trap  is  difficult  to  isolate. 

This  is  the  case  for  the  PI  phosphor.  The  rate  limiting  process  can  be 
identified  unambiguously  if  there  is  only  one  trap  such  as  in  B725  or  if 
the  traps  are  widely  separated  such  as  in  P39. 

For  the  two  low  temperature  traps  in  P39,  the  trap  depths  may  be  ob¬ 
tained  by  representing  the  decay  curve  as  a  summation  of  exponential  terms 
and  plotting  the  time  constants  as  functions  of  temperature.  Using  only 
two  exponential  terms  at  each  temperature,  we  obtain  the  time  constants  and 
plot  the  result  in  Fig.  3.  The  agreement  between  the  thermal  glow  curve 
and  decay  time  experiments  is  quite  good,  since  the  trap  depths  and  the  fre¬ 
quency  factors  agree  within  one  order  of  magnitude.  This  agreement  is  re¬ 
markable  in  view  of  the  fact  that  the  frequency  factor  is  exponentially  re¬ 
lated  to  the  temperature.^  However,  other  powders  show  less  favorable 
correlation  as  a  result  of  the  non-exponential  decay  due  to  tunneling  or 
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rUec) 


F igure 


3:  Plot  of  the  first  and  second  decay  time  constants  vs  inverse 

temperature.  The  energy  and  S-factor  obtained  from  the  glow 
curve  method  are  E^=0.33A  eV,  Si=2.4xl0^  1/sec  and  E2=0.63  eV 
and  S2=9x10^  1/  sec. 


bimolecular  mechanisms. 
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E.  ENERGY  AND  ELECTRIC  FIELD  DEPENDENCE  OF  Si-S102  INTERFACE  STATE 
PARAMETERS  BY  OPTICALLY  ACTIVATED  ADMITTANCE  EXPERIMENTS 

(T.  C.  Poon,  H.  C.  Card,  S.  M.  So,  W.  Hwang) 

(Principal  Investigator:  W.  Hwang  (212)  280-3115) 


Si-SiO^  interface  states  have  been  the  subject  of  considerable  experi¬ 
mental^^  ^ ^  and  theoretical^^  study,  and  understanding  of  the  origin 

of  these  states  depends  upon  additional  experimental  information  beyond  that 

(10) 

obtained  by  conventional  admittance  measurements.  In  a  previous  paper, 
we  have  described  a  new  experimental  method  by  which  to  obtain  the  interface 
state  parameters  in  a  metal  oxide  semiconductor  (MOS)  structure,  which  in¬ 
volves  measurements  of  the  ac  admittance  as  a  function  of  frequency  under 
a  background  of  optical  illumination.  This  technique  provides  for  an  ex¬ 
perimental  determination  of  the  capture  cross  sections  of  interface  states 
for  both  electrons  and  holes,  as  well  as  for  the  energy  dependence  of  the 
capture  cross  sections  and  interface  state  densities. 

Consider  an  MOS  structure  as  shown  in  Fig.  1  under  a  bias  voltage  V,  and 
exposed  to  optical  illumination  of  photon  energy  in  excess  of  the  silicon 
energy  gap.  Depending  on  the  illumination  level,  minority  carriers  (holes) 
are  generated  optically  in  the  region  near  the  surface  of  the  semiconductor 
(x  »  0).  Under  the  conditions  where  the  surface  concentration  of  holes  exceeds 
that  of  electrons,  i.e.,  p(0)  >  >n(0),  the  loss  process  described  by  Nicollian 
and  Goetzberger is  dominated  by  the  interaction  between  the  minority  car¬ 
riers  and  interface  states.  Hence  conductance  measurements  of  the  ac  losses 
involved  in  the  charge  exchange  give  information  regarding  those  interface 
states  at  the  energy  where  the  minority-carrier  quasi-Fermi-level  E^  inter¬ 
sects  the  surface.  Here  we  have  assumed  that  the  rate  of  photogeneration  and 
supply  of  minority  carriers  to  the  surface  is  sufficient  so  as  not  to  inter¬ 


fere  with  the  rate  of  charge  exchange  with  interface  states. 


Since  the  basic 
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Figure  1:  Energy-band  diagram  for  a  MOS  structure  on  an  n-tvpe  semicon¬ 
ductor  under  a  given  voltage  bias  V(a)  in  darkness,  (b)  under 
optical  illumination  with  hv>Eg.  E^n  and  Efp  represent  the 
electron  and  hole  Fermi  levels,  Es  the  energy  of  an  interface 
state,  and  Jsc  and  Jsv  are  the  currents  associated  with  charge 
exchange  of  interface  states  with  the  conduction  and  valance 
bands  (as  described  in  Ref.  10). 


phenomena  of  the  conductance  technique  due  to  Nicollian  and  Goetzberger  are 


r 


therefore  preserved,  the  time  constant  and  capture  cross  section  associated 

with  the  minority-carrier  interaction  can  be  determined  as  well  as  the  inter- 

,  .  .  (10) 

face  state  density. 

To  better  understand  the  admittance  of  the  optically  excited  MOS  struc¬ 
ture,  we  now  solve  Poisson's  equation  for  the  case  of  low-level  injection 
caused  by  the  optical  illumination.  A  horizontal  quasi-Fermi-level  E^p  is 
assumed  here  and  this  will  be  justified  in  a  later  section.  Low-level  injec¬ 
tion  implies  that 

nno  =  ni6Xp[(Efn  ‘  Ei)/kT]  *  Nd  (1) 


and 


P*o  -  n.exp[(Ei  -  Efp)/kT],  (2) 

jif  yt 

where  n  and  P  are  the  zero-field  (tfj  =  0)  electron  and  hole  concentrations 
no  no  s 

n^^  is  the  intrinsic  concentration  corresponding  to  an  intrinsic  level  Ei ,  and 

N.  is  the  doping  density.  Poisson's  equation  states  that 
a 


32^  _  C(x) 
3  2  e 


X  s 


(3) 


where 

£  is 
s 


ip  is  the  potential  in  the  semiconductor  at  a  distance  x  from  the  surface 
the  dielectric  constant  of  the  semiconductor,  and 


£(x)  =  q (n 


P*  +  P  -  n  ) 
no  n  n 


(4) 


is  the  space-charge  density.  By  substituting  the  basic  equations 

* 

n  =  n  exp(q<l/kT) ,  ( 5) 

n  no 

P  =  P  exp(-qj)/kT)  ,  (6) 

n  no 


and  using  Eq.  (4)  in  Eq.  (3),  the  second-order  differential  equation  can  be 
solved.  The  results  for  the  surface  electric  field  and  charge  density  are 


IS. I  ‘  Vi'  aV’ 


(7) 


Q  sc  -  ^2EskT^LD)G  «V  PnO/NV  AEf^ 


(8) 


where 


G  (V  PnO/Nd*  AEf> 


(9) 


Ljj  -  (2eskT/q2Nd)1/2,  (10) 

where  L_  is  referred  to  as  the  extrinsic  Debye  length,  P  is  the  zero-field 
tJ  no 

minority-carrier  concentration  in  darkness,  ig  is  the  surface  potential  of  the 
semiconductor,  and 


4Ef  '  Etn  -  Efp  <n> 

corresponds  to  the  separation  or  "splitting”  of  the  majority-and  minority- 

carrier  Fermi  levels  due  to  the  optical  illumination  and  the  resultant  en- 

* 

hancement  of  p(0)  .  If  the  surface  potential  under  ilium  ;.on  is  \ji  ,  the 

s 

change  in  surface  potential  due  to  illumination  is  ti 

*  i  -  i  ;  (12) 

s  s  s 

the  corresponding  change  in  surface  charge  is 
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(13) 


AQ  =  C  AU>  , 
s  ox  s 

since  Ai|j  =  -  AV ^  for  a  given  bias  voltage  V  between  metal  and  semiconductor, 

and  AQ  =  -  Q  ,  where  C  is  the  oxide  capacitance  obtained  from  the  accumula- 
s  ox  ox 

tion  capacitance,  and  V  is  the  potential  difference  across  the  oxide  layer. 
Hence  the  surface  charge  density  under  illumination  is 

<  "  ^s  +  AQs  ■  (14) 

* 

The  value  of  AE.  can  be  determined  graphically  from  the  new  values  of  ty 
r  s 

and  Q  .  Once  this  parameter  is  obtained,  the  surface  concentration  of 

s 

minority  carriers  can  readily  be  found  from 

P(0)  =  Pn()  exp  [  (AEf  -  q**)/kT].  (15) 

Notice  that  the  effect  on  AQg  of  a  change  in  charge  density  in  interface  states 
AQis  is  not  included  in  our  analysis.  It  turns  out,  as  we  shall  discuss  in 
more  detail  in  a  later  section,  the  contribution  of  this  charge  can  be  neglected 
when  the  interface  state  density  is  sufficiently  low,  as  it  is  in  the  present 
samples. 

The  silicon  wafers  used  in  these  investigations  were  (lll)-oriented ,  polished 
n-type  epitaxial  layers  with  film  resistivity  of  7.83  Qcm  grown  on  heavily  doped 
n+  substrates.  These  slices  were  cleaned  in  hot  organic  solvents,  and  in  a 
solution  of  in  the  ratio  of  5:1.  Deionized  water  of  10-MC  purity 

was  used  to  rinse  the  samples,  before  and  after  the  surfaces  were  etched  in 
dilute  hydroflouric  acid  (HF) .  The  samples  were  dried  in  N^  gas,  and  oxidized 
in  dry  at  1200°C  in  the  presence  of  hydrogen  chloride  (HC1).  Hi gh-tempera- 
ture  (1200°C)  annealing  for  10  min  with  high-puritv  gas  was  performed  im¬ 
mediately  after  the  oxidation  process.  The  SiO?  thicknesses  were  measured  bv 
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ellipsometry  prior  to  metallization.  Thin  A1  electrodes  ( t  100  $  thick  and 
1  mm  in  diameter)  were  evaporated  through  an  out-of-contact  mask  at  a  base 
pressure  of  10  ^  Torr.  Postmetalization  annealing  with  forming  gas  (35% 

+  65%  N^)  at  450°C  for  a  duration  of  10  min  was  employed  prior  to  the  re¬ 
moval  of  oxide  layers  and  evaporation  of  Ohmic  contacts  at  the  back  surfaces 
of  the  samples. 

The  frequency  dispersion  of  the  ac  admittance  (capacitance  and  conduc¬ 
tance  vs  dc  bias  voltage)  was  measured  between  32  Hz  and  200  KHz  by  phase- 
sensitive  detection  techniques,  using  a  PAR  Model  124  lock-in  amplifier.  Ac¬ 
curacy  is  checked  by  a  Hewlett  Packard  HP  4274  A  impedance  meter.  A  slow 
ramp  (£  10  mV/s)  was  employed  for  the  dc  bias,  and  data  was  recorded  on  a 
X-Y  recorder,  while  at  the  same  time  being  digitized  and  stored  in  a  DEC-PDP 
11/03  minicomputer.  This  system  has  been  described  in  previous  studies .  (^^-2) 
The  optical  illumination  is  provided  by  a  tungsten  lamp  through  an  optical 
conductor  so  as  to  maintain  the  measuring  sample  at  room  temperature.  The 
received  optical  intensity  is  controlled  by  adjusting  the  distance  d  of  the 
source  above  the  sample  in  order  to  maintain  a  constant  spectral  content  of 

the  illumination.  Resultant  relative  illumination  levels  (proportional  to 
2 

1/d  )  and  the  corresponding  quasi-Fermi-level  separations  are  given  in 
Table  j  for  the  five  illumination  intensities  used  in  this  experiment. 

The  resultant  experimental  dependence  of  interface  state  density  on 
energy  is  shown  in  Fig.  2.  The  full  circles  correspond  to  dark  conductance 
measurements,  the  open  circles  to  conductance  measurements  under  optical  il¬ 
lumination.  N.  vs  E  has  also  been  determined  from  the  low-f requencv  capaci- 
is 

tance  characteristic  in  darkness,  and  excellent  agreement  with  Fig.  2  was 

obtained,  in  justification  of  the  present  measurements.  The  agreement  between 

various  measurements  of  N.  has  also  been  described  in  our  earlier  papers. 

is 
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Figure  2:  Energy  distribution  of  interface  state  density 
obtained  from  dark  and  optically  illuminated 
conductance  measurements . 
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Table  I 


Data  for  the  five  illumination  levels 
used  in  the  present  experiments. 


Illumination 

Levels 

Relative 

Illumination 

Levels 

Ef  ' 
eV 

Dark 

0 

0.000 

1 

1 

0.070 

2 

7.7 

0.137 

3 

25.0 

4 

68.8 

0.245 

5 

156.3 

0.270 

Figure  3  shows  the  observed  dependence  on  energy  of  the  capture  cross 

sections  for  electrons  and  for  holes,  a  and  a  ,  respectively,  for  the  inter- 

n  p 

face  states  of  the  devices.  A  strong  dependence  of  o^  on  energy 

is  obtained,  and  is  found  to  be  several  orders  of  magnitude  larger  than 

a  .  The  three  curves  for  a  correspond  to  illumination  levels 
n  P 

shown  in  Table  I  and  we  note  that  the  same  o  is  not  obtained,  for  a  given 

P 

energy  with  different  illumination  levels,  for  reasons  discussed  below.  The 

results  for  a  and  a  have  been  obtained  from  data  such  as  that  shown  in 
P  n 

Fig.  2  for  various  values  of  applied  bias  voltage  V.  Let  us  first  consider 
the  dependence  of  a  on  energy.  In  the  present  technique,  we  can  obtain 

p 

0^  as  a  function  of  energy  by  using  several  different  combinations  of  illumina 
tion  intensity  and  bias  voltage,  since  both  of  these  affect  the  position  in 


the  energy  gap  of  the  hole  quasi-Fermi  level  E 


fp‘ 

The  capture  cross  sections  o^  of  interface  states  for  minority  carriers 

(holes)  are  shown  to  be  orders  of  magnitude  larger  than  those  for  electrons 

0  indicative  of  an  acceptor-like  coulomb  capture  process,  and  also  to  exhibit 
n 

a  strong  dependence  upon  energy  in  the  lower  half  of  the  silicon  energy  gap, 

for  Al-Si07  (111)  structures  with  oxide  thickness  of  =800  £.  This  apparent 

dependence  on  energy  actually  arises  from  a  dependence  on  the  electric  field 

(13) 

at  the  semiconductor  surface  and  not  from  the  energy  of  the  state  itself. 
Optically  activated  admittance  measurements  are  an  effective  means  of  separat¬ 
ing  these  dependences.  The  physical  theory  underlining  the  field  dependence 
of  the  capture  cross-section  is  being  investigated. 
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F.  STATISTICAL  FLUCTUATIONS  OF  DOPANT  IMPURITIES  IN  ION- IMPLANTED  BIPOLAR 

TRANSISTOR  STRUCTURES  AND  THE  MINIMUM  DEVICE  DIMENSIONS  FOR  VLSI 

SYSTEM  RELIABILITY 

(P.  R.  Prucnal,  W.  Hwang,  H.  C.  Card) 

(Principal  Investigator:  W,  Hwang  (212)  280-3115) 

With  the  trend  towards  increasing  functional  density  of  very  large  scale 
integrated  circuits,  there  are  a  number  of  physical  phenomena  that  will  re¬ 
sult  in  limitations  to  the  successful  operation  of  these  circuits  and  systems. 
An  excellent  review  of  these  phenomena  has  been  presented  by  Keyes^^  and 

several  authors  have  treated  one  or  more  of  the  limiting  mechanisms  with  vary- 

(2)— (13) 

Ing  degrees  of  generality.  Among  the  more  serious  limitations,  i.e. 

those  which  are  encountered  first  as  the  physical  device  and  circuit  dimen¬ 
sions  are  decreased,  are  maximum  power  dissipation,  metal  electromigration, 
edge  leakage  currents,  and  statistical  fluctuations  of  dopant  impurities. 

Other  factors  such  as  punch  through  of  depletion  regions,  hot  carrier  injec¬ 
tion  into  insulators,  and  carrier  velocity  saturation  are  the  result  of  operat¬ 
ing  the  smaller  devices  under  different  conditions  of  electric  field  and  cur¬ 
rent  density  than  their  larger  counterparts,  and  are  not  essential  results 

,  i  a  „  u  ,  (5) (13) 

of  a  truly  scaled  technology. 

We  consider  only  one  of  the  above  limitations  to  the  continual  reduction 
in  device  dimensions,  the  problems  associated  with  the  statistical  nature  of 
the  positions  of  the  dopant  impurities.  Furthermore,  we  restrict  ourselves 
to  the  particular  case  of  ion  implanted  impurities  in  a  bipolar  transistor 
structure,  but  we  feel  that  the  results  obtained  for  this  structure  in  terms 
of  minimum  allowable  dimensions  will  be  approximately  correct  for  other  devices. 

As  the  size  of  any  characteristic  physical  feature  of  the  bipolar  tran¬ 
sistor  (for  example,  the  dimensions  of  the  implanted  emitter  or  base  regions 
in  the  plane  of  the  silicon  suri  .ce,  or  the  width  of  the  base  region  in  the 
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direction  normal  to  the  plane  of  the  surface)  is  reduced,  the  statistical 


fluctuations  among  devices  of  the  number  of  implanted  impurities  in  the  active 

volume  becomes  important  for  two  reasons.  First,  the  impurity  concentrations 

in  a  device  that  has  been  scaled  down  in  size  by  a  factor  k  is  increased  by 

this  same  factor  lc.  (13)  g^nce  an  dimensions  are  scaled  by  k,  the  active 

3 

volume  of  the  device  is  1/k  of  the  larger  device,  so  that  the  number  of  im- 

2 

purities  within  the  device  is  reduced  by  1/k  .  If  this  number  of  impurities 

is  treated  as  Poisson  random  variable,  since  the  standard  deviation  is  then 

the  square  root  of  the  mean,  the  ratio  of  the  standard  deviation  to  the  mean 

2 

is  increased  by  k.  Second,  there  will  be  a  corresponding  increase  of  k  in 
the  number  of  components  (all  of  whose  behavior  influences  that  of  the  system) 
per  chip,  if  the  conservative  estimate  is  made  that  the  chip  size  itself  does 
not  increase.  This  second  consideration  is  of  consequence  because  one  can 
tolerate  a  larger  standard  deviation  with  a  smaller  number  of  devices  for  a 
given  yield  of  the  chip.  In  order  to  calculate  the  yield,  one  must  postulate 
a  tolerable  limit  to  the  minimum  or  maximum  number  of  impurities  per  device 
in  order  that  the  device  (and  hence  the  system  on  the  chip  of  which  it  is  a 
part)  will  function  within  the  specifications.  We  are  primarily  interested 
in  the  minimum  dimensions  for  adequate  yield,  and  since  both  the  statistical 
fluctuations  are  increasing  and  the  fluctuations  we  can  tolerate  are  decreasing 
with  reducing  size,  the  exact  choice  of  allowed  fluctuation  of  any  electrical 

(  14) 

parameter  for  the  device  will  not  greatly  affect  the  results. 

One  final  point  should  be  made  regarding  the  operation  of  all  devices  on 
the  chip.  With  VLSI  systems  it  will  not  be  possible  to  test  all  devices  or 
circuit  functions,  and  "bad”  devices  are  likely  to  escape  unnoticed,  until 
exhaustive  use  of  the  system  detects  the  fault.  Because  of  this,  it  will  be 
necessary  to  be  very  conservative  in  determining  the  allowed  probability  of 


a  faulcy  device  on  the  chip.  One  will  have  to  aim  for  very  high  yields  when 


considering  any  limitation  such  as  this  one  where  the  origin  of  a  failure 
arises  from  a  random  process.  As  an  example,  we  describe  the  effects  of  the 
random  fluctuations  among  devices  in  the  number  of  dopant  impurities  on  the 
behavior  of  ion- implanted  bipolar  transistor  structures. 

Let  us  assume  that  we  can  tolerate  a  100  X%  variation  in  the  current  gain 
3  of  the  transistors  assembled  on  a  given  chip.  This  means  that  the  varia¬ 
tions  in  the  base  Gummel  number  must  be  less  than  100  for  any  device  on 
the  chip,  if  no  additional  variations  in  6  are  introduced  by  fluctuations  in 
the  emitter  Gummel  number.  Fluctuations  in  the  base  Gummel  number  are 

induced  by  variations  in  the  total  number  of  impurities  n  in  the  base  region, 

2 

which  has  mean  value  G^L  »  ,  where  L  is  the  device  dimension  and  is  the 

mean  number  of  base  impurities.  The  probability  .that  a  device  fails  is  then 
equal  to  the  probability  that  n  deviates  from  N,  by  more  than  100  X%.  Using 

D 

the  well-known  Tschebechev  inequality ,  ^ ^  the  probability  th3t  a  given 
device  of  dimension  L  fails  is  then 


P1  -Pr{|„-Nbl  — 

A  N, 


The  probability  P  of  one  or  more  device  failures  on  a  chip  containing 
chip 

N  ,  .  devices  is,  from  Eq.  (1)  and  the  binomial  distribution, 
chip 


P  <  1 
chip  - 


Now  as  the  device  size  L  changes  by  a  factor  k,  Eqs.  (1)  and  (2)  become 


p  <  - ±_ 

*2N, 


,  ■>  2 
A  'k  It 
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p  .  .  <  1  - 
chip  — 


1 


2  2 
kl  N, 


N  U-1 
chip 

k2 


(4) 


In  Eq.  (4)  we  see  that  by  both  decreasing  the  average  number  of  impurities 

2  2 

in  any  given  device  by  k  ,  and  increasing  the  number  of  devices  by  k  ,  the 

probability  of  failure  P^  of  one  or  more  devices  on  a  chip  increases  dra¬ 
matically  with  a  reduction  in  dimensions.  To  remedy  this  problem,  an  N-bit 
parity  check  may  be  used.  In  this  way,  the  probability  of  error  in  an  N-bit 

word  is  approximately  the  binomial  probability  P  that  two  devices  in  N  fails: 

N 


N  2  N-2 

PN  *  f(N-l)P‘  (1-P1)N 


(5) 


Furthermore, 
chip  is 


N 

the  probability  P^^ 


of  one  or  more  N-bit  words  failing  on  a 


(6) 


N 

In  Fig.  1  we  show  the  dependence  of  P^,  pc^p  and  ^chip  on  t*ie  lateral 

dimension  L  of  the  implanted  bipolar  devices,  for  various  word  lengths  N, 

where  A  “  0.20.  As  the  number  of  bits  involved  in  a  parity  check  N  is  reduced, 

the  minimum  dimension  L  for  a  given  probability  of  chip  failure  decreases, 

or  the  probability  of  success  for  a  given  L  increases.  From  results  such  as 

those  of  Fig.  1,  we  are  able  to  plot  the  dependence  of  the  minimum  dimension 

N 

L  upon  N,  for  various  values  of  chip  yield  Y  ■  l-P^^p  (fig.  2) .  As  we  pointed 
out  above,  the  design  must  be  rather  conservative  (design  for  high  yield)  when 
considering  against  a  random  device  fault  which  may  not  easily  be  discovered 
during  testing  or  field  operation.  As  N  decreases,  the  fraction  of  bits 
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PROBABILITY  OF  FAILURE,  (P) 


chip  yields. 


unavailable  for  information  processing  of  course  decreases .  In  the  worst  case 


N  =  2  and  only  every  second  bit  is  used  for  storage  of  information,  so  that 
the  useful  chip  area  is  effectively  halved.  For  this  case,  the  minimum  L  for 
a  yield  of  99%  from  Fig.  2  is  -  1.5  pirn. 

The  random  fluctuations  in  the  number  of  ion  implanted  impurities  in  the 
base  region  of  bipolar  transistor  structures  is  a  serious  limitation  to  the 
continued  reduction  in  device  size  demanded  by  very  large  scale  integration  of 
these  devices.  These  fluctuations  have  a  doubly  stochastic  nature  by  virtue 
of  uncertainties  in  the  implantation  parameters  such  as  the  standard  devia¬ 
tion  of  the  implanted  profiles,  so  that  the  simple  Poisson  statistics  are  not 
sufficient  to  determine  the  minimum  device  size  for  tolerable  chip  yields. 

These  doubly-stochastic  effects  will  be  studied  in  the  future.  This  minimum 
device  size  expressed  as  the  minimum  lateral  dimension  (emitter  size)  is 
expected  to  be  on  the  order  of  1  <  L  <  5  Um  with  a  more  exact  value  given  above 
and  dependent  on  the  chip  organization  for  parity  checks. 
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III.  ENERGY  TRANSFER  AND  RELAXATION  IN  SMALL  POLYATOMIC  MOLECULES 


A.  VIBRATIONAL  ENERGY  TRANSFER  IN  LASER  EXCITED  COF  :  INFRARED  FLUORESCENCE 

FROM  THE  INTERMEDIATE  MODE  v  *  1 

4 

(R.  K.  Bohn,  K.  H.  Casleton,  Y.  V..  C.  Rao ,  G.  W.  Flynn) 

(Principal  Investigator:  G.  W.  Flynn  (212)  280-A162) 

In  a  previous  study^  COF^  was  shown  to  absorb  CO^  infrared  laser 
radiation  in  the  mode  (C-F  symmetric  stretch)  .  The  overtones  of  the 
mode  were  found  to  be  rapidly  populated  by  the  ladder  climbing  process  (1) . 

2COF2  (v2)  -  COF 2  (2v2)  +  COF2(0)  (1) 

Since  the  2v2  level  is  Fermi  mixed  with  the  carbonyl  stretch  state,  strong 

1R  emission  is  observed  from  2v2  -*•  0  at  X  •v  5y  •  For  certain  C02  lines, 

the  population  of  the  2v2  state  can  be  directly  populated  by  hot  band  pumping 

of  the  v2  -*  2v2  transition. ^ ^  Because  of  the  rapid  coupling  of  and 

2v2  via  process  (1),  which  requires  an  average  of  30  gas  kinetic  collisions 

to  reach  steady  state,  fluorescence  from  2v2  simply  reflects  the  changes 

in  population  of  the  v2  state  for  processes  slower  than  (1). 

There  are  several  features  of  vibrational  energy  transfer  in  C0F2 

which  distinguish  it  from  other  polyatomic  molecules,  mostly  derivatives  of 

(2) 

CH^,  that  have  been  investigated  in  detail.  Since  CO ?2  is  completely 
composed  of  heavy  atoms,  most  of  the  fundamental  vibration  frequencies  are 
relatively  low  frequency  rather  than  clustered  into  regions  of  high  and 
low  frequencies.  The  large  moments  of  inertia  of  C0F2  are  expected  to 
reduce  the  importance  of  rotational  motions  in  the  vibrational  energy 

(3) 

transfer  processes  and  to  allow  for  direct  collisional  energy  transfer 
between  the  fundamentals  (v  =  1  levels)  of  most  modes.  Relatively  strong 
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fluorescence  may  be  observed  from  at  least  four  different  states  (2v^, 


i 


at  5u ,  v,  at  8u,  v,  at  13u,  and  v.  and/or  v.  at  (16-17u).  Thus,  the 
4o  3  5 

populations  of  these  states  may  be  independently  monitored  as  a  function 

of  time  after  pulsed  laser  pumping.  COF^  is  a  tetra-atomic  molecule  with 

relatively  dense  energy  levels  whose  energy  transfer  properties  should  more 

closely  approximate  those  of  larger  molecules.  Stretch-stretch  vibrational 

energy  transfer  occurs  rapidly  in  CH^ ,  11+3  collisions,  for  the  symmetric- 

-1  (4) 

asymmetric  stretch  transfer  which  is  endothermic  by  about  100  cm 

Since  COF^  contains  no  hydrogen  atom,  determination  of  stretch-stretch  energy 

transfer  rates  will  indicate  whether  heavy  atom  molecules  behave  differently. 

The  8p  fluorescence  from  the  state  of  COF^  is  particularly  interest¬ 
ing  because  its  growth  after  pulsed  laser  pumping  is  characterized  by  two 
exponential  rises.  This  characteristic  is  rather  unusual  and  has  not  been 
frequently  observed  in  IR  fluorescence  studies  of  vibrational  energy  transfer. 
The  consequences  of  this  double  rise  for  understanding  energy  transfer  in 
COF2  are  the  major  subject  of  this  study. 

The  experimental  apparatus  for  laser-induced  fluorescence  was  composed 
of  a  Q-switched  CO^  laser  previously  described  in  detail.  ^  In  most  of 
the  experiments  the  P(22)  or  the  R(14)  lines  of  the  10. 6u  band  were  used 
for  excitation.  Laser  pulses  were  repeated  at  200  Hz  with  vl  mJ  pulse  ^ 
and  a  width  of  M.  ps  (FWHM) .  Typically,  >80%  of  the  laser  power  was  carried 
by  the  exciting  line.  The  stainless  steel  fluorescence  cell  had  an  outgas 
rate  of  <10  mtorr  hr  ^ .  Pressure  was  measured  with  a  corrosion-resistant 
capacitance  manometer  (MKS  Baratron). 

Carbonyl  fluoride  samples  were  obtained  from  both  PCR  Research  Chemicals 
(Gainesville,  FL)  and  Svnthatron  Corp.  (Parsippanv,  NJ)  which  contained 


^10%  CC^  impurity.  Samples  were  purified  gas  chromatographically  using  a 

1.5  ft  x  1/A  in.  stainless  steel  column  filled  with  Carbosieve  S  (Supelco, 

Inc.,  Bellefonte,  PA).  Another  column  apparently  identical  to  that  used 

in  an  earlier  purification^^  did  not  separate  CO^  from  COF^ . 

Fluorescence  perpendicular  to  the  laser  beam  was  observed  at  8y  using 

a  7.5  -  8.5p  band  pass  filter  and  a  Au:Ge  (77K)  detector.  At  5p  fluorescence 

was  monitored  with  a  MgF^  exit  window,  a  A.6p  long  pass  filter  and  either 

the  Au:Ge  detector  or  an  InSb  (77K)  detector.  The  fluorescence  signal  was 

observed  with  a  transient  digitizer/digital  signal  averager  combination  and 

stored  on  floppy  disks  for  subsequent  analysis.  Data  were  analyzed  using 

a  modification  of  a  non-linear  least  squares  program. ^  Fluorescence  data 

n  -Fjt 

were  fit  to  the  function,  Z  A  e  +  B  where  B  is  a  background,  A.  is 

i=l  1  1 

the  amplitude  of  the  ith  term,  r^  is  the  corresponding  rate,  and  n  can  be 

chosen  as  one  or  two.  A.,  r.  and  B  can  be  varied  but  B  was  held  fixed  at 

1  1 

the  pre-laser  pulse  value  in  our  analyses.  Fluorescence  curves  composed 

of  more  than  two  exponential  terms  were  analyzed  by  fitting  the  slowest 

exponential  terms,  subtracting  the  calculated  fits,  and  fitting  the  difference 

curves  to  further  sums  of  exponential  terms. 

An  abbreviated  vibrational  energy  level  diagram  for  COF^  is  shown  in 

(8) 

Fig.  1.  The  CC>2  laser  pumps  ground  state  molecules  into  the  level. 

Emission  from  the  level  occurs  at  8p.  This  fluorescence,  shown  in  Fig.  2, 

is  describable  by  two  exponential  rise  terms  and  two  exponential  fall  terms. 

One  of  these  is  a  very  slow  but  large  amplitude  decay  which  is  due  to  the 

flow  of  heat  from  the  center  of  the  fluorescence  cell  to  its  walls  at  a 

(9) 

rate  controlled  by  the  thermal  conductivity  of  the  gas.  The  rise  terms 
are  characterized  by  the  rate  constants  156  +  18  ms  "*'torr  ^  and  29  +  5  ms  * 
torr  ^ .  The  reported  uncertainties  are  two  standard  deviations.  The 
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amplitude  of  the  29  ms  torr  process  is  about  one-half  that  of  the 
156  ms  "'"torr  ‘  process.  The  primary  fall  of  the  fluorescence  is  described 
by  the  rate  constant  5.6  +  0.2  ms  ^torr  The  thermal  conductivity  rate, 
which  was  not  measured,  was  of  the  order  of  100  sec  ^ . 

An  earlier  study^^  demonstrated  that  the  overtone  of  the  pumped 
state,  2 v^,  Fermi  mixed  with  v^,  is  populated  rapidly  by  the  ladder  climbing 
process  (1).  Therefore,  the  5u  fluorescence  from  this  state  reflects 
population  changes  in  the  state  since  2v2  and  v2  are  so  rapidly  coupled. 
The  decay  of  the 5p  fluorescence  is  describable  by  three  exponential  terms, 
one  of  which  is  a  very  slow  but  small  amplitude  thermal  decay.  The  rate 
constants  for  the  two  major  decays  are  34  +  2  ms  ^  torr  ^  and  7.4  +  1.2  ms  ^ 
torr  ^ .  The  amplitude  of  the  faster  process,  which  was  reported  in  the 
earlier  study, ^  is  about  7  times  that  of  the  slower  process. 

There  seems  little  doubt  that  the  measured  rate  constants  of  5.6  +  0.2 
and  7.4  +  1.2  ms  ^torr  1  in  the  8y  and  5u  fluorescence,  respectively,  simply 
reflect  the  same  kinetic  eigenvalue  observed  in  two  different  states. 
Furthermore,  this  is  the  rate  most  closely  associated  with  "vibration- 
translation/rotation"  (V-T/R)  relaxation  of  the  vibrational  manifold  since 
it  is  the  slowest  rate  and  one  that  is  immediately  followed  by  thermal 
heating  of  the  system  (which  is  suggestive  of  overall  equilibration  of 
all  degrees  of  freedom  at  a  temperature  above  ambient).  In  C0F2  some 
caution  must  be  exercised  in  ascribing  this  rate  to  V-T/R  relaxation  of 
just  the  very  lowest  levels  (v^,  v^)-  Because  of  the  relatively  balanced 
coupling  between  different  modes  and  the  moderately  slow  intermode  energy 
transfer  for  this  molecule,  the  slowest  or  "V-T/R"  eigenvalue  is  expected 
to  contain  substantial  contributions  from  one  or  more  intermode  vibration- 
vibration  (V-V)  processes . 
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The  rates  of  29  +  5  and  34  +  2  ms  ^"torr  ^  observed  in  the  8y  and  5p 
fluorescence,  respectively,  also  appear  to  reflect  the  same  kinetic 
eigenvalue  appearing  in  two  different  states.  The  8p  rate  is  particularly 
difficult  to  measure  accurately  because  it  is  the  third  eigenvalue  stripped 
off  the  fluorescence  curve  and  has  only  about  half  the  total  signal  amplitude 
In  addition  the  8y  fluorescence  signals,  though  good  for  the  infrared,  are 
noticeably  more  noisy  than  those  at  5y.  Thus  we  assume  that  these  two 
rates  are  the  same  kinetic  eigenvalue  appearing  with  different  amplitudes 
in  the  two  different  fluorescences.  There  are  other  compelling  reasons, 
as  noted  below,  for  believing  that  the  major  decay  rate  for  \>2  l°ss  must 
appear  in  the  rise  of  v^. 

In  proposing  a  mechanism  to  explain  the  observed  data  for  these  two 
fluorescing  states,  the  following  points  are  very  important:  (1)  must 
be  filled  from  either  the  pumped  state  or  the  2\>^  overtone,  Fermi 
mixed  set.  Given  the  energy  gap  between  2v2>  and  plus  the  relatively 
small  populations  in  these  high  "overtone"  states,  filling  from  is  by 
far  the  most  reasonable  mechanism.  (2)  Except  for  the  "V-T/R"  eigenvalue, 
the  2v2>  states  decay  only  at  a  single  rate  of  approximately  34  ms  "'‘torr  ^ 
As  noted  above,  due  to  the  rapid  coupling  of  v2  to  2v2  (and  v^)  via  (1), 
this  decay  surely  represents  flow  of  energy  (population)  out  of  the  v2 
level  and  into  the  other  vibrational  modes.  The  5p  decays  merely  reflect 
the  behavior  of  the  lower  level  (3)  Despite  the  solid  conclusions  of 

the  above  two  comments,  appears  initially  to  fill  faster  (156  ms  '''torr 
than  v2  decays  (34  ms  ^torr  ’*’)  .  Since  the  v2  reservoir  must  be  the  source 
of  energy  for  v^,  the  fast  filling  of  is  extremely  anomalous. 

The  above  features  taken  together  suggest  that  the  fast  rise  (156  ms  ^ 
torr  of  actually  corresponds  to  the  loss  of  population  from  even 


though  the  156  ms  torr  rate  appears  in  the  fluorescence  rise.'  In 
addition  the  slow  decay  of  (34  ms  ^torr  also  appears  in  the  rise 
as  would  be  expected  if  filled  from  v^.  As  will  be  shown  below  in  the 
discussion  of  kinetics,  a  state  being  filled  slowly  (e.g.  34  ms  ^torr 
but  emptied  rapidly  (e.g.  156  ms  ^torr  can  exhibit  exactly  this  behavior. 
While  such  a  situation  is  kinetically  unusual,  it  is  by  no  means  unique. 

For  example,  when  the  stretching  modes  of  SO^  are  laser  pumped,  they  in 
turn  fill  the  bending  mode  more  slowly  than  it  is  emptied  by  V-T/R 
relaxation. 

Two  exponential  rise  terms  in  the  population  of  a  state,  in  this 
instance,  could  also  occur  if  the  state  is  filled  from  two  different 
reservoirs  at  two  different  rates.  Since  COF^  is  initially  excited  to  only 
a  single  state,  it  is  conceivable  that  energy  is  transferred  in  such 

a  way  that  two  population  reservoirs  develop  which  each  fill  the  state 
at  different  rates.  However,  the  energy  levels  near  the  state  are  sparse 
enough  that  this  mechanism  seems  unlikely.  Furthermore,  one  of  the  popula¬ 
tion  reservoirs  must  ultimately  have  its  source  in  the  state,  which 
decays  with  a  rate  constant  of  34  ms  ^torr  ^ .  Thus  there  seems  to  be  no 
way  to  generate  another  population  reservoir  so  quickly  that  it  fills  the 

v,  state  with  a  rate  constant  of  156  ms  ^torr 
4 

(12) 

In  1910,  Lowry  and  John  discussed  the  kinetics  of  successive, 
reversible  first-order  reactions  (2) 

k12  k23 

A:  A2  A3  (2) 

k21  k32 
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Lowry  and  John  emphasized  that  if  k^2  >  k^,  ^•e*  t^le  rate  constant  for 
filling  from  the  initial  state  is  greater  than  the  rate  constant  for 
filling  A^  from  the  final  state  A^,  then  the  population  of  the  intermediate 
species,  A 2>  passes  through  a  maximum.  What  they  did  not  emphasize  is 
that  if  k^9  >  then  the  population  of  A^  rises  asymptotically  to  an 

equilibrium  value  and  that  the  rise  is  a  sum  of  two  exponential  terms. 

Thus,  if  one  identifies  the  species  A^  with  the  initially  pumped  state,  v^, 
the  intermediate  species  with  the  fluorescing  state,  v^,  and  the  final 
species  A^  with  a  state  or  states  such  as  +  v^,  2v^,  and/or  2v^  (see 
Fig.  1)  then  the  rise  of  the  state  is  expected  to  be  described  by  two 
exponential  terms  if  k^2  >  k^2«  ^'^■s  mechanism  is  also  completely  compatible 

with  the  rate  constant  34  +  2  ms  ^orr  \  equaling  the  slower  rise  of  the 
intermediate  state,  29  +  5  ms  ^torr  \  Thus,  we  conclude  that  the  mechanism 
in  which  the  state  is  an  intermediate  between  the  initially  prepared 
v2  state,  and  a  final  state  is  most  compatible  with  the  observation. 
Reexpressing  mechanism  (2)  for  COF^ 

COF  (v,)  +  COF  +  278  cm"1  (3) 

2  4  2 

COF2  (x)  +  COF2  +  (4  to  79  cm-1)  (4) 

where  state  x  probably  refers  to  either  2v^  (AE  “4  cm  1) ,  2v^  (AE  *  79  cm  1) 
or  v3+v3  (AE  »  43  cm  1) .  Although  the  Lowry- John  model  also  predicts  that 
the  initially  pumped  state  (v2)  will  decay  via  two  exponential  terms,  the 
amplitude  of  the  faster  decay  would  be  so  small  in  the  present  case  as  to 
be  experimentally  unobservable. 


cof2  (v2)  +  cof2 


24 


42 


cof2  (v4)  +  cof2 


4x 


'  x4 
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The  relaxation  mechanism  (3,4)  is  incomplete  since  it  does  not  recount 

for  the  other  fluorescences  observed  in  COF^  nor  does  it  include  V-T/R 

relaxation  which  allows  all  of  the  state  populations  to  return  to  equilibrium 

values.  It  is  this  process  which  is  largely  involved  in  the  'v  6  ms  ^torr  1 

exponential  decays  observed  in  both  the  8u  and  5u  fluorescence.  If  w^ 

choose  for  the  elementarv  rate  constants  the  values  k„,  =  20  ms  ^torr  ^ 

24 

(800  collisions)  and  k^  =  50  ms  "''torr  ^  (325  collisions),  we  obtain 

160  and  36  ms  ^torr  ^  as  the  fluorescence  rise  rates  for  the  v,  state  with 

4 

the  amplitude  of  the  faster  term  twice  that  of  the  slower  term.  The 

effects  of  other  processes  which  must  occur  in  the  complete  energy  transfer 

pathway  will  modify  the  values  of  the  elementary  rate  constants  in  order 
to  maintain  agreement  with  the  observed  rate  parameters.  For  example, 
deactivation  of  COF^x)  states  must  occur  at  a  near  gas  kinetic  rate  via 
processes  such  as 

COF2(2v5)  +  C0F2(O)  +  2COF2(v5) 

C0F2(?v3)  +  COF2(0)  t  2COF2(v3)  (5) 

COF2(v3  +  v5)  +  COF2(0)  £  C0F2(v3)  +  COF^v^ 

In  addition  there  is  some  evidence  in  other  fluorescence  emissions  that 

the  relaxation  of  v2  via 

COF  (v  )  +  C0F.  2  C0F  (v£)  +  C0F,  +  199  cm"1  (6) 

2  2  2  2  0  2 

may  be  competitive  with  deactivation  through  (process  3). 
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Conclusions 


(1)  The  asymmetric  stretch  mode  of  COF^  fills  from  the  symmetric 
stretch  level  via  a  relatively  inefficient  collisional  mechanism  which  is 
translationally  endothermic. 

(2)  Efficient  loss  of  population  from  occurs  rapidly  through  the 

overtone  combination  levels  2v^»  2v^,  which  are  nearly  resonant 

with  v, . 

4 

(3)  The  slow  filling  of  coupled  with  its  rapid  emptying  creates  a 
kinetic  anomally  in  which  the  filling  and  emptying  rates  both  appear  in  the 
rise  of  fluorescence  from  this  state. 

(4)  The  initially  pumped  mode  of  COF^  appears  to  decay  via  two 

competing  mechanisms  (the  endothermic  path  through  and  an  exothermic 

path  through  v,). 

o 

(5)  The  rate  of  collision  induced  intermode  energy  transfer  in  COF^ 
is  quite  slew  (£  500  collisions)  for  a  four  atom  polyatomic.  The  features 
which  seem  to  lead  to  this  situation  are  the  relatively  wide  spacing  of 
vibrational  mode  frequencies  and  the  presence  of  all  heavy  atoms  in  the 
molecule.  Thus  the  amplitude  of  vibrational  motions  and  the  ability  of 
rotations  to  take  up  energy  are  reduced  compared  to  hydrogen  or  deuterium 
containing  species. 

(6)  The  overall  vibration-translation/rotation  relaxation  process 
for  pure  COF^  requires  approximately  2500  collisions  on  the  average. 

*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  NSF-CHE77-11384  and  by  the  Department  of  Energy  under  Contract 
DE-AS02  -  ER-78-S-0 2-4940 . 
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B.  TEMPERATURE  DEPENDENCE  OF  SO?  ENERGY  TRANSFER  RATES* 

(J.  L.  Ahl,  G.  W.  Flynn) 

(Principal  Investigator:  G.  W.  Flynn  (212)  280-4162) 

The  development  of  new  molecular  lasers  is  dependent  upon  our  under¬ 
standing  of  the  V-V  and  V-R,T  relaxation  processes  occurring  in  vibrationallv 
excited  gases.  There  have  been  many  measurements  of  the  relaxation  rates 
of  polyatomic  molecules  at  room  temperature  and  theories  to  explain  these 
results,  however  a  proper  test  of  our  understanding  of  these  processes 
requires  a  knowledge  of  the  temperature  dependence  of  the  relaxation  rates. 
Though  a  number  of  such  studies  have  been  conducted  using  the  ultrasonic 
absorption/dispersion  technique,  the  proper  interpretation  of  the  raw  data 
is  open  to  question. ^ The  more  reliable  technique  of  laser  induced 
fluorescence  has  been  used  to  measure  temperature  dependences  only  in  a  few 
cases,  and  then  on  hydrogen  containing  molecules. 

We  are  currently  studying  the  vibrational  relaxation  of  SO^  over  the 
temperature  range  210  K  to  500  K.  Sulphur  dioxide  has  three  infrared  active 
vibrational  modes,  the  stretch  at  1151  cm  \  the  bend  at  518  cm  , 
and  the  stretch  at  1361  cm  .  Using  a  Q-switched  CO^  laser,  we  excite 
the  band  and  observe  the  fluorescence  decay  of  the  band.  Previous 
room  temperature  studies^^  show  a  rapid  equilibration  of  the  and 

bands,  followed  by  a  slower  equilibration  with  the  band.  This  suggests 
the  possibility  of  an  optically  pumped  laser  operating  on  either  the  -*• 
or  -*•  Vj  transition.  Observation  of  the  temperature  dependences  of  these 
rates  will  indicate  the  proper  operating  parameters  for  such  a  laser. 

In  addition  to  studying  the  relaxations  of  pure  S0o,  we  are  also 
studving  the  temperature  dependence  of  the  rare  gas  induced  relaxation  of 
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SO,,.  The  experimental  results  can  he  compared  to  a  number  of  semiclassical 
and  quantum  calculations  on  this  system. 

In  our  preliminary  studies,  we  have  obtained  results  on  the  pure  SO^ 
system  over  a  temperature  range  from  218  K  to  365  K.  The  raw  data  has  been 
analyzed  assuming  a  kinetic  system  consisting  of  only  two  rates,  the  V-V 
equilibration  of  and  v^,  and  the  (slower)  equilibration  of  the  v^, 
pair  with  the  2v^  state.  The  2v^  state  is  assumed  to  relax  very  rapidly 
to  the  state,  and  thus  is  not  significantly  populated.  From  the  observed 
decay  rates,  we  are  able  to  obtain  the  above  noted  microscopic  rate  constants. 
The  temperature  dependence  of  the  slow  rate  corresponding  to  the  v^,  to 
2v^  equilibration,  agrees  with  an  SSH  theory  calculation  by  Rao  .  ^  However, 
the  fast  rate  constant,  corresponding  to  the  to  equilibration  shows 
non-SSH  behavior  at  low  temperatures.  Preliminary  results  suggest  that  the 
rate  constant  reaches  a  minimum  at  approximately  room  temperature,  and 
increases  at  lower  temperatures.  Further  experiments  are  underway  to  check 
these  results  and  to  determine  the  rare  gas  dependence  of  these  rate 
constants . 

*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  NSF-CHE77-11384 . 
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c. 


LASER  EXCITED  INFRARED  FLUORESCENCE  IN  OXALYL  FLUORIDE 


(T.  H.  Allik,  C.  W.  Flynn) 

(Principal  Investigator:  G,  W.  Flynn  (212)  280-4162) 

Time-resolved  laser  excited  fluorescence  experiments  have  appeared 

in  the  literatures  since  the  mid  1960's.^^  Although  numerous  articles 

have  been  written  on  the  topic  of  vibrational  relaxation  in  simple 

polyatomic  molecules,  little  is  understood  concerning  the  energy  transfer 

mechanism  of  molecules  with  "hundreds"  of  vibrational  states.  We  have 

studied  one  such  case,  the  molecule  oxalyl  fluoride,  trans-C^O^F,, . 

Oxalyl  fluoride  has  twelve  fundamental  frequencies  and  the  number 

-1  (2) 

of  vibrational  states  alone  below  2500  cm  is  >  5000.  The  molecule 

-1  (3) 

has  a  low  lying  torsional  mode  v^,  whose  frequency  is  127  cm  .  Ab 

initio  calculations  for  the  barrier  to  internal  rotation  in  C^O^F,  suggest 
-1  (4) 

a  value  of  1830  cm  .  An  interesting  situation  occurs  between  approxi¬ 

mately  2000  and  8000  cm  ^ .  The  mode  is  in  the  "quasicontinuum"  whereas 
most  of  the  other  fundamentals  are  in  discrete  states.  This  is  not  the 
case  for  "simple"  polyatomics! 

The  antisymmetric  C-F  stretching  mode  of  oxalyl  fluoride  was 
excited  with  a  Q-switch  CC^  laser  operating  at  the  R  branch  of  the  9.6u 
band  (energy  =  3/4  mj/pulse) .  Laser  induced  fluorescence  was  observed  from: 
the  antisymmetric  carbonyl  stretching  mode  at  5.3p,  the  combination 
state  at  8.0u  and  states  between  ll-14u.  Figure  1  shows  the 

vibrational  energy  states  and  the  density  of  vibrational  states  for  C0FC0F. 
In  all  three  cases  the  fluorescence  was  observed  to  be  quite  similar  with 
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h600 


a  fast  rise  followed  by  a  double  exponential  decay.  Figures  2,  3  and  4 
show  the  fluorescence  from  the  5  and  8  u  states.  The  fast  rise  is  within 
the  laser  pulse  indicating  an  equilibration  of  energy  in  <  2  gas  kinetic 
collisions.  The  slow  decay  is  attributed  to  the  decay  of  blackboav 
radiation  from  "hot"  thermalized  gas  in  the  center  of  the  cell.  The 
fast  decay  rate  for  the  5.3u  state  was  analyzed  by  fitting  semilog  plots 
of  fluorescence  intensity  versus  time.  These  rates  were  then  measured 
as  a  function  of  ^2°2^2  Pressure>  over  a  range  from  50  to  250  millitorr. 
Figure  5  shows  a  Stern-Vollmer  plot  of  rate  vs.  pressure  which  provides  the 
desired  rate  constant.  The  other  states  (8y  and  ll-14p)  show  the  same  order 
of  magnitude  decay  rate. 

Since  we  are  not  looking  at  fluorescence  from  the  pump  mode,  the  very 
fast  rise  time  indicates  that  a  rapid  V-V  equilibration  is  occurring.  This 
fast  V-V  process  can  be  either  collisional  or  collisionless  in  nature.  If 
it  is  collisional,  the  molecules  are  undergoing  an  extremely  efficient 
energy  exchange  process.  If  it  is  collisionless,  oxalvl  fluoride  is 
exchanging  quanta  between  the  modes  intramolecularly .' 

The  mechanism  for  the  V-V  energy  transfer  is  however  not  yet  clear. 

Two  possible  pathways  may  be  occurring: 

I.  Fluorescence  from  a  Combination  State 

Excitation  +  hv  -*■  C2O2^2^V10^  ^ 


Ladder 

Climbing 


2C202F2(v10>  ;  C2O2F2(2v10)  *  C202F2<0) 


(2) 


FLUORESCE 

:.4  uSr  " 

K  SHO 


0.057  TORR 
C  FULL  SCALE 
S 


FIGURE  2 


5.3u  FLUORESCENCE 
PRESSURE  0.216  TORR 


1 


0.48P  SEC  FULL  SCALE 
44  K  SHOTS 


FIGURE  3 
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8. On  FLUORES 


PRE 

40. 


URE 


64  K  SHO 


FIGURE  4 


V-V  Equilibration 
(Collisional  (M) 
or  Colllsionless) 


C2°2F2(2v10)  +  (M)  *■  C2°2F2(V9  +  V  +  (M)  +  34  CTD 


Fluorescence 


C2°2F2(V9  +  V  ""  C2°2F2(v5)  +  hv 


1 

(3) 

(A) 


or 


Hot  Band 
Excitation 


V-V 

Equilibration 


Fluorescence 


II.  Hot  Band  Absorption 

C2°2F2(V  +  hv  ^  C2°2F2(V3  +  V10)  (5) 

C2°2F2(V3  +  V10)  +  (M)  *  C2°2F2(v9)  +  (M)  +  37  cm~1 

(6) 

C2°2F2CV  "■  C2°2F2(0)  +  hv  (7) 


Nevertheless,  the  cold  gas  filter  technique^3^  should  be  able  to  differentiate 
between  the  two.  These  experiments  are  currently  underway  and  should  yield 
interesting  results. 

The  fast  decay,  V-T  rate,  which  corresponds  to  approximately  8  gas 
kinetic  collisions,  can  be  attributed  to  the  presence  of  the  large  number 
of  low  lying  modes.  These  modes  provide  a  near  resonant  path  for  the  initial 
excitation  to  be  converted  into  translation,  with  the  maximum  defect  of 
127  cm'1. 

As  seen  in  the  case  of  acetyle  fluoride, ^  the  relaxation  of  molecules 
with  a  dense  vibrational  manifold  having  states  within  kT  of  the  ground 
state  appears  to  be  extremely  fast.  These  systems,  where  discrete  states 
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exchange  energy  with  states  in  the  quasicontinuum,  provide  good  models  for 
relaxation  of  more  complex  molecules. 


(1)  For  a  brief  history,  see  J.  T.  Yardley,  "Introduction  to  Molecular 
Energy  Transfer"  (Academic  Press,  New  York,  1980),  Chapter  3. 

(2)  Assuming  a  harmonic  oscillator  with  a  3  cm  ^  anharmonicity  in  v., . 

(3)  J.  L.  Hencher  and  G.  W.  King,  J.  Molec.  Spectrosc.  IJS*  168  (1965). 

(4)  J.  Tyrrell,  J.  Am.  Chem.  Soc.  98_,  5456  (1976), 

(5)  R.  McNair,  S.  Fulghum,  G.  Flynn,  M.  J.  Feld,  and  B.  Feldman,  Chem. 
Phys.  Lett.  ^8,  241  (1977). 

(6)  Columbia  Radiation  Laboratory  Progress  Report  It 29,  March  .,1,  1979, 
pp.  65-66. 
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D.  COLLISION  ASSISTED  LASER  EXCITATION  OF  SFfe  IN  A  SUPERSONIC  NOZZLE* 

(M.  I.  Lester,  D.  R.  Coulter,  L.  M.  Casson,  R.  B.  Bernstein,  G.  W.  Flynn) 

(Principal  Investigator:  G.  W.  Flynn  (212)  280-4162) 

Previous  experiments  have  shown ^  that  internal  energy  can  be  deposited 
in  a  beam  of  SF,  by  irradiation  with  a  low-power  cw  CO.  laser  in  the  "colli- 

D  L 

(2)  (3) 

sional  region"  at  the  nozzle  exit.  '  Greatly  enhanced  excitation  has 

(A) 

now  been  achieved  by  irradiating  inside  the  nozzle.  Both  arrangements 

take  advantage  of  molecular  collisions  to  overcome  inherent  weak  absorption 
in  polyatomic  molecule  beams.  Collisions  assist  by  broadening  absorption 
lines,  promoting  rotational  hole  filling,  and  providing  up-pumping  through 
(V-V)  energy  transfer.  Following  irradiation,  excited  beam  molecules  move 
rapidly  into  the  collision-free  region,  freezing  much  of  the  absorbed  energy 
in  vibrational  modes. 

Collisions,  which  are  an  inherent  part  of  the  excitation  and  supersonic 
expansion  processes,  may  also  lead  to  some  relaxation  of  the  absorbed  energy. 
The  partitioning  of  the  absorbed  energy  in  the  molecular  beam  between  internal 
and  translational  degrees  of  freedom  can  be  examined  as  a  function  of  col¬ 
lision  rate.  Further,  the  laser  experiments  can  be  compared  to  the  expan¬ 
sion  of  a  thermally  equilibrated  system  prepared  by  bulk  heating  of  the 
nozzle  source.  The  extent  to  which  laser  excitation  with  collisions  assist¬ 
ing  in  the  absorption  process  produces  a  non- equilibrium  distribution  which 
differs  from  bulk  heating  can  thus  be  determined. 

The  previous  apparatus^  has  been  provided  with  an  infrared-transparent, 
capillary  nozzle  beam  source.  ^  The  nozzle  is  a  triangular  cross-section, 

10  mm  long  capillary  formed  by  clamping  a  flat  ZnSe  plate  on  a  polished  block 


93 


r 


of  stainless  steel  containing  a  "V"  shaped  channel  (0.2  mm  deep  x  0.3  mm  wide). 

The  flowing  gas  within  the  capillary  is  irradiated  ("double  pass")  by  a  cw 

2 

CO^  laser  at  10.6  Um  [P(16),  5.5-8. 5W;  area  £a.  1  cm  ]  focused  with  a  0.4m  f .£. 

lens  at  any  given  point  along  the  capillary  upstream  or  downstream  of  the  exit. 

The  electron  bombardment  ionizer-quadrupole  mass  filter  monitors  the  number 

density  n  of  SF.  in  the  beam  and  provides  time-of-f light  (TOF)  measurements 
b 

of  velocity  distributions.  The  energy  flux  (power,  W)  carried  by  the  beam  is 
measured  with  the  pumped  liquid  He-cooled  bolometer. For  many  different 
irradiation  positions,  mass  spectrometer  and  bolometer  measurements  of  the 
mechanically  chopped  (25  Hz)  molecular  beam  were  recorded  (via  lock-in  detec¬ 
tion)  with  and  without  cw  laser  excitation.  TOF  measurements  were  made  at 
fewer  positions  of  laser  irradiation;  thus  velocity  distributions  were  inter¬ 
polated  at  intermediate  positions.  Separate  bolometer  experiments  have  shown 
that  laser  scatter  along  the  capillary  axis  is  undetectable  except  when  the 
laser  irradiates  the  exit  surface  of  the  nozzle. 

For  the  heated  nozzle  experiments,  the  ZnSe  plate  of  the  infrared- 
transparent  capillary  nozzle  beam  source  is  replaced  with  an  identical  plate 
of  highly  polished  copper.  The  entire  capillary  nozzle  assembly  (copper 
plate,  stainless  steel  block  and  clamp)  is  resistively  heated  using  heating 
tape.  The  nozzle  temperature  is  monitored  with  a  chromel-alumel  thermocouple 
which  is  pos.  .  .oned  within  the  copper  plate.  Mass  filter,  bolometer  and 
TOF  measurements  were  recorded  as  a  function  of  nozzle  temperature  in  the 
range  from  290  to  550  K.  Since  heating  and  cooling  of  the  nozzle  may  take 
several  hours,  accurate  measurements  require  that  the  AC  components  of  the 
bolometer  and  mass  filter  be  stable  on  this  time-scale. 

Fig.  1  shows  the  following  results  from  the  laser  excitation  experiments 
plotted  vs.  laser  irradiation  position:  the  relative  change  in  bolometer 
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Figure  1:  Plotted  as  a  function  of  laser  irradiation  position  along  the 
infrared-transparent  nozzle  are:  (a)  the  relative  change  in 
bolometer  (Wjyw0 )  anc*  (b)  ®ass  filter  (n^/r^)  signals  (open 
circles:  bolometer  data  corrected  for  "laser  scatter''^)); 

(c)  the  stream  velocity  vg  and  (d)  the  width  parameter  as  from 
TOF  velocity  analysis;  (e)  the  calculated  average  total  (T) , 
internal  (I)  and  kinetic  (K)  energy  per  beam  molecule  (open 
points  indicate  use  of  interpolated  velocity  distributions) ; 
and  (f)  the  ratio  of  the  increment  in  I  (i.e.,  AE*nt)  to  that 
in  T  (AEtot)  due  to  laser  excitation. 


(H  /H  )  and  mass  filter  (n^/n^)  signals  upon  laser  excitation,  and  the  TOF 

results  (v  ,ct  ) .  (All  of  the  results  shown  are  for  a  pressure  of  205  Torr  of 
s  s 

SF^  in  the  oven.)  Irradiating  the  beam  downstream  of  the  exit  yielded  small 
effects  similar  to  those  observed  with  the  original  nozzle . ^ ^ ^  Ir¬ 
radiation  inside  the  nozzle  induced  a  large  enhancement  (ca .  sixfold)  in 
bolometer  output  accompanied  by  an  increase  in  velocity  of  the  beam  molecules 
and  an  increase  in  their  number  density.  These  effects  are  seen  to  be  maxi¬ 
mized  ca.  1  mm  upstream  of  the  exit.  Laser  irradiation  of  pure  Ar  or  of  SF 
-  o 

at  wavelengths  far  from  resonance  yielded  no  enhancement  of  bolometer  signal. 

Under  all  irradiation  conditions  with  or  without  SF,  in  the  capillary  the 

b 

temperature  rise  of  the  stainless  steel  nozzle  was  merely  a  few  degrees,  even 
over  the  course  of  many  minutes  of  continuous  irradiation. 

The  total,  translational  and  internal  energy  per  molecule  in  the  beam, 
calculated  from  the  laser  excitation  data  as  in  Ref.  (1),  are  also  shown  in 
Fig.  1.  Irradiation  inside  the  nozzle  yields  an  increase  in  both  the  average 
kinetic  energy  (K)  and  the  average  internal  energy  (I) .  At  the  optimal  ir¬ 
radiation  position  the  increment  in  average  total  energy  (T)  due  to  laser 
excitation  in _ca.  9.5  x  10  ^  J  molec  ^  (4760  cm  \  0.6  eV) ,  equivalent  to 
ca.  five  10.6  um  photons.  Of  this  absorbed  energy  _ca.  80%  is  retained  as 
internal  excitation. 

Normal  heating  experiments  yielded  similar  qualitative  trends  as  observed 
with  laser  excitation.  As  the  nozzle  was  heated,  both  the  velocity  of  the 
beam  molecules  and  their  number  density  increased.  At  the  highest  temperatures 
the  velocity  distributions  determined  from  TOF  were  comparable  in  magnitude 
to  those  observed  at  the  maximum  excitation  conditions  with  the  laser.  However 
in  contrast  to  the  sixfold  enhancement  of  bolometer  signal  detected  with  laser 
excitation,  thermal  heating  produces  a  bolometer  signal  increase  ^W55Q^/,W7qg^) 


of  approximately  a  factor  of  two. 


The  distribution  of  energy  in  the  beam  molecules  between  translational 

and  internal  degrees  of  freedom  following  thermal  heating  of  the  nozzle  is 

calculated  as  above.  The  results  at  room  temperature  and  at  the  maximum 

nozzle  temperature  of  550K  are  shown  in  Table  I.  Also  displayed  are  the 

energy  distributions  obtained  under  optimum  laser  excitation  conditions.  Note 

that  the  contribution  to  kinetic  energy  (K)  are  similar  for  laser  excitation 

and  thermal  heating  at  550K.  However,  the  internal  energies  are  dramatically 

different.'  Laser  excitation  yields  a  SF,  vibrational  temperature  (Ttr__)  of 

O  VId 

^  740K  whereas  thermal  heating  produces  a  nominal  vibrational  temperature  of 
310K.  The  vibrational  temperatures  achieved  by  thermal  heating  are  much  lower 
than  the  nozzle  temperature  indicating  substantial  relaxation  during  the  super 
sonic  expansion  process.  Saturation  of  the  bolometer  during  the  laser  experi¬ 
ments  may  also  have  occurred.  If  so,  this  would  mean  that  the  laser  pumped 
molecules  have  an  even  higher  internal  temperature.  This  point  is  presently 
under  investigation. 


Table  I 

Energy  Distributions  in  the  SF,  Molecular  Beam 

D 


THERMAL 

290K 

THERMAL 

550K 

LASER  EXCITATION 
1  MM  UPSTREAM 

KINETIC  (K) 

920 

1770 

2000 

INTERNAL  (I)  (TylB) 

410  ( 260K) 

680  (310K) 

4100  (740K) 

TOTAL  (T) 

13  30 

2450 

6100 

All  energies  in  cm 
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The  laser  results  may  be  interpreted  via  a  collision-assisted  absorp¬ 
tion  and  up-pumping  mechanism.  During  expansion  the  pressure  and  thus 
collision  rate  decreases  along  the  capillary.  As  the  point  of  irradiation 
is  moved  "up  stream,"  the  collision  frequency  in  a  given  irradiation  volume 
increases,  resulting  in  increased  absorption. The  high  internal  excita¬ 
tion  achieved  by  irradiating  inside  the  nozzle  also  implies  that  many 

vibrational  states  are  collisionally  coupled  to  the  v  pumped  state  of  SF,  . 

d  b 

However,  when  the  collision  frequency  is  further  increased,  V-T  relaxation 

and  energy  transfer  to  the  walls  limits  the  energy  retained  in  SF,  vibration. 

b 

The  increase  in  kinetic  energy  with  laser  excitation  indicates  that  at 

least  some  V-T  relaxation  has  occurred  before  SF,  molecules  reach  the 

b 

collision-free  region  of  the  molecular  beam.  This  observation  is  consistent 

with  expectations  based  upon  the  time  scale  for  relaxation  in  the  bulk  gas 

phase^  at  a  pressure  in  the  range  of  50-75  Torr,  comparable  to  that 

estimated  for  the  SF,  in  the  capillary  immediately  upstream  of  the  nozzle 
o 

exit.  Irradiating  far  upstream  (in  the  high  pressure  region)  should  simulate 
bulk  heating,  yielding  equilibration  among  all  degrees  of  freedom. 

As  demonstrated  in  Table  I,  laser  excitation  inside  the  nozzle  produces 
a  level  of  excitation  and  a  partitioning  of  the  absorbed  energy  which  differs 
significantly  from  thermal  heating  of  the  nozzle.  Thus,  laser  excitation 
even  with  some  collisions  can  provide  a  distinctly  non-equilibrium  distribu¬ 
tion  of  energy  between  internal  and  translational  energies.  Further,  the 
present  results  suggest  that  thermal  heating  even  to  temperatures  at  which 

SF,  would  thermallv  decompose,  would  not  produce  vibrational  excitation 
b 

at  the  level  achieved  with  laser  irradiation  inside  the  nozzle. 

The  present  excitation  technique  provides  a  highly  energetic  molecular 
beam  whose  energy  content  can  be  rapidly  modulated,  synchronous! y  with  the 
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laser.  Problems  such  as  surface  reactions  and  thermal  decomposition 
associated  with  conventional  heated-oven  sources  are  avoided.  This  excita¬ 
tion  method  should  be  of  general  applicability  to  many  molecules  and  is  of 
potential  use  in  reactive  scattering  studies  of  vibrationally  excited 
species^^  and  multiphoton  excitation  diagnostics.^^ 


*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  NSF-CHE77-11384  and  by  the  Department  of  Energy  under  Contract 
DE-AS02-ER-78-S-02-4940 .  Partial  equipment  support  was  provided  by  the 
National  Science  Foundation  under  Grant  NSF-CHE77-24343 . 
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E.  EXCIMER  LASER- INDUCED  PHOTOFRAGMENTATION  AND  VIBRATIONAL  ENERGY 

DISPOSAL  IN  MOLECULES* 

(J.  0.  Chu,  G.  W.  Flynn) 

(Principal  Investigator:  G.  W.  Flynn  (212)  280-4162) 

In  contrast  to  the  numerous  studies  of  collisional  vibrational  relaxa¬ 
tion  of  neutral  molecules,*^  there  have  been  few  investigations  of  the 
reverse  process  of  collisional  vibrational  excitation.  This  can  in  part 
be  attributed  to  the  expected  small  excitation  cross  sections  for  such 
energy  exchange  events.  In  addition,  bulb  techniques  must  be  abandoned  in 
favor  of  molecular  beams  in  order  to  obtain  collisional  energies  above  the 
endothermic  threshold.  On  the  other  hand,  beam  experiments  involving 
neutral  collision  partners  suffer  from  difficulties  associated  with  resolving 
individual  vibrational-rotational  transitions.  Only  one  report  of  a 

molecular  beam  experiment  has  been  made  in  which  a  single  vibrationally 

(2) 

excited  state  under  single-collision  conditions  has  been  resolved. 

(3) 

A  new  bulk  phase  technique  has  recently  been  developed  for  studying 
the  process  of  translational  to  vibrational  (T-V)  energy  transfer  and 
translationally  enhanced  reactive  events.  The  method  involves  excimer 
laser  photolysis  of  bound  molecules  to  produce  translationally  "hot"  (or 
superthermal)  atoms  (or  radical  fragments)  in  a  gas  sample,  followed  by 
direct  time-resolved  detection  of  infrared  fluorescence  from  vibrationally 
excited  product  species.  This  technique  essentially  provides  a  sensitive 
and  direct  probe  of  the  translational-to-vibrational ,  rotational  (T-V,R) 
energy  transfer  dynamics  that  occur  during  each  highly  energetic  collisional 
event.  Furthermore,  infrared  luminescence  detection  of  T-V  excited  products 
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can  readily  be  made  wavelength  state  specific,  and  subsequent  vibrational 
(possibly  rotational)  distributions  can  be  obtained  from  the  direct  wavelength- 
resolved  emission  profile.  In  the  case  of  generating  "hot"  highly  reactive 
radical  fragments,  dynamical  competition  between  the  direct  T-V  energy 
transfer  and  the  reactive  exchange  reactions  can  be  studied  in  detail. 

We  have  initiated  a  series  of  experiments  to  study  the  efficiency  of 
vibrational  excitation  by  energetic  heavy  atom  fragments.  In  particular 
we  have  begun  a  study  of  the  C1^/C0^  and  ICl/CC^  systems.  Schematically 
the  experiments  can  be  represented  by 

Cl^  +  (193  nm)hv  -*•  2C1  "HOT"  Excimer  Laser  Photolysis 

Cl  "HOT"  +  C02  -  C02+(v3)  +  Cl  (T-V)  Energy  Exchange 

C02+  -*■  C02  +  hv3(4.3p)  Infrared  Emission 

In  the  IC1/C02  system,  laser  photodissociation  probably  produces  excited 
2 

"hot"  Cl  (  P^/2^  atoms  which  also  undergo  rapid  reaction  with  IC1  to  produce 
(4) 

Cl2.  We  have  clearly  observed  vibrational  excitation  of  C02  in  these 

experiments  (see  Fig.  1)  and  are  in  the  process  of  analyzing  data.  In  the 
IC1/S02  system,  we  have  also  observed  a  definite  (T-V)  energy  transfer  con¬ 
tribution  in  addition  to  the  SO(v  «  3,4)  near  resonant  transfer  to  SO^v^), 
see  Fig.  2. 

In  addition  to  producing  translationallv  energetic  fragments,  laser 
photodissociation  can  also  produce  fragments  with  excess  vibrational  energies. 
By  studying  the  energy  disposition  produced  in  the  vibrationallv  excited 
fragments,  we  can  probe  and  observe  the  photodissociation  process.  In  our 
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FIGURE  1,  C02(v3)  FLUORESCENCE  CURVES 
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FIGURE  2.  S02(v3)  FLUORESCENCE  CURVE 
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experiments,  we  have  obtained  preliminary  results  on  the  vibrational  energy 
disposal  in  the  photofragmentation  of  (HCN)^-  In  particular,  infrared 
emission  has  been  observed  from  different  vibrational  modes  of  HCN  produced 
via  photofragmentation,  see  Fig.  3. 

A  schematic  of  the  experimental  apparatus  is  shown  in  Fig.  4.  A  pulsed 
Lumonics  TE-860S  excimer  laser  operated  either  at  193  nm  (ArF)  or  248  nm 

2 

(KrF)  was  used.  All  experiments  were  performed  at  low  fluences  (<  3  mJ/cm  ) 
in  order  to  avoid  multiple  photon  events.  Infrared  emission  from  the  system 
is  observed  with  either  a  InSb,  AuGe  or  a  HgCdTe  infrared  detector.  These 
detectors  only  viewed  the  emission  through  the  appropriate  infrared  filters. 
The  output  from  the  detector  is  amplified  and  then  digitized  with  a  transient 
Bromiation  model  8100  recorder.  Digitized  signals  are  summed  in  a  signal 
averager  (Nicolet  1170/Northern  575A)  and  then  stored  in  a  PDP  11-34  computer. 
A  conventional  glass  vacuum  flow  system  is  used  to  provide  a  fresh  gas 
sample  before  each  laser  pulse.  All  gas  samples  are  used  after  multiple 
freeze-pump-thaw  cycles. 

We  have  obtained  preliminary  results  on  the  vibrational  energy  disposal 
in  the  photofragmentation  of  (HCN)^  via  infrared  emission  detection  of  the 
different  vibrational  modes  of  HCN.  Results  have  also  been  obtained  in  the 
study  of  the  relaxation  process  for  electronically  excited  SC^  molecules 
produced  via  excimer  laser  pumping  (see  Fig.  5).  In  addition  laser  photo¬ 
dissociation  has  been  used  to  obtain  heavy  Cl  atoms  with  substantial  excess 
translational  energies.  These  energetic  Cl  atoms  are  found  to  produce 
vibrational  excitation  of  stable  C0^  molecules  during  collision  encounters. 

The  relative  efficiency  of  heavy  vs  light  "hot”  atoms  in  producing  vibra- 
tionally  excited  molecules  is  being  investigated. 
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FIGURE  3.  HCN  FLUORESCENCE  CURVES 


Figure  4:  Experimental  Apparatus  Diagram 
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FIGURE  5.  (ELECTRONIC  RELAXATION)  S02(v3)  FLUORESCENCE  CURVE 
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F. 


PHOTOPHYSICS  AND  PHOTOCHEMISTRY  IN  LOW  TEMPERATURE  SOLID  MATRICES* 


(P.  Beeken,  M.  Mandich,  G.  W.  Flynn) 

(Principal  Investigator:  G.  W.  Flynn  (212)  280-4162) 

Investigation  of  trapped  molecular  species  is  a  growing  area  of 
interest  for  both  molecular  physicists  and  spectroscopists  alike.  In  the 
rare  gas  matrix,  molecules  are  isolated  in  small  cages  where,  unable  to 
translate  and  collisionally  interact,  they  may  be  studied  under  a  variety 
of  excitation  conditions.  A  knowledge  of  the  species  of  interest  as  well 
as  an  understanding  of  the  interactions  with  the  host  lattice  can  be 
obtained. 

Experiments  involving  the  transfer  of  electronic  to  vibrational  energy 
(E  -►  V)  have  been  an  active  area  of  interest  for  several  years  in  the  gas 
phase.  ^  ^  E-*V  studies  can  lead  to  a  better  understanding  of  photo¬ 
chemical  mechanisms,  to  the  development  of  new  pumping  systems  for  molecular 

lasers,  and  can  provide  tests  of  theoretical  data  for  nonadiabatic  energy 

(2) 

transfer  theories.  Because  of  the  many  careful  investigations  in  the 

gas  phase,  the  system  would  be  an  excellent  place  to  begin  an 

inquiry  into  E-*-V  transfer  in  the  matrix. 

There  are  many  questions  to  be  asked  concerning  E-+V  energy  transfer 
in  the  rare  gas  solid.  One  very  important  consideration  is  the  fate  of 
the  halogen  Br^  after  excitation  with  green  light.  When  Br.,  is  excited 
with  green  light  in  the  gas  phase  it  may  dissociate  into  the  constituent 
atoms.  In  the  solid,  however,  the  molecule  is  held  in  a  rigid  cage  and 
cannot  undergo  the  same  process.  We  need,  therefore,  to  investigate  the 
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various  relaxation  pathways  available  to  bromine  in  the  matrix  and  proceed 
to  look  for  possible  transfer  to  other  molecules  such  as  CO^. 

Mixtures  of  freeze-pump-thawed  Br ^  (Alfa  Ventron  99.998%)  and  the 
inert  gas  are  mixed  in  a  glass  and  stainless  steel  vacuum  line  and  sprayed 
on  a  cold  window  using  a  pulsed  deposition  technique.  A  typical  backing 
pressure  of  400  torr  is  used  within  a  10  ml  expansion  volume  to  deposit, 
with  10  to  12  valve  cycles,  roughly  50  ml-torr  of  gas.  The  matrix  is 
maintained  on  a  Csl  window  held  at  temperatures  varying  upwards  from  12K 
by  an  Air  Products  Displex  202A  closed  cycle  Helium  refrigerator  fitted 
with  a  cowling  to  permit  the  mounting  of  windows  transparent  to  excitation 
and  fluorescence  frequencies  (see  Figure  1) . 

Two  sources  of  excitation  are  used  in  these  studies:  a  Chromatix 
CMX4  linear  flash  pumped  dye  laser,  and  a  Molectron  DL200  pumped  dye 
laser.  The  principal  dyes  used  in  this  series  of  investigations  were 
obtained  from  Exciton  Corporation  and  are  LD490,  Coumarin  500,  Rhodamine  6G, 
and  Fluorescin.  The  laser  intensity  was  measured  with  a  Molectron  J3-02DW 
joulemeter  while  fluorescence  is  gathered  by  a  variety  of  detectors. 

Visible  fluorescence  is  obtained  by  an  RCA  C31034A  photomultiplier  and  far 
red  to  near  IR  light  is  seen  by  an  Hamanatsu  R316  S-l  phototube.  Spectra 
were  resolved  using  a  boxcar  integrator  and  an  Instruments  SA  monochromator 
with  a  grating  ruled  with  600  groves/mm  and  blazed  at  700  nm.  IR  fluorescence 
beyond  the  range  of  the  phototube  is  measured  with  an  InSb  detector  using 
various  interference  filters  to  resolve  the  wavelength.  Time  dependent 
signals  are  taken  using  a  transient  recorder  signal  averager  combination. 

Argon  Matrix:  The  emission  spectrum  of  Br^  (1:500  mixture)  from  a  500  nm 
excitation  appears  in  Fig.  2  and  Fig.  3.  It  should  be  noted  that  the 
asymmetry  of  the  isotopic  triplets  is  indeed  a  real  phenomenon  and  is  not  an 
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Figure  1:  Schematic  of  Experimental  Apparatus. 

V  *  Manual  Value;  SV  ■  Solenoid  Value;  Mono  = 
Monochromator;  PM  «  Photomultiplier;  HV  *■  High 
Voltage  Power  Supply;  InSb  -  Indium  Antimonide 
IR  Detector;  F  =  Filter;  Trans.  Dig.  *  Transient 
Digitizer /Signal  Averager. 


Ill 


-i  rv  ♦.*  I'.  a'  1  z  I . 


3r  :  Ar  1:500 
Slow  Decay  10.  2.  msec 

Fast  Decay  10- .  +  ■-  9.  msec' 
Taken  with  ar.  F-l  ?>'  at  9v>H 


Br,:  At  1:500 
Slow  Decav  .19  + ; 
Fast  Decav  15.9 
Taken  with  an  InSr 
above  1.0  microns 


artifact  of  the  recording  process.  This  is  a  point  which  will  be  alluded 
to  later  in  the  discussion.  The  time  dependence  of  the  various  peaks, 
including  different  regions  of  the  triplets,  have  been  measured  and  appear 
in  Table  I.  The  lifetimes  of  the  various  emissions  did  not  vary,  within 
experimental  error,  as  both  temperature  and  concentration  were  changed. 

A  sample  fluorescence  appears  in  Fig.  3.  IR  emission  is  also  observed 
above  1.0  microns  but  was  not  quite  intense  enough  to  resolve  with  the 
monochromator.  The  IR  emission  is  not  seen  when  a  red  excitation  source 
was  used. 

Krypton  Matrix:  Although  no  spectrum  has  been  obtained  from  the  Br^/Kr 
mixture  (1:1000)  the  time  dependence  of  the  signals  reveal  no  anomalous 
behavior  with  respect  to  the  Ar  data  aside  from  a  slight  shift  in  lifetime. 

As  it  was  with  the  Ar  sample  the  time  characteristics  show  no  appreciable 
sensitivity  to  changes  in  temperature.  Work  is  currently  in  progress  to 
obtain  the  emission  profile  and  gather  frequency  dependent  lifetimes. 

Xenon  Matrix:  It  is  in  this  matrix  that  the  Br^  sample  (1:1000)  sho..s 
a  significant  deviation  from  the  Ar  data.  The  spectrum's  most  significant 
feature  (Fig.  4)  is  the  complete  loss  of  the  triplet  structure  so  pronounced 
in  the  Ar  spectrum.  In  addition,  the  time  dependent  information  illustrated  in 
Table  I  show  a  considerable  deviation  from  the  expected  trend  in  going  from 
Ar  to  Kr.  The  one  aspect  common  to  the  previous  studies  is  that  the  life¬ 
times  are  indeed  independent  of  temperature. 

Other  work  has  been  reported  on  spectroscopy  of  bromine  in  rare  gas 
matrices.  a  red  (614  nm)  excitation  source  was  used  to  obtain  a 

lifetime  of  7.3  usee  as  well  as  an  emission  and  excitation  profile  of  Br,.^^ 
In  a  series  of  studies  on  trapped  halogen  molecules  the  emission  spectra 
were  assigned  to  transitions  from  the  ground  vibrational  level  of  the 
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excited  electronic  state  to  the  various  vibrational  levels  of  the  ground 
electronic  well.  The  work  concluded  by  assigning  the  emission  solelv  to  the 
B3it(0+u)  state  (see  Fig.  5).  In  a  similar  study  isolated  samples  of  Cl2  were 
irradiated  with  a  source  that  would  correspond  to  a  gas  phase  dissociation 
limit .  Here  the  investigators  reported  that  the  \(2u)  was  the  only  source 
of  the  emission  in  the  trapped  chlorine.  We  have  investigated  the  trapped 
halogen,  Br,,,  under  conditions  which  closely  resemble  the  Cl2  work;  we  irradiate 
with  a  source  which  corresponds  to  the  gas  phase  dissociation  limit  of  Br . 

When  we  examined  the  time  resolved  data,  we  found  toward  the  red  of  the 
isotopic  triplets  two  exponential  decays.  The  blue  side  of  the  triplet 
yields  only  one  decay.  Closer  examination  of  the  spectra  does  suggest 
that  there  is  a  new  resonance  growing  into  the  emission  profile  toward 
the  red  peak  of  the  triplet.  Toward  the  IR  this  peak  resolves  itself 
separately  from  the  isotopic  lines. 

We  examined  lifetimes  of  bromine  further  to  the  infrared,  although 
the  sensitivity  of  our  detection  system  would  not  allow  us  to  use  a  mono¬ 
chromator.  We  did  discover  a  decay  rate  very  close  to  the  slow  decay 
measured  in  the  visible/near  IR  region  of  the  spectrum.  The  Franck-Condon 
profiles^3  as  given  in  Fig.  6  suggest  that  emission  from  the  IR  would 
reflect  solely  the  lifetime  of  the  A  state  while  emission  in  the  visible 
could  most  likely  be  due  to  B  state.  The  region  sampled  in  between, 
therefore,  should  exhibit  characteristics  of  both  emissions.  Our  emission 
profile  coupled  with  the  time  dependent  data  obtained  from  the  trapped 
bromine  strongly  suggests  that,  indeed,  two  states  are  responsible  for  the 
emission. 

When  bromine  is  isolated  in  xenon  quite  an  unusual  profile  is  obtained 
and  very  short  lifetimes  are  observed  (see  Table  I  and  Figure  4  for  the 
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fbfENTIAL 


Bromine  Potential  Diagram  with  gas  phase  absorption 
coef  f ic ient . 


spectrum) .  The  lifetime  in  the  IR  region,  as  well  as  elsewhere,  yielded 

an  unexpected  deviation  from  the  trend  seen  in  going  from  Ar  to  Kr  to  Xe. 

The  most  obvious  feature  of  the  spectra  is  the  complete  loss  of  the  structure 

so  readily  obtained  in  the  Ar  sample  and  the  loss  of  a  second  decay.  We 

believe  that  the  lifetime  seen  in  the  IR  is  the  lifetime  observed  further 

into  the  visible  region  of  the  spectrum. 

The  primary  problem  at  this  point  is  to  try  to  determine  what  pathways 

are  available  to  Br^  if  it  cannot  dissociate  as  it  would  in  the  gas  phase. 

The  Ar  data  suggests  that  the  sample  may  indeed  dissociate  within  the  cage 

(4) 

and  reform  in  an  excited  state.  Previous  work  has  shown  through  selective 

isotopic  excitation  that  very  little  electronic  to  electronic  energy  transfer 

occurs  even  between  the  nearly  resonant  isotopic  neighbors.  This  implies 

that  if  the  bromine  molecule  were  excited  into  one  electronic  state,  it 

would  not  be  likely  to  transfer  the  excitation  to  another  state  of  another  Br^. 

The  dissociation/recombination  idea  is  strongly  supported  by  a  lack  of  a 

rise  time  in  any  of  the  molecular  fluorescence  signals.  In  Xe  the  excited 

Br^  may  be  forming  an  excimer  of  the  form  Br^-Xe  which  would  broaden  the 

emission  profile  while  leaving  the  resonance  spacings  close  to  the  value 

expected  for  Br^.  The  close  interaction  of  the  bromine  and  xenon  might 

then  account  for  the  extremely  short  lifetime. 

Bromine  under  green  excitation  yields  an  emission  profile  from 

different  emitting  states.  The  spectra  coupled  with  lifetime  studies 

3  + 

indicates  that  the  primary  source  of  the  emission  is  from  the  B  tt (0  u)  and 
3 

A  r(lu)  states.  The  data  strongly  suggest  that  the  mechanism  for  production 
of  population  in  all  these  states  may  be  a  dissociation/recombination 
reaction,  where  the  new  molecule  can  reform  in  a  variety  of  excited  states. 

A  different  mechanism  of  interaction  is  indicated  by  the  results  from  the 
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Xe  spectra;  the  excited  molecule  may  be  forming  an  excimer  with  its 
host.  With  an  increased  understanding  of  Br^  in  these  assorted  rare  gas 
hosts  we  can  begin  to  evaluate  possibilities  of  transfer  of  energy  to  other 
species  within  the  same  matrix  material. 

Our  initial  impetus  for  pursuing  the  investigation  into  the  nature  of 
the  excited  bromine  was  to  study  a  potential  mechanism  for  electronic  to 
vibrational  energy  transfer  in  a  rare  gas  matrix.  The  Br^-CO^  system  has 
been  well  characterized  in  the  gas  phase  and  our  ultimate  goal  is  to 
investigate  the  possibility  for  transfer  of  the  electronic  excitation  to 
vibrational  degrees  of  freedom. 

*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  NSF-CHE77-11384  and  by  the  Department  of  Energy  under  Contract 
DE-ASO2-ER-78-S-02-4940. 

(1)  A.  Hairi  and  C.  Wittig,  J.  Chem.  Phys.  617,  4454  (1977). 
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IV.  ENERGY  TRANSFER  AND  RELAXATION  IN  ALKALI  METALS 

A.  NUCLEAR  SPIN  POLARIZATION  OF  Xe  by  OPTICAL  PUMPING* 

(B.  Suleman,  T.  McClelland,  N.  D.  Bhaskar,  W.  Happer) 

(Principal  Investigator:  W.  Happer  (609)  452-5584) 

129 

We  have  successfully  polarized  Xe  by  spin  exchange  with  optically 
pumped  rubidium.  The  polarization  was  detected  by  adiabatic  reversal  of 
the  xenon  spins,  as  described  in  the  preceding  progress  report.  Work 
on  this  project  is  continuing  at  Princeton  University. 

*This  research  was  also  supported  by  the  Air  Force  Office  of  Scientific 
Research  under  Grant  AFOSR-7 9-0082 . 


B.  SPECTROSCOPY  OF  THE  5D  STATE  OF  CESIUM* 


(A.  Sharma,  N.  D.  Bhaskar,  W.  Happer) 

(Principal  Investigator:  W.  Happer  (609)  452-5584) 

We  have  discovered  that  the  cesium  atoms  within  our  glass  sample  cells 

lase  continuously  on  many  transitions  in  the  infrared  spectral  region, 

e.g.  7P  to  6S  at  3 . 095u  and  2.93p;  7S  to  6P  at  1.469p  and  1.359p;  7P  to 

5D  at  1.376p  and  1.360u.  The  atoms  are  excited  with  cw  laser  radiation 
o  o 

at  4593A  and  4555A  and  the  lasing  occurs  without  any  feedback  mirrors. 

We  believe  this  is  the  first  reported  example  of  mirrorless  cw  lasing. 

Work  on  this  project  is  continuing  at  Princeton  University. 

*This  research  was  also  supported  by  the  Army  Research  Office  (Durham) 
under  Contract  DAAG29-80-C-0043 . 
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c.  INFRARED  AESORPTION  OF  ALKALI  MOLECULES* 

(A.  Vasilakis,  N.  D.  Bhaskar,  W.  Happer) 

(Principal  Investigator:  W.  Happer  (609)  452-5584) 

We  have  completed  our  studies  of  the  newly  discovered  infrared 
absorption  bands  of  dense  alkali  vapors  of  Cs,  Rb,  K  and  Na.  All  elements 
show  very  similar  structure  between  lp  and  2p,  structure  which  may  be 
due  to  absorption  from  triplet  molecules.  Preliminary  attempts  to  observe 
these  bands  in  fluorescence  failed  to  show  any  signal.  We  will  soon 
repeat  these  experiments  with  a  much  better  infrared  detector,  intrinsic 
germanium  with  a  sensitivity  at  least  200  times  better  than  the  PbS  detector 
we  used  in  our  first  experiments.  We  have  also  constructed  a  high  pressure 
alkali-vapor  jet  which  will  allow  us  to  attain  much  better  solid  angles  of 
viewing  than  in  our  earlier  experiments.  This  work  is  continuing  at 
Princeton  University. 


*This  research  was  also  supported  by  the  Army  Research  Office  (Durham) 
under  Contract  DAAG29-80-C-0043. 


D.  GROUND  STATE  POLARIZATION  BY  TWO-PHOTON  OPTICAL  PUMPING  OF  ATOMIC 
VAPORS* 

(N.  Tran,  N.  D.  Bhaskar,  W.  Happer) 

(Principal  Investigator:  W.  Happer  (609)  452-5584) 

We  have  detected  strong  two  photon  excitation  of  rubidium  according 
to  the  scheme  discussed  in  the  previous  progress  report.  Substantial 
(i.e.  10%  to  20%)  spin  polarization  is  also  detected.  However,  our 
analysis  of  the  data  is  complicated  by  the  occurrence  of  strong  excimer 
absorption  which  also  leads  to  spin  polarization  of  a  comparable  order 
of  magnitude.  This  work  is  continuing  at  Princeton  University. 

*This  research  was  also  supported  by  the  Naval  Air  Development  Center, 
Warminster,  PA. 
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E.  TEMPERATURE  AND  DENSITY  VARIATION  OF  THE  INFRARED  BANDS  OF  CsXe ; 

CESIUM  POLYXENIDE  BANDS* 

(R.  Novak,  N.  D.  Bhaskar,  W.  Hamper) 

(Principal  Investigator:  W.  'Happer  (609)  452-5584) 

This  work  has  been  completed  and  the  results  have  been  published  . 

*This  research  was  also  supported  by  the  Army  Research  Office  (Durham)  under 
Contract  DAAG29-80-C-0043 . 

(1)  R.  Novak,  N.  D.  Bhaskar,  and  W.  Happer,  J.  Chem.  Phys.  _71,  4052 
(1979). 
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FARADAY  ROTATION  AND  MAGNETIC  CIRCULAR  DICHROISM  OF  ALKALI-NOBLE 
GAS  SYSTEMS* 


(M.  Islam,  W.  Happer) 

(Principal  Investigator:  W.  Happer  (609)  452-5584) 

Extensive  measurements  of  the  magnetic  circular  dichroism  cf  the 
CsAr  excimer  bands  and  Cs^  dimer  bands  have  been  made  during  the  past 
year.  We  have  been  able  to  gain  for  the  first  time  unambiguous  information 
about  the  angular  momentum  and  angular  momentum  coupling  schemes  of  the 
electronic  states  of  excimer  molecules. 

*This  research  was  also  supported  by  the  Army  Research  Office  (Durham) 
under  Contract  DAAG29-80-C-0043 . 


G.  PARTICULATE  FORMATION  IN  CESIUM-DEUTERIUM  CAS  CELLS* 

(J.  Camparo,  J.  Pietras,  N.  D.  Bhaskar,  W.  Happer) 

(Principal  Investigator:  W.  Happer  (609)  452-5584) 

We  have  used  spin  exchange  magnetic  resonance  spectroscopy  to 
unambiguously  identify  free  hydrogen  atoms,  free  deuterium  atoms  and 
free  electrons  in  cells  containing  mixtures  of  hydrogen  or  deuterium 
and  cesium  vapor  which  is  optically  pumped  by  the  4593^  second  resonance 
line  of  the  cesium  atom.  The  results  will  be  published  in  a  May  1981 
issue  of  Physical  Review  Letters.  This  work  is  continuing  at  Princeton 
University. 

*This  research  was  also  supported  by  the  Air  Force  Office  of  Scientific 
Research  under  Grant  AFOSR-7 9-0082 . 
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H.  SPIN  DESTRUCTION  IN  COLLISIONS  BETWEEN  CESIUM  ATOMS* 

(J.  Camparo,  J.  Piecras,  N.  Bhaskar,  W.  Happer) 

(Principal  Investigator:  W.  Happer  (609)  452-5584) 

This  work  has  been  completed  and  the  results  were  published. ^ 
Closely  related  work  on  "light  narrowing,"  a  narrowing  of  the  magnetic 
resonance  lines  of  optically  pumped  alkali  vapors  by  intense  pumping 
light  has  also  been  completed  and  a  paper  summarizing  the  results  will 
be  published  in  Physical  Review  in  May  or  June  of  1981. 

*This  research  was  also  supported  by  the  Air  Force  Office  of  Scientific 
Research  under  Grant  AF0SR-79-0082 . 

(1)  N.  D.  Bhaskar,  J.  Pietras,  J.  Camparo,  and  W.  Happer,  Phys.  Rev. 
Letters  44 ,  930  (1980) . 
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I. 


SPIN  EXCHANGE  AND  RELAXATION  IN  Na-NOBLE  GAS  MIXTURES* 


(B.  Suleman,  N.  D.  Bhaskar,  T.  McClelland,  W.  Happen) 

(Principal  Investigator:  W.  Happer  (609)  452-5584) 

The  spin  destruction  rates  of  sodium  in  xenon  have  been  measured 
and  found  to  be  even  greater  than  the  rates  of  rubidium  in  xenon. 

Van  der  Waals  molecules,  NaXe,  are  also  found  to  play  an  important 
role  in  the  relaxation.  Interium  results  have  been  published^  and 
work  on  this  project  is  continuing  at  Princeton  University. 

*This  research  was  also  supported  by  the  Air  Force  Office  of  Scientific 
Research  under  Grant  AFOSR-79-0082. 


V.  PICOSECOND  ENERGY  TRANSFER  AND  PHOTOFRAGMENTATION  SPECTROSCOPY 


A.  PICOSECOND  LASER  STUDIES  OF  PHOTODISSOCIATION  AND  ENERGY  RELAXATION 

IN  SHORT-LIVED  FRAGMENTS* 

(C.  Dupuy,  G.  M.  Korenowski,  M.  McAuliffe,  W.  M.  Hetherington  III, 

K.  B.  Eisenthal) 

(Principal  Investigation:  K,  B.  Eisenthal  (212)  280-3175) 

Molecules  of  the  diazo  family  C  =  N^  are  of  widespread  utility  in 
microelectronics  and  semiconductor  fabrication,  and  in  the  general  use  of 
photosensitive  materials  for  information  recording  and  lithography . ^ ^ 
For  many  molecules  of  this  type  the  dominant  channel  for  electronic  energy 
relaxation  following  ultraviolet  photoexcitation  (or  electron  beam  excita¬ 
tion)  is  bond  rupture  generating  nitrogen  and  an  important  short  lived  and 
reactive  intermediate  called  a  carbene.^^  We  can  write  the  photo¬ 
dissociation  as 


Diazo  Compound  Carbene 


where  R^  and  R£  represent  the  groups  attached.  There  has  been  a  long  standing 

and  lively  interest  in  the  structural  and  chemical  properties  of  carbenes  due 

(3)-(I8) 

to  their  novel  reactions  and  their  distinctive  electronic  structures.  ' 

The  lowest  lying  singlet  and  triplet  states,  either  of  which  can  be  the  ground 

(3)  (4) 


state  as  determined  by  the  attached  groups  R^  and  R^,  are  close  in  energy. 


This  proximity  of  singlet  and  triplet,  each  having  a  different  chemistry, 
makes  the  kinetics of  interconversion  a  key  aspect  of  the  properties  of 


(3) (4) 


these  important  chemical  intermediates . 


In  order  to  elucidate  the  dynamics  of  the  photodissociation  of  the  parent 
diazo  compound  and  the  subsequent  energy  relaxation  and  reaction  of  the  carbene 
generated,  an  ultraviolet  picosecond  laser  pulse  is  used  for  photoexcitation. 

The  parent  molecule  was  diphenyldiazomethane.  It  was  dissolved  in  3-methyl 
pentane  in  acetonitrile  at  25°C  and  irradiated  at  266  nm.  The  sequence  of 
energy  decay  steps  following  photoexcitation  of  diphenyldiazomethane,  ultimately 
leading  to  diphenyl  carbene  in  its  triplet  ground  state  can  be  "written  as: 


Excited  Singlet 
Diphenyl diazome thane 


Ground  State 
Triplet 

Diphenyl carbene 


where  Ph  is  the  phenyl  group 


<2>~ 


In  our  initial  work 


(19) 


we  estab- 


3  * 

lished  that  a  low  yield  of  the  excited  state  triplet  carbene  DPC  is  produced 

on  one  photon  excitation  of  the  parent  compound.  Knowing  the  energy  of  the 

3  * 

exciting  photon  and  the  energy  of  DPC  we  conclude  that  the  energy  of  the 


carbon-nitrogen  bond  in  diphenyldiazomethane,  heretofore  unknown,  is  less 
than  2eV.(20) 
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In  our  work  reported  here^^  we  have  used  a  streak  camera  to  detect  the 
3  * 

rise  time  of  the  DPC  to  be  <15  ps.  Thus  the  bond  rupture  and  excited  state 
* 

interconversion,  k^,  in  c^e  carbene  must  occur  in  less  than  15  ps.  To  ob- 

3 

tain  the  kinetics  of  formation  of  the  ground  state  triplet  carbene,  DPC,  we 

used  a  laser  induced  fluorescence  method.  Measurement  of  the  triplet  fluores- 

3  * 

cence  intensity  from  DPC  as  a  function  of  the  time  separation  between  the 

excitation  and  probe  pulses  yielded  a  value  of  110  ps  for  the  rise  time  of 

3  3  * 

DPC  (Fig.l).  Since  we  independently  measured  the  lifetime  of  DPC  to  be  4  ns, 

3 

we  conclude  that  the  formation  route  of  DPC,  which  occurs  in  110  ps,  cannot 
3  * 

be  from  DPC  .  The  ground  state  triplet  must  therefore  come  primarily  from 
■*DPC.  The  value  of  k^  is  thus  110  ps. 

From  our  collaborative  work  with  Professor  N.J.  Turro  we  know  the  value 
of  the  reverse  process,  namely  k  ,  is  4.1  x  10^  sec  3 .  We  therefore  find 

lb 

that  the  equilibrium  constant  for  the  process 


ST 


DPC 


*  3 


DPC 


TS 


is 


K  =  -9--1  x-V--gec-Tl  =  2.2  x  103 

4.1  x  10  sec 


We  calculate  the  free  energy  difference  at  25°C  to  be  4.6  kcal/mole.  After 
correcting  for  the  entropy  difference  arising  from  the  multiplicity  differences 
of  the  singlet  and  triplet  states  we  find  an  energy  separation  between  the 
ground  triplet  and  first  singlet  state  to  be  3.9  kcal/mole. 


*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  NSF-CHE79-23291  and  by  the  Air  Force  Office  of  Scientific  Research 
under  Grant  AF0SR-77-3407B  (expired  9/30/80)  and  AF0SR-81-0009 . 
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Figure  1:  Formation  kinetics  of  DPC.  The  laser  induced 

fluorescence  signal  as  a  function  of  the  time  delay 
between  excitation  and  probe  pulses.  The  solid 
curve  is  a  calculated  fit. 
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B.  PICOSECOND  LASER  STUDIES  OF  THE  DYNAMICS  OF  CHAIN  RELAXATION* 

(M.  K.  Crawford,  Y.  Wang,  K.  B.  Eisenthal) 

(Principal  Investigator:  K.  B.  Eisenthal  (212)  280-3175) 

The  dynamics  of  chain  molecules  in  solution  is  currently  a  subject  of 
great  interest  and  activity. ^  ^  In  the  present  study  we  have  used  intra¬ 
molecular  electron  transfer  between  the  two  end  groups  of  a  chain  molecule 

(9) 

to  probe  chain  dynamics.  We  have  shown  in  earlier  work  that  intramolecular 

exciplex  formation  in  nonpolar  solvents  is  governed  by  the  rate  of  conforma¬ 
tional  changes.  Thus  by  monitoring  the  fluorescence  decay  of  the  initially 
excited  group  and  the  rise  of  the  exciplex  fluorescence  various  chain  relaxa¬ 
tion  processes  can  be  examined. 

In  the  present  study  anthracene  and  N , N-dimethylaniline  are  located 
at  each  end  of  the  chain  molecule  consisting  of  three  CH^  groups.  A- (CH^) ^-DMA. 
It  is  therefore  the  end  to  end  relaxation  mode  that  we  are  studying  in  this 
;-jo'",ecule.  Using  a  picosecond  laser  pulse  to  excite  the  anthracene  part  we 
monitor  the  fluorescence  decay  of  anthracene  and  the  exciplex  fluorescence 
rise  with  a  picosecond  streak  camera  (Fig.  1).  In  nonpolar  solvents  we  have 
found  that  the  end  to  end  relaxation  of  a  four  bond  chain  is  exponential, 
i.e.  the  rate  of  relaxation  is  independent  of  time.  This  is  the  first  experi¬ 
mental  observation  of  the  dynamics  of  end  to  end  motion  in  short  chain  mole¬ 
cules.  There  have  been  a  number  of  theoretical  and  computer  simulation  treat¬ 
ments  of  the  end  to  end  motions  of  polymer  chains,  some  as  short  as  ten 
bonds.  We  are  presently  engaged  in  analyzing  these  theories  in 

terms  of  our  experimental  findings. 

*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  NSF-CHE79-23291  and  by  the  Air  Force  Office  of  Scientific  Research 
under  Grant  AFOSR-77-3407B  (expired  9/31/80)  and  AFOSR-81-0009 . 
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Figure  1:  Decay  of  the  anthracene  moiety  measured  at  410  run 
and  the  formation  of  the  exciplex  measured  at  two 
wavelengths,  480  nm  and  520  run.  The  dots  are 
experimental  points  and  the  curves  are  theoretical 
calculations . 
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c.  ELECTRONIC  ENERGY  RELAXATION  BY  CHARGE  TRANSFER  PROCESSES* 

(M.  K.  Crawford,  Y.  Wang,  K.  B.  Eisenthal) 

(Principal  Investigator:  K.  B.  Eisenthal  (212)  280-3175) 

Charge  transfer  (CT)  interactions  between  an  excited  state  molecule  and 
surrounding  molecules  provide  an  extremely  important  and  common  pathway  for 
electronic  energy  relaxation. In  nonpolar  media  a  new  and  transient 
species  called  an  excited  CT  complex  is  produced,  whereas  in  polar  media  the 
formation  of  ions  can  be  observed.  From  our  previous  studies^^  on  the 

system  anthracene  which  in  its  excited  state  acts  as  an  electron  acceptor  A* 
and  N, N-dimethylaniline  as  the  donor  D,  we  have  found  that  an  electron  jump 
from  D  to  A*  over  a  reaction  distance  of  8  X  in  less  than  10  ^  sec  (10  ps)  . 

A*  +  D  .lQ.li.Isec  >  (A-  -  D+)* 


Excited 

Anthracene  N, N-dimethylaniline  Excited  CT  Complex 

By  connecting  the  anthracene  and  N, N-dimethylaniline  with  three  CH0 
groups  we  have  imposed  geometric  restrictions  on  the  relative  orientational 
motions  of  A  and  D  while  still  maintaining  their  close  proximity  (  <4  X) . 
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In  this  way  we  have  been  able  to  determine  the  role  of  geometry  on  the  electron 
transfer  process .  (H)  ye  jlave  fountj  that  in  nonpolar  media  the  formation 

of  the  excited  CT  complex  for  the  takes  nanoseconds  to  achieve 

versus  picoseconds  (10  ^  sec)  for  the  free  A  and  D.  We  have  therefore  learned 
chat  there  are  specific  geometries  required  to  effect  electron  transfer.  For 
the  molecule  A-CQ^j-D  the  formation  time  is  determined  by  rotations  of  A  and 
D  about  the  CH^  groups  into  the  appropriate  electron  transfer  conformation. 

The  issue  we  then  proceeded  to  examine  was  the  role  of  the  solvent  di¬ 
electric  constant  on  the  electron  transfer  process  and  its  interplay  with  the 
geometric  restrictions  in  the  A— (CHjK-D  molecule.  To  this  end  we  have  studied 
A-CCH^I^-D  in  the  strongly  polar  solvent  acetonitrile  (e  =  37)  using  picosecond 
fluorescence  as  a  probe  method.  Following  picosecond  laser  excitation  of  the 
anthracene  moiety  at  351  ran  (third  harmonic  of  a  Nd:  phosphate  glass  system) 
the  decay  kinetics  of  the  excited  anthracene  singlet  was  monitored  at  411  nm 
and  the  formation  kinetics  of  the  excited  CT  complex  was  monitored  in  the 
500-600  ran  region  (Fig.  1). 

Although  the  CT  complex  is  formed  "instantaneously",  we  find  that  the 
excited  anthracene  moiety,  which  Is  monitored  at  411  nm,  has  a  lifetime  of 
7  +  1  ps  (Fig.  2).  Since  the  CT  rise  is  faster  than  the  anthracene  decay  we 
conclude  that  these  observed  emissions  correspond  to  independent  processes. 

The  very  fast  CT  rise  suggests  that  only  those  molecules  in  a  favorable 
ground  state  conformation  at  the  time  of  excitation  lead  to  the  excited  com¬ 
plexes.  For  those  molecules  in  unfavorable  extended  conformations  a  new  non- 
radiative  decay  channel  is  responsible  for  the  observed  7  ps  lifetime  of  the 
excited  anthracene  moieties.  This  channel,  which  becomes  accessible  in  highly 
polar  acetonitrile,  is  direct  formation  of  ion  pairs.  For  extended  molecules, 
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1.2  PICOSEC/ CHANNEL 


Figure  2:  Anthracene  moiety  fluorescence  in  acetonitrile.  The 

laser  pre-pulse  is  both  a  pulse  width  display  and  a  time 
marker  for  signal  averaging.  Solid  curve  is  a  theoretical 
fit  to  the  experimental  points. 
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the  excited  anthracene  decays  directly  to  the  ion  pair  since  the  internal  ro¬ 
tational  motions  necessary  for  the  formation  of  excited  CT  complexes  are  too 
slow  to  compete.  The 'appearance  of  this  new  channel  for  energy  dissipation 
in  polar  solvents  is  a  result  of  solvent  stabilization  of  the  ion  pair.  On 
the  other  hand,  in  the  nonpolar  solvents,  e.g.  2-methylbutane ,  it  is  probable 
that  the  energy  of  A  is  higher  than  that  of  A*-(CH,,)  -D  in  extended 

conformations.  The  following  scheme  summarizes  the  results  in  very  polar 
med  ia : 

A-D  A* 


where  represents  the  sandwich  conformation,  A-D,  an  extended  conforma- 

—  -f-  —  +  * 

tion,  A  -D  the  solvent  stabilized  extended  ion  pair,  and  (A  D  )  the 
s  s  v — * 

excited  CT  complex. 

We  are  extending  our  studies  to  other  systems  to  test  our  conclusions 
on  the  nature  of  dielectric  and  geometric  effects  on  charge  transfer  phenomena 


in  liquids. 
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VI.  GENERATION  AND  CONTROL  OF  RADIATION 


A.  SPONTANEOUS  AND  INDUCED  COHERENT  RADIATION  GENERATION  AND  CONTROL  IN 
ATOMIC  VAPORS* 

(K.  P.  Leung,  R.  Kachru,  E.  Whittaker,  T.  W.  Mossberg,  S.  R.  Hartmann) 
(Principal  Investigator:  S.  R.  Hartmann  (212)  280-3272) 

Because  of  the  inherent  doppler-free  nature  of  the  photon  echo  we 

were  able  to  investigate  the  sub-doppler  regime  of  the  sodium  D  line  where 

we  saw  for  the  first  time  the  effect  of  velocity  changing  collisions . ^ 

Our  first  experiments  used  He  gas  to  perturb  the  Na  resonance  and  were 

reported  in  last  year's  progress  report.  In  the  absence  of  velocity  changing 

collisions,  the  broadening  cross  section  is  independent  of  the  photon- 

echo  laser-excitation-pulse-separation-time  t.  We  found  instead  that  o  ^ 

2 

had  the  form  o  r,  -  a  +  at  .  It  has  been  shown  that  in  the  collision 
ef  f  o 

(2) 

kernel  approximation  the  general  expression  for  is  given  by 

°eff  *  a0  +  (c»  "  aa)  t1  ~  “  /T  dt  S<kt>l  (1) 

0 


where  a  -  a  at  t  s  0,  a  s  a  at  t  =  k  is  the  wave  vector  associated  with 
o  00 

the  resonance  line  and  g  is  the  Fourier  transform  of  the  collision  kernel 

g(v) .  The  collision  kernel  gives  the  probability  that  a  collision  will 

result  in  a  velocity  change  Aw  where  v  is  the  initial  velocity.  Since  g(v) 

is  supposed  to  be  an  even  function  of  v  it  follows  that  for  small  t  the 

2 

cross  section  a  has  the  form  a  =  a  +  ar  .  The  constant  a  is  equal  to 

o 

12  2  2 

(a  -  a  )—  k  (v^)  where  (  v  )  is  the  second  moment  of  the  collision  kernel, 
oo  o  o 

We  have  extended  our  photon  echo  relaxation  measurements  to  the  other 
noble  gases  where  we  find  similar  evidence  of  velocity  changing  collisions. 
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We  have  used  these  results  to  measure  g  the  phase  changing  collision 

cross  section, and  we  find  reasonable  agreement  with  more  standard  line 

(3) 

broadening  measurements.  See  Table  I  where  we  compare  broadening 

constants  in  terms  of  FWHM  with  the  results  obtained  by  other  workers. 

There  are  serious  problems  associated  with  standard  line  broadening  measure¬ 
ments  and  the  echo  technique  which  allows  work  at  low  perturber  pressures 
is  singularly  free  of  them.  Most  notable  among  the  advantages  of  the  echo 
technique  is  the  absence  of  optical  pumping  effects  and  the  freedom  from 
high  frequency  stability  requirements. 

We  have  just  completed  a  series  of  two-pulse  photon  echo  relaxation 
measurements  on  the  Li  D  line  perturbed  by  the  noble  gases.  Here  the  effect 
of  velocity  changing  collisions  is  more  dramatic  and  the  cross  section  has 
the  interesting  dependence  on  t  shown  in  Fig.  1.  The  solid  lines  are 
theoretical  fits  using  a  gaussian  kernel  for  all  perturbers  except  He  where 
we  use  a  Lorentzian.  An  argument  based  on  the  uncertainty  principle 
insures  that  we  have  come  close  to  the  asymptotic  limit  where  ae£j  becomes 
constant.  These  measurements  are  readily  analyzed  by  writing  (1)  in  the 
asymptotic  form 

o  ,,  *  a  +  (a  -  a  )  [1  -  ( rr/ 2)  g(o)/kx]  (2) 

eft  o  00  o 

so  that  a  plot  of  °effT  vs  T  should,  for  large  enough  t,  yield  a  straight 

line  whose  slope  gives  ct^  and  whose  intercept  gives  (o^  -  o  ) (ir/2)g(o) /k. 

This  plot  is  shown  in  Fig.  2.  If  we  use  the  fit  of  Fig.  1  to  determine 

a  we  obtain  g(o)  from  our  knowledge  of  (o  -  a  ) (n/2)g(o) /k.  We  have 
O  00  o 

obtained  a  simple  interpretation  for  o  and  o  which  we  find  can  be  represent- 

o 

ed  by 
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Table  I 


Comparison  of  the  results  obtained  in  the  present  experiment  with 
previous  work.  The  broadening  constants  represent  the  FWHM  in 
units  of  10~9  rad  S-1  cm^ . 


Broadening  constants 

.This  work 

Previous 

(a) 

work 

(b) 

Theory 

(c) 

He  Dx 

3.57(0.23) 

3.55(0.09) 

3.80(0.2) 

3.76 

d2 

3.90(0.25) 

4.30(0.03) 

4.38(0.4) 

3.94 

Ne  D 

2.27(0.15) 

2.79(0.05) 

2.28(0.2) 

D2 

2.37(0.16) 

2.53(0.03) 

2.88(0.2) 

Ar  Dx 

4.33(0.28) 

5.56(0.04) 

5.54(0.4) 

°2 

4.31(0.27) 

4.22(0.10) 

4.54(0.4) 

Kr  Dx 

4.51(0.29) 

4.83(0.09) 

D2 

4.33(0.29) 

4.73(0.11) 

Xe  D1 

5.00(0.32) 

5.83(0.12) 

5.26(0.4) 

D2 

5.08(0.33) 

4.94(0.11) 

4.64(0.5) 

aR.  H. 

,  Chatham,  A. 

Gallagher, 

and  E.  L. 

Lewis , 

J.  Phys.  B:  Atom.  Molec.  Phys.  22*  L7  (1980) . 

^N.  Lwin,  D.  G.  McCartan,  and  E.  L.  Lewis, 

J.  Phys.  B:  Atom.  Molec.  Phys.  9_,  L161  (1976); 
Astrophys.  J.  213  599  (1977). 

CA.  D.  Wilson  and  Y.  Shimoni,  J.  Phys.  B: 

Atom.  Molec.  Phys.  _8>  2415  (1975). 
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1 


=  Re  /  (jfc  -  l)-ds/jk-k  (3) 

Sl 

where  Che  index  Z  is  either  0  or  •.  In  this  expression  j  represents 
the  amplitude  of  the  dipole  current  density  t|>*pVi p  -  Vtp*ptp  where  p  is  the 
dipole  moment  operator.  Then  corresponds  to  a  plane  wave  current  and 
j  is  the  current  when  a  plane  wave  is  scattered  by  a  perturber.  For 
l  ■  “  the  area  of  integration  corresponds  to  an  infinitesimal  solid  angle 
directed  along  k  and  displaced  infinitely  far  away.  This  yields  as  the 
average  of  the  total  scattering  cross  sections  for  the  ground  and  excited 
states.  For  1  =  0  velocity  changes  are  not  important  and  encloses  the 
perturber.  Then  oq  becomes  identical  with  the  phase  changing  cross  section 
as  obtained  by  Baranger .  ^ 

The  wide  variety  of  seemingly  unrelated  reordering  phenomena  which 
exists  today  has  prompted  us  to  develop  a  simple  diagrammatic  method  for 
their  analysis. With  this  method  one  treats  ordinary  two-pulse  photon 
echoes,  raman  echoes,  tri-level  echoes,  ^  grating  echoes,^  etc.  on 
the  same  basis. 

A  new  method  of  observing  coherent  transient  has  been  developed.  A 
probing  laser  pulse  is  used  to  modify  the  polarization  properties  and/or 
propagation  direction  of  a  coherent  transient  signal.  The  excitation 
pulses  can  then  be  blocked  by  spatial  filters  or  passive  polarization  elements, 
and  the  echo  can  be  observed  without  the  use  of  fast  optical  shutters.  This 
technique  should  prove  widely  applicable  and  especially  important  in  the 
study  of  picosecond  processes  where  sufficiently  fast  optical  shutters  are 
difficult  to  obtain.  Preliminary  experiments  have  been  performed  on  the 
resonance  lines  of  atomic  Li  and  T1  which  demonstrate  the  practicality  of  the 
technique. 
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Our  previous  measurements  of  stimulated  echo  relaxation  in  Na  vapors 

have  stimulated  us  to  develop  a  picosecond  capability  in  order  to  examine 

the  behavior  of  the  collision  kernel  at  short  times  where  our  data  base 

was  deficient.  We  also  will  use  this  apparatus  for  examining  two-pulse 

echo  relaxation  in  the  short  time  regime.  This  is  especially  important  in 

2 

the  case  of  Li  where  our  data  to  date  does  not  show  the  x  behavior  for  a 
which  is  to  be  expected  at  short  time.  As  of  this  writing  we  have  produced 
subnanosecond  pulses  by  using  an  actively  mode  locked  dye  laser  pumped  by 
a  cw  argon  laser.  A  pulsed  laser  is  being  used  to  pump  a  dye  laser  to 
selectively  amplify  one  of  the  picosecond  pulses. 

In  a  separate  experimental  setup  we  are  developing  a  pulsed  dye  laser- 
amplifier  system  using  a  YAG  laser  as  the  pump  source.  With  this  system 
we  plan  to  extend  our  measurement  of  Rydberg  atom  relaxation  properties 
into  regimes  previously  inaccessible  because  of  lack  of  excitation  power. 

The  first  experiment  we  plan  to  perform  will  be  to  extend  our  work  in  Na- 
vapor  using  the  tri-level  echo  technique  to  beyond  principle  quantum  numbers 
of  40.  We  expect  to  see  a  dramatic  Increase  in  collision  cross  section 
when  S  and  D  level  spacing  becomes  small  enough  so  that  £  and  n  changing 
collisions  become  probable. 

*This  research  was  also  supported  by  the  Office  of  Naval  Research  under 
Contract  N00014-78-C-0517 . 
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B.  RELAXATION  AND  EXCITATION  TRANSFER  OF  OPTICALLY  EXCITED  STATES  IN  SOLIDS* 


(K.  Chiang,  E.  Whittaker,  T.  J.  Chen,  S.  R.  Hartmann) 
(Principal  Investigator:  S.  R.  Hartmann  (212)  280-3272) 


3  3 

Measurements  of  photon  echo  modulation  on  the  P.-  H.  transition  of 

0  4 

3+ 

the  Pr  ion  in  a  crystal  of  LaF^  proved  to  be  singularly  successful  in 

3  3 

obtaining  spectroscopic  and  relaxation  data  for  all  aspects  of  the  P  -  H^ 

transition.  In  particular,  it  is  possible  to  obtain,  in  addition  to  the 

homogeneous  optical  linewidth,  a  complete  description  of  the  hyperfine 

3  3 

Hamiltonians  associated  with  both  the  PA  and  H,  states  as  well  as  the 

u  4 

associated  nuclear  transition  linewidths . ^  As  discussed  in  the  last 

progress  report,  there  are  experimental  difficulties  connected  with  carrying 

13  34- 

out  an  echo  modulation  experiment  on  the  H^  transition  of  Pr  in 

LaF^  and  the  data  which  we  presented  reflected  this  fact.  It  was  therefore 
decided  to  proceed  in  two  separate  ways.  We  would  improve  upon  our  experi¬ 
mental  technique  with  the  hope  that  we  could  make  the  technique  successful 

(2) 

in  this  application  and  we  would  separately  perform  a  PENDOR  experiment 
to  unambiguously  determine  the  ^Dj  excited  state  splittings  which  were  in 
doubt.  Both  approaches  were  successful. 

A  PENDOR  experiment  works  by  exploiting  the  sensitivity  of  the  three- 
pulse  stimulated  photon  echo  to  hyperfine  transitions  induced  by  a  resonant 
radiofrequency  magnetic  field.  The  resonance  condition  is  detected  by  a 
change  in  the  echo  signal  amplitude.  The  timing  sequence  is  shown  in  Fig.  1. 
The  first  and  second  90°  laser  excitation  pulses  were  produced  by  a  single 
N2  laser  pumped  dye  laser  whose  output  is  split,  one  half  going  into  the 
first  pulse  and  the  other  half  directed  to  a  white  cell  where  it  is  delayed 
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to  form  the  second  pulse.  The  third  90°  laser  excitation  pulse  is  produced 

by  a  separately  triggered  laser  pumped  dye  laser.  These  lasers  were  the 

same  as  were  described  in  the  last  progress  report. 

Stimulated  photon  echoes  were  produced  on  the  transition  connecting 

3  1 

the  lowest  levels  of  the  crystal-field  split  levels  of  the  H.  and 

4  z 

states.  These  levels  are  electronic  singlets  that  are  further  split  into 

three  doubly  degenerate  nuclear  levels  by  the  combined  influence  of  the 
3+ 

Pr  nuclear  quadrupolar  and  second  order  magnetic  dipolar  interactions. 

(3)  3  1 

Previously  measured  nuclear  splittings  for  both  the  and  states 

are  small  compared  to  the  10  Ghz  bandwidth  of  the  laser  excitation  pulses. 

Figure  2  shows  the  percentage  decrease  in  the  stimulated  photon  echo 

intensity  as  a  function  of  the  applied  RF  frequency.  The  dips  at  3.7  Mhz 

and  4.65  Mhz  correspond  to  the  energy  splittings  in  the  state,  and  agree 

(3) 

very  well  with  Erickson's  reported  values  of  3.7  Mhz  and  4.7  Mhz.  Our 

measured  linewidths  of  350  Khz  (FWHM)  are,  however,  larger  than  the  200  + 

50  Khz  linewidth  that  he  reported.  In  Fig.  2,  the  dotted  and  solid  lines 

are  Lorentzian  and  Gaussian  fits,  respectively,  to  the  data  points  with  a 

(FWHM)  linewidth  of  350  Khz.  The  quality  of  our  data  prevents  our  choosing 

one  over  the  other,  but  since  a  large  part  of  the  inhomogeneity  of  the 

19 

nuclear  hyperfine  levels  is  due  to  the  perturbing  effects  of  eleven  F 

neighbors,  one  would  expect  the  linewidth  to  be  Gaussian  rather  than  Lorentzian. 

(1)  3  3 

Fourier  analysis  of  the  H^-  Pq  transition  has  shown  that  the 

3 

inhomogeneous  linewidth  of  the  nuclear  transitions  in  the  Pq  state  is 

3 

approximately  20  Khz  while  the  inhomogeneous  linewidth  of  the  ground 

3 

state  is  approximately  200  Khz.  For  the  Pq  state,  J=0,  and  the  linewidths 

141 

are  due  primarily  to  the  static  and  quadrupolar  interactions  of  the  Pr 
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nucleus  with  the  neighboring  F  nuclei  and  the  electric  field  gradient  at 
141  3 

the  Pr  site  respectively.  For  the  state,  J=4,  and  the  close  spacing 

3 

of  the  crystal  field  split  levels  of  the  multiplet  gives  rise  to  an 

enhanced  moment  that  interacts  strongly  with  the  surrounding  nuclear  neighbors 
and  leads  to  an  increased  linewidth.  In  like  manner  the  levels  would  be 

3 

expected  to  be  broadened  beyond  the  20  Khz  linewidth  of  the  state. 
According  to  a  second  moment  calculation,  the  contribution  to  the  inhomo¬ 
geneous  linewidth  with  the  unenhanced  Pr-F  nuclear  dipolar  interaction  is 
13  Khz.  When  we  take  into  account  the  experimentally  determined  enhanced 
nuclear  gyromagnetic  ratio  *v  3  Khz/gauss  for  the  state^  we  obtain 
an  inhomogeneous  linewidth  of  approximately  40  Khz.  The  contribution  from 
other  sources  such  as  lifetime  broadening  ('v  600  hz)  and  phonon  processes 

(3) 

800  Hz)  are  negligible.  Since  we  typically  used  rf  pulses  20  secs 
long,  the  contribution  to  the  linewidth  due  to  the  finite  duration  of  the 
rf  field  is  8  Khz.  The  rather  large  measured  linewidth  that  must  be 
accounted  for  is  possibly  due  to  rf  power  broadening.  These  results  will 
be  published  shortly. 

We  have  also  taken  several  steps  to  improve  the  results  of  the  echo 

13  3+ 

modulation  experiment  on  the  at  the  LaF^:Pr  system.  The  first 

step  was  to  install  intracavity  etalons  in  the  Hansch  type  dye  lasers, 

thereby  reducing  laser  bandwidth  to  4  Ghz,  but  increasing  the  power  per 

unit  bandwidth  somewhat.  This  step  enhanced  our  signal-to-noise  by  a 

factor  of  three.  We  also  obtained  crystal  samples  whose  length  was  three 

times  longer  than  our  original  samples.  Since  echo  intensity  is  proportional 
2 

to  (aL)  ,  where  a  is  the  small  signal  absorptivity  and  L  the  interaction 


length,  this  step  enhanced  the  signal  by  about  a  factor  of  ten.  Despite 


these  improvements  we  found  that  the  overall  echo  envelope  decayed  to 

less  than  the  detection  noise  level  beyond  a  pulse  separation  of  one 

3+ 

microsecond,  for  a  sample  doped  with  1%  Pr  .  However,  a  dramatic  improve- 

3+ 

ment  was  noted  when  a  Pr  sample  with  lower  concentration  was  used.  At  0.1% 
3+ 

Pr  concentration,  a  measurable  signal  persists  out  to  pulse  separations 
of  nine  microseconds,  when  the  sample  is  maintained  at  a  temperature  of 
2.5K. 

We  have  also  made  a  number  of  improvements  in  the  computer  assisted 
data  taking  process,  which,  while  not  directly  improving  signal  size,  do 
effectively  improve  the  quality  of  our  measurements  by  increasing  data 
taking  efficiency.  In  particular,  two  new  computer  operated  peripherals 
have  been  installed  in  the  system.  The  first  is  a  computer  controllable 
electrical  attenuator  which  is  used  to  maintain  the  modulated  signal  within 
the  linear  range  of  the  data  acquisition  electronics.  The  computer  is  thus 
able  to  adjust  the  signal  level  automatically.  The  second  device  is  a  pair 
of  digitally  programmable  pulse  generators,  used  to  set  the  laser  pulse 
separation  and  Pockel's  cell  gate.  These  devices  are  of  a  "random  access" 
nature  allowing  the  pulse  separation  to  be  varied  over  several  orders  of 
magnitude  between  a  single  laser  shot.  This  allows  us  to  periodically 
measure  the  echo  magnitude  at  a  single  fixed  reference  time  without  suffering 
a  loss  in  the  data  taking  rate. 

The  reference  measurement  allows  us  to  compensate  for  echo  amplitude 
drift  during  the  course  of  a  single  data  taking  run,  as  well  as  providing 
a  normalization  factor  which  facilitates  averaging  data  sets  from  different 
runs.  The  result  of  all  these  improvements  can  be  seen  in  the  data  presented 
in  Fig.  3,  which  shows  the  echo  envelope  versus  pulse  separation  for  the 
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Figure  3:  Photon  echo  amplitude  versus  pulse  separation  for 
the  transition  at  Pr3+:LaF_. 
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13  3+ 

D2_  transition  of  .01%  Pr  rLaF^.  Each  point  is  the  average  of  100  or 
more  echoes  from  data  sets  taken  on  different  days. 

The  high  quality  of  this  data  allows  us  to  extract  a  significant  amount 

of  spectroscopic  information  about  the  system  under  investigation.  Fig.  4 

shows  the  discrete  Fourier  transform^  at  the  data  shown  in  Fig.  3.  Clearly 

3 

shown  are  the  resonances  due  to  the  hyper fine  structure  in  the  ground  and 
1  3+ 

excited  states  of  Pr  in  LaF^.  These  resonances  are  in  close  agreement 
with  the  values  deduced  in  previous  experiments ^  ^  including  our  PENDOR 
experiment  outlined  above.  In  addition  to  the  resonance  frequencies  we 
are  also  able  to  deduce  the  inhomogeneous  linewidths  of  the  nuclear  transi¬ 
tions.  The  values  for  the  ground  state  result  are  in  close  agreement  with 
our  previous  measurement.  The  deduced  linewidths  of  <  100  Khz  (our  present 
resolution  limit  due  to  finite  span  of  time  domain  data)  is  a  new  upper 
limit  on  the  excited  state  linewidths,  and  are  consistent  with  recently 
published  theoretical  discussions  of  the  inhomogeneous  nuclear  linewidths. 

Development  work  over  the  last  year  has  given  us  a  very  versatile 
apparatus  for  investigating  photon  echo  phenomena  in  solids.  We  intend  to 
pursue  a  vigorous  program  of  study  with  this  setup.  In  particular,  we 
intend  to  carry  out  modulation  experiments  similar  to  the  one  described 


(6) 


above  in  a  variety  of  other  ion-host  systems,  most  likely  beginning  with 

Er3+:LaF^.  In  addition,  we  hope  to  utilize  the  excited  state^  and  tri- 

level  echo'  '  techniques  for  the  first  time  in  solids.  A  particularly 

interesting  set  of  levels  for  use  in  multilevel  experiments  would  be  the 

3Pq  and  3Sq  levels  in  Pr^tLaF^.  Magnetic  interactions  between  Pr3+  and 

neighboring  nuclei  have  been  shown^^  to  be  the  dominant  line-broadening 

3  1 

mechanism.  These  effects  should  be  very  small  for  an  ion  in  a  Pq- 
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Figure  4:  Discrete  Fourier  Transform  of  data  shown  in  Fig.  3. 

The  peaks  at  3.72  MHz  and  4.79  MHz  are  the  excited 
state  hyperfine  resonances.  The  peak  at  8.51  MHz 
is  a  sum  frequency  of  the  excited  state  resonances 
and  obscures  a  ground  state  resonance  at  8.48  MHz. 
Also  present  are  ground  state  resonances  at  16.7  MHz 
and  25.0  MHz. 
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superposition  as  both  states  are  electronic  singlets,  having  only  nuclear 
quadrupole  interactions.  Use  of  the  ^Sq  state  is  an  exciting  prospect,  too, 
because  of  its  close  proximity  to  4fSd  configuration  states,  which  may 
induce  previously  unseen  effects  on  the  photon  echoes.  These  experiments 
should  be  under  way  in  the  near  future. 

*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  NSF-DMR80-06966. 
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SIGNIFICANT  ACCOMPLISHMENTS  AND  TECHNOLOGY  TRANSITION  REPORT 

Significant  Accomplishments 

1.  Solid  State  Electronics  and  Device  Fabrication 

We  have  fabricated  Schottky  barrier  diodes  with  a  silicon-aluminum- 
silicon  structure  by  a  laser  irradiation  process.  This  novel  technique  pro¬ 
duces  diodes  with  a  large  barrier  height  (>  0.95  eV)  and  good  rectifying 
properties,  and  it  is  likely  to  be  useful  in  multi-layer  integrated-circuit 
applications . 

By  making  use  of  the  large  diffusity  of  aluminum  in  polysilicon,  we  have 
found  a  low  temperature  process  for  making  a  poly  Si-Si  diode.  With  our 
approach,  both  single-crystal  and  poly crystalline  p-n  junctions  can  be  fab¬ 
ricated  on  the  same  IC  chip  giving  added  flexibility  in  the  design  of  in¬ 
tegrated  structures. 

2.  Physical  Properties  of  Solid  State  Devices 

We  have  recently  developed  a  theoretical  model  for  the  statistical 
fluctuations  of  dopant  impurities  in  ion-implanted  bipolar  transistor  struc¬ 
tures  and  the  minimum  device  dimensions  for  VLSI  system  reliability. 

3.  Energy  Transfer,  Photophysics  and  Photochemistry 

The  evolution  of  spin  polarization  due  to  collisions  between  ground- 
2 

state  atoms  in  ^1/2  states  is  dominated  by  spin  exchange.  Recently  we  have 
shown  that  collisions  between  alkali  atoms  also  destroy  the  electronic  spin 
at  a  rate  which  is  about  IX  of  the  spin  exchange  rate.  This  spin  destruc¬ 
tion  is  probably  due  to  magnetic  dipole  spin-spin  interactions,  and  it  has 
important  consequences  for  various  optically  pumped  devices. 
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While  investigating  spin  exchange  collisions  in  dense,  laser-pumped 
alkali  vapors  we  noticed  that  a  pronounced  narrowing  of  the  magnetic  resonance 
linewidth  occurred  when  the  laser  beam  intensity  increased.  In  all  other 
experiments  known  to  us,  an  increase  in  the  optical  intensity  causes  the 
magnetic  resonance  linewidth  to  broaden.  More  careful  experimental  and 
theoretical  studies  of  this  light  narrowing  phenomenon  showed  that  it  was 
due  to  the  interplay  of  the  optical  pumping  with  the  peculiar  properties 
of  spin  exchange  collisions.  The  phenomenon  may  be  of  importance  in  some 
optically  pumped  devices. 

By  combining  picosecond  transient  absorption  and  emission  methods  we 
have  succeeded  in  developing  a  model  for  the  interplay  between  solvent  po¬ 
larity  and  geometric  structure  on  excited  state  electron  transfer  processes. 
We  have  found  that  in  low  dielectric  solvents  the  electron  transfer  between 
two  organic  molecules  separated  by  three  methylene  groups  can  only  occur  for 
very  limited  configurations.  However  in  the  high  dielectric,  solvent  rapid 
electron  transfer  to  ion  radical  states  occurs  indicating  that  geometric 
factors  are  no  longer  crucial.  We  have  used  these  studies  to  establish  that 
the  dynamics  of  short  chain  molecules,  or  polymer  segments,  can  be  obtained 
from  measurements  of  these  electron  transfer  processes. 


Technology  Transition  Reports 

1.  Multifrequency  heterodyne  techniques  developed  by  M.  C.  Teich  under 

JSEP  support  (Research  Unit  12)  are  currently  under  consideration  for  possible 
application  in  over-water  pilot  rescue  missions  where  the  use  of  a  light¬ 
chopping  mechanism  is  interdicted  because  of  its  interruptive  nature 
(Naval  Weapons  Center) . 

2.  The  Neyman  Type-A  characterization  of  photomultiplier  fluorescence 


noise  developed  by  M.  C.  Teich  under  JSEP  support  (Research  Unit  13)  is 
being  used  to  design  the  star  scanner  for  the  optical  guidance  system  in 
the  NASA/JPL  Galileo  orbiter  mission  scheduled  to  explore  Jupiter  in  1988. 
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